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INTRODUCTION 

This  manual  describes  some  standard,  general  methods  to  analyze 
strong-motion  accelerograph  records  of  highway  bridge  motions.  These 
analysis  methods  are  appropriate  for  typical  highway  grade  separation 
bridges.  These  bridges  may  be  steel  or  concrete,  simple  spans  or 
continuous,  straight  or  curved,  but  are  restricted  to  being  moderate  in 
length.  The  assumption  that  all  bridge  support  motions  are  identical 
requires  that  bridge  lengths  be  moderate. 

Methods  are  described  for  the  analysis  of  analog  records  visually 
without  the  need  of  a  high-speed  digital  computer  and  for  the  analysis  of 
digitized  data  with  the  use  of  a  computer.  The  objectives  of  both  types  of 
analysis  are  (1)  to  find  force  levels  throughout  the  bridge  for  a  particular 
earthquake,  and  (2)  to  find  parameters  of  a  mathematical  model  which  best 
describes  the  measured  response.  There  are  many  secondary  objectives  that 
may  be  considered  also,  but  for  brevity  only  these  two  primary  objectives 
are  considered  here. 

The  analyst  is  assumed  to  be  familiar  with  structural  engineering  and 
to  have  some  knowledge  of  structural  dynamics.  Appendix  B  includes  a 
section  on  structural  dynamics  in  which  those  concepts  that  are   needed  for 
the  specific  analysis  procedures  outlined  here  are  developed.  A  greater 
understanding  of  structural  dynamics  may  be  obtained  from  standard  texts 
(Biggs,  1964;  Clough  and  Penzien,  1975;  Hurty  and  Rubinstein,  1964), 

A  few  words  of  caution  are  in  order  about  the  use  of  this  manual. 
There  are  few  uncertainties  in  obtaining  force  levels  in  bridges  from 
strong-motion  records.  Force  levels  are  computed  from  an  extension  of  the 
expression 

(Force)  =  (Mass) (Acceleration) 


Masses  can  be  estimated  quite  accurately  from  drawings,  and  accelerations 
can  be  measured  accurately  at  specific  points.  There  are  always 
uncertainties  in  describing  complete  bridge  accelerations  from  the 
accelerations  of  a  few  points,  but  with  the  judicious  selection  of 
instrument  locations  these  uncertainties  can  be  minimized.  Therefore,  force 
level  analyses  can  be  performed  with  confidence. 

On  the  other  hand,  when  strong-motion  records  are  analyzed  to  obtain 
model  properties,  many  uncertainties  exist.  In  this  manual  properties  are 
sought  for  only  a  linear  mathematical  model.  This  model  is  chosen  because 
it  is  commonly  used  in  structural  dynamics,  because  general  solution  methods 
exist  to  solve  the  resulting  set  of  equations,  and  because  that  solution  is 
relatively  easy  to  obtain.  This  model  is  not  selected  because  it  happens  to 
be  the  best,  although  fortunately  it  is  quite  good  for  non-damaging  levels 
of  motion.  It  may  not  be  appropriate  for  high-level  damaging  motions  or  for 
motions  involving  significant  foundation  compliance. 

It  is  beyond  the  scope  of  this  manual  to  discuss  higher  order 
mathematical  bridge  models,  which  are  appropriate  for  nonlinear  or  inelastic 
response.  It  would  also  be  \/ery   difficult  since  there  is  no  single,  general 
higher  order  model;  there  are  many.  This  whole  field  of  nonlinear  higher 
order,  model  identification  from  strong-motion  records  is  in  an  early  stage 
of  development.  Even  the  methods  suggested  here  to  identify  model 
properties  for  the  linear  model  are  not  well  established  at  this  time. 
Rapid  advancement  is  expected  in  the  next  few  years,  particularly  with 
regard  to  the  automatic  analysis  of  digitized  records  on  high  speed  digital 
computers. 


I.  STRONG-MOTION  INSTRUKENTATION  AND  RECORDS 

A  strong-motion  record  is  defined  here  as  the  record  of  structural 
motion  caused  by  a  relatively  large  earthquake.  This  strong  motion  may  be 
recorded  in  a  number  of  ways.  It  may  be  a  light  trace  on  film,  a  line  on  a 
strip  of  paper,  or  an  analog  or  digital  signal  on  magnetic  tape.  In  all 
cases  the  motion  is  recorded  as  a  function  of  time. 

The  motion  recorded  is  that  for  a  point  moving  along  a  single  axis  and 
may  have  dimensions  of  displacement,  velocity,  or  acceleration.  Most  often 
acceleration  is  recorded.  The  motion  is  recorded  as  an  electrically  or 
mechanically  generated  signal  which  is  proportional  to  the  motion  within 
some  acceptable  error. 

The  adjective  "strong"  is  included  to  differentiate  these  recordings 
from  recordings  taken  with  other  common  seismographs  used  to  determine 
earthquake  magnitudes  and  epicentral  distances.  Such  seismographs  are 
capable  of  detecting  earthquakes  almost  anywhere  in  the  world,  and  are 
necessarily  very   sensitive.  Strong  motions  would  cause  these  seismograph 
recordings  to  be  off  scale.  In  general,  strong  motions  could  cause  damage 
to  structures,  and  are  recorded  near  earthquake  epicenters.  The  actual 
level  of  ground  motion  above  which  motions  are  "strong"  is  not  well  defined. 
Typically,  ground  motions  with  a  peak  acceleration  greater  than  0.01  ^ 
trigger  strong-motion  accelerometers,  those  greater  than  0.05  £  are  of 
structural  interest,  and  those  greater  than  0.10  £  are  potentially  damaging. 

In  the  United  States,  strong-motion  records  are  normally  recorded  on 
accelerographs.  The  accelerographs  are  of  two  general  types:  analog  or 
digital.  In  either  case,  acceleration  is  sensed  by  accelerometers  and 
recorded  as  a  continuous  function  of  time.  In  the  analog  instrument,  it  is 
recorded  on  light-sensitive  film  or  paper  and  in  the  digital  instrument,  it 


is  recorded  in  digital  form  on  magnetic  tape.  Both  instrument  types  are 
powered  by  batteries,  are  triggered  into  operation  by  the  strong  motion 
itself  (i.e.,  motion  having  an  acceleration  equal  to  or  greater  than  a 
pre-set  triggering  level,  commonly  0.01  £,  and  are  available  either  as 
self-contained  triaxial  instruments  or  remote  accelerometer/central  recording 
instruments.   In  the  self-contained  variety,  all  major  components--accelero- 
meters,  recorder,  and  trigger— are  housed  in  one  container,  whereas  in  the 
remote  recording  type,  these  elements  can  be  physically  separated  but  are 
interconnected  by  low-voltage  data  cable. 

A  typical  analog  film  record  from  a  three-component  self-contained 
accelerograph  is  shown  in  figure  1  (records  from  remote  accelerometer/central 
recording  systems  are  similar  although  they  normally  contain  more  data 
traces).  The  record  contains  time  ticks,  two  fixed  traces,  three  accelera- 
tion data  traces  and  a  WIJVB  time  code.  The  time  ticks  denote  film  speed 
and  normally  appear  every   half  second;  each  fixed  trace  is  a  reference 


-  .   .    .     .        ^WWVB  Time  Code 
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Figure  1.-  Copy  of  a  typical  three-component  strong-notion  accelerograph  record, 


trace  that  is  used  to  extract  long-period  errors  from  the  data  due  to  film 
distortion  or  the  film  drive  mechanism;  and  the  WWVB  time  code  is  used  to 
identify  when  the  event  occurred.  The  most  important  elements  of  the  record, 
the  acceleration  data  traces,  have  values  of  amplitude  that  are  approximately 
proportional  to  real  values  of  recorded  acceleration.  For  most  instruments 
in  present  use,  the  sensitivity  of  the  acceleration  traces  is  on  the  order 
of  1.90  cm/g.  The  sensitivity  is  calibrated  during  or  before  installation 
and  should  appear  on  the  record  or  accompany  it.  Other  relevant  information 
that  should  accompany  the  record  includes  orientation,  location,  natural 
frequency,  damping  of  each  accelerometer,  film  speed,  instrument  serial 
number,  and  location  of  the  recorded  earthquake.  If  such  information  does 
not  accompany  a  record,  most  of  it  should  be  available  either  from  the 
instrument  owner,  the  manufacturer,  or  the  agency  or  organization  maintaining 
the  instrument. 

Records  from  a  typical  digital  instrument  are  normally  recorded  in 
digital  form  on  one-  or  four-track  cassette  tape  and  are  not  immediately 
available  in  analog  form.  Analog  records  are  made  from  a  digital  instrument 
by  processing  the  digital  data  with  a  digital-to-analog  converter  and 
plotter.  Hardware  of  this  type  has  been  developed  by  the  firms  that 
manufacture  the  digital  instruments  and  should  be  available  either  through 
these  firms  or  the  agency  or  organization  that  maintains  the  instrument. 

In  order  to  convert  both  the  analog  film  and  digital  cassette  records 
into  a  digital  form  that  can  be  processed  on  high-speed  digital  computers, 
additional  processing  is  required.  Analog  records  are  converted  into 
digitized  form  on  systems  known  as  digitizers.  In  such  a  system,  a  handheld 
cursor  or  automatic  scanning  device  is  used  to  follow  the  trace  and  assign 
numerical  coordinate  values  (normally  Cartesian)  to  all  points  of  interest. 


Such  points  are  norraally  selected  at  equal  Intervals  or  at  all  peaks, 
valleys  and  points  of  Inflection.  The  nwierical  values  are  recorded  on 
computer-compatible  cards  or  on  iMgnetlc  tape.  Data  frcMis  a  digital 
strong-ntotlon  Instrument  is  Immediately  available  in  digitized  form  on 
cassettes  but  requires  processing  for  conversion  into  a  computer-compatible 
format.  The  processing  is  handled  by  hardware  designed  by  the  instrument 
manufacturer.  Its  availability  should  be  known  by  the  strong-motion 
instrument  manufacturer  or  the  agency  or  organization  maintaining  the 
instrument. 

After  conversion  Into  a  high-speed-coraputer-compatible  format,  the 
digitized  data  can  be  corrected  to  extract  errors  Introduced  during  the 
recording  and  digitization  process  and  to  account  for  the  dynamic 
characteristics  of  the  recording  Instrument.  Procedures  for  correcting  the 
data  are  vrell  established  (Brady  and  Perez,  1976;  ?11klofsky  and  Mancini, 
1977;  Trifunac  and  Lee,  1973)  and  are  presently  performed  routinely  by  the 
U.S.  Geological  Survey.  When  the  processing  is  completed,  the  digitized 
data  are  then  ready  for  analysis.  Procedures  for  standardized  analysis  of 
three-component  ground  motion  records  are  also  njell  estdblished  (Brady  and 
Perez,  1976;  Trifunac  and  Lee,  1973)  and  are  routinely  performed.  The 
analysis  of  a  set  of  earthquake  records  from  a  particular  bridge,  however, 
is  not  yet  a  routine  procedure.  The  methodology  described  below  has  been 
designed  for  that  purpose. 


II.  FORHULATIOfI  OF  A  HATHEMATICAL  mOEl 
FOR  THE  DYNAHIC  BEHAVIOR  OF  BRIDGES 

As  indicated  above,  the  two  primary  objectives  in  the  analysis  of 
strong-motion  records  are:  (1)  to  verify,  modify,  or  improve  the 
mathematical  models  used  to  design  a  bridge  to  resist  earthquake  motion  and 
(2)  to  find  force  levels  throughout  a  bridge.  If  a  number  of  analyses  lead 
to  a  consistent  set  of  modeling  assumptions,  then  this  set  of  assunq)tions 
can  be  used  vfith  confidence  in  the  design  of  new  bridges.  It  is  hoped  that 
the  analyses  verify  rather  than  modify  the  initial  design  assumptions. 

A  mathematical  model  of  the  bridge  is  formulated  from  the  initial 
design  of  the  bridge.  Nianerous  assumptions  are  made  in  the  formulation  of 
this  model  and  these  assumptions  are  to  be  verified  from  the  measured 
records.  The  bridge  response  motions,  or  output  motions,  are  computed  using 
the  mathematical  bridge  model  and  a  given  design  earthquake.  From  these 
output  motions  force  and  stress  levels  at  selected  points  in  the  bridge  are 
computed.  These  force  and  stress  levels  are  then  compared  with  allowable 
levels.  If  the  computed  levels  exceed  the  allowable  levels,  then  the  bridge 
design  must  be  modified  and  reanalyzed  after  a  new  mathematical  model  has 
been  formulated. 

Experimentally  it  is  desirable  to  measure  the  input  ground  motions  and 
a  sufficient  number  of  bridge  response  motions  to  describe  the  complete 
bridge  motion  to  within  some  acceptable  error.  From  the  output  motions 
alone,  the  force  and  stress  levels  throughout  the  bridge  may  be  computed  and 
compared  with  corresponding  acceptable  levels.  The  force  levels  for  a 
specific  earthquake  are  not  of  general  interest  because  that  particular 


earthquake  or  any  scale  multiple  of  it  is  never  likely  to  repeat  itself.  Of 
far  greater  general  interest  is  to  identify  the  mathematical  model  that 
would  have  duplicated  the  measured  output  motions  for  the  measured  input 
motions.  A  comparison  of  this  model  and  the  model  used  to  design  the  bridge 
will  test  the  assumptions  made  to  formulate  the  model  for  design. 

In  general,  mathematical  models  of  dynamic  bridge  behavior  are 
formulated  during  the  design  process  to  predict  what  distortions  are 
expected  relative  to  the  ground  when  a  bridge  is  subjected  to  an  earthquake. 
It  is  important  to  recognize  that  these  distortions  are  in  addition  to  the 
absolute  displacements  of  the  bridge  following  the  ground  motions.  If  there 
are  no  distortions  (the  bridge  is  essentially  rigid)  there  will  still  be 
significant  earthquake  forces  in  the  bridge  equal  to  the  products  of  bridge 
masses  and  ground  accelerations.  We  hope  to  learn  from  the  dynamic  model  of 
the  bridge  whether  or  not  these  distortions  are  significant  relative  to  the 
absolute  ground  displacements.  In  some  cases  the  distortions  are  \/ery   small 
relative  to  the  absolute  displacements  of  the  bridge.  For  such  cases,  the 
accuracy  of  the  dynamic  model  is  unimportant  because  the  contribution  to  the 
total  inertial  load  from  the  flexible  distortions  is  so  small.  For  certain 
loads,  the  bridge  may  "resonate"  producing  yery   large  distortions  from  a 
relatively  small  input.  Such  vibrations  rarely  result  from  earthquakes.  For 
\/ery   flexible  structures,  distortions  one  to  four  times  the  absolute  ground 
displacements  may  be  expected  (Newmark  and  Rosenblueth,  1971;  Wiegel,  1970). 
For  relatively  rigid  bridges,  distortions  much  less  than  the  absolute  ground 
motion  are  likely. 

The  importance  of  the  mathematical  model  of  dynamic  behavior  depends 
upon  the  bridge  being  analyzed.  Therefore,  there  is  no  absolute  measure  of 


the  accuracy  of  mathematical  models  for  all  bridges,  but  it  must  be 
evaluated  for  each  class  of  bridges. 

A  few  fundamental  assumptions,  generally  made  in  the  formulation  of  a 
dynamic  mathematical  model,  will  be  made  in  this  report.  It  is  first 
assumed  that  a  linear  elastic  model  may  describe  adequately  the  dynamic 
distortions  of  the  bridge.  Forces  are  assumed  to  be  linearly  related  to 
bridge  displacements  and  velocities  relative  to  the  ground,  and  linearly 
related  to  absolute  bridge  accelerations.  The  second  fundamental 
assumption  is  that  all  supports  of  the  bridge  are  subjected  to  the  same 
ground  motion  at  the  same  time.  For  this  assumption  to  hold,  the  bridge 
length  must  be  small  with  respect  to  the  wave  lengths  of  the  ground  motion 
waves  which  produce  the  most  significant  bridge  response,  presumably  at  the 
lower  natural  frequencies  of  the  bridge.  This  assumption  seems  reasonable 
for  typical  highway  bridges  on  firm  foundation  material,  but  may  be 
erroneous  if  the  foundation  material  varies  significantly  from  one  end  of 
the  bridge  to  the  other. 

Both  of  these  assumptions  are  made  not  necessarily  because  they  will 
lead  to  the  most  accurate  model,  but  because  they  will  lead  to  mathematical 
equations  that  can  be  solved  with  comparative  ease.  Once  these  assumptions 
are  made  the  mathematical  model  is  well  defined  and  the  model  parameters  may 
be  determined  from  the  response  of  the  structure.  If  the  model  is  not 
assumed  to  be  linear,  the  problem  is  not  well  defined  and  there  is  no  end  to 
the  possible  forms  that  the  model  may  take. 

Ultimately,  these  two  modeling  assumptions  must  be  evaluated  by 
analyzing  actual  strong-motion  records.  If  model  parameters  can  be  found 
with  accuracy  such  that  a  mathematical  model,  based  upon  the  above  two 


assumptions,  can  duplicate  the  measured  response  well,  then  it  can  be 
deduced  that  the  two  assumptions  are  reasonable.  If  model  parameters  cannot 
be  found,  then  the  assumptions  are  poor.  It  is  beyond  the  scope  of  this 
report  to  suggest  higher  order  models  that  may  be  evaluated.  For 
information  on  higher  order  models  the  reader  should  consult  Clough  and 
Penzien  (1975)  and  Hurty  and  Rubinstein  (1964). 

Within  the  context  of  these  two  assumptions,  the  mathematical  model 
will  take  one  of  three  forms  depending  upon  the  coordinate  system  chosen  to 
describe  the  bridge  distortions. 

The  first  and  most  common  form  of  mathematical  model  is 

[M]{u}+[C]{u}+[K]{u}  =  -{M^>%-^VV^^z^'^zg  ^^^ 

where 

{u}  is  a  vector  of  bridge  displacements  relative  to  the 

ground  at  a  finite  number  of  preselected  points, 

[M]  is  a  matrix  of  masses  lumped  at,  and  corresponding 

to,  the  physical  coordinates  chosen  to  describe  the 
bridge  distortions, 

[C]         is  a  damping  matrix, 

[K]  is  a  stiffness  matrix, 

^xa'  ^vn'  ^za      ^^^   mutually  orthogonal  ground  accelerations, 

(•)=d/dt      is  differentiation  with  respect  to  time,  and 

{M  }j  {M.,}*  Wy}  are  vectors  of  participation  factors  corresponding 

x     y     z 

to  the  ground  motions  U  ,  U  ,  U  ,  respectively. 

xg       yg       zy 

The  participation  factor  M.   is  the  quasi -static  force  in  the  ith  coordinate 

direction  produced  by  a  unit  acceleration,  \}^(t)=].     Similar  definitions 

xy 


hoi 


d  for  the  participation  factors  M  •  and  M  . 
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This  mathematical   model   is  called  a  lumped-mass  model    in  physical 
coordinates.     It  can  be  used  to  identify  mathematical  models  for  short 
simple-span  bridges  and  yery  long  multi-span  continuous  bridges.     It  is 
called  a  lumped-mass  model   because  the  mass  of  the  bridge  is  typically 
lumped  (for  analysis  purposes)  at  a  selected  number  of  points.     The 
coordinates  are  physical   because  they  are  translations  and  rotations  at 
selected  points  which  are  measureable. 

For  many  bridges,  particularly  short  continuous  bridges,  it  may  be  more 
appropriate  to  describe  the  bridge  distortions  in  generalized  coordinates. 
This  is  the  second  form  of  mathematical   model.     In  this  coordinate  system, 
bridge  distortions  are  expanded  in  terms  of  a  set  of  r  assumed  distorted 
shapes. 

u(s,t)   =     2  $  (s)a.(t)  (2) 

i=l   ^         ^ 

where 

u(s,t)  is  a  vector  of  bridge  distortions  as  a  function  of 

time,  t,  and  s,  the  distance  along  the  centerline  of 
the  bridge  element  (bridges  typically  can  be  modeled 
as  a  collection  of  one-dimensional   elements), 

$.(s)  is  a  vector  description  of  the  ith  assumed  distorted 

shape  and 

a-(t)  is  a  time-dependent,  generalized  displacement 

associated  with  the  ith  assumed  distorted  shape. 

At  any  time  the  vector  of  distortions  at  a  location  s  is  given  by 

u(s,t)  =  u^(s,t)i+u^(s,t)j+u^(s,t)k  (3) 

where  T,  j,  and  k  are  unit  vectors  in  the  mutually  perpendicular  directions 
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X,  y,  z  respectively  parallel  to  the  recorded  ground  accelerations  d^At), 
U  (t),  U  (t).  The  vector  description  of  the  assumed  distorted  shapes, 
ii).(s),  similarly  can  be  written  in  terms  of  its  orthogonal  components  in  the 
X,  y,  and  z  directions.  Note  that  the  generalized  displacements,  a. (t), 
cannot  be  measured  directly,  but  must  be  computed  from  displacements  in 
physical  coordinates,  u-(t),  by  the  inverse  of  the  transformation 

{u(t)}  =  [i|^]{a(t)}  (4) 

The  matrix  element  i|^.  •  is  the  displacement  in  the  ith  physical   coordinate 
direction  for  the  jth  assumed  distorted  shape. 

The  mathematical  model   describing  the  dynamic  bridge  distortions  can 
now  be  given  by 

[m]{a}+[c]{a}=[k]{a}  =  -{mx>Uxg-{my}Uyg-{mpU^g  (5) 

The  inertial   terms,  integrated  over  all  elements  of  the  bridge,  are 

m^.j  =  /M(s)5.(s).5j(s)ds  (6) 

m^^.  =  /M(s)5.(s).ids  (7) 

niyi  =  /M(s)5i(s).jds  (8) 

f"zi   "  /Ks)ip^-(s).kds  (9) 

and  typically  are  computed  from  material   properties,  the  geometry  of  the 
bridge  being  analyzed,  and  the  assumed  deflected  shapes.     The  quantity  M(s) 
is  an  inertia  per  unit  length  of  bridge  (a  mass  when  $(s)  is  a  displacement 
and  a  rotational   inertia  when  J(s)   is  a  rotation).     Analytical  expressions 
for  the  damping  coefficients,  c-.  and  stiffness  coefficients  k-.,  can  also 
be  formulated.     However,  these  are  the  unknowns  to  be  obtained  from  the 
strong-motion  records. 

Although  this  formulation,  in  terms  of  generalized  coordinates,  a-(t). 
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looks  identical  to  the  model  formulation  in  terms  of  physical  coordinates, 
u-(t),  there  is  one  significant  advantage  for  most  short  continuous  bridges. 
To  achieve  the  same  degree  of  accuracy  describing  structure  distortions,  it 
will  require  fewer  coordinates  a. (t)  than  u-(t).  Since  one  transducer  is 
needed  to  record  each  distortion,  in  either  coordinate  system,  the  system 
with  the  fewer  coordinates  required  will  require  fewer  transducers.  Another 
way  of  interpreting  this  advantage  is  that  the  coordinate  system  requiring 
the  fewest  number  of  coordinates  to  describe  the  bridge  distortions  for  a 
given  accuracy  will  describe  the  bridge  distortions  most  accurately  for  the 
same  number  of  transducers. 

For  all  bridges,  distortions  can  be  described  most  accurately,  with  the 
fewest  number  of  coordinates,  if  normal  coordinates  are  used  (Clough  and 
Penzien,  1975),  This  is  the  third  form  of  the  mathematical  model.  Here 
bridge  distortions  are  expanded  in  terms  of  normal  mode  shapes 

u(s.t)  =  i^i5^.(s)n^(t)  (10) 

where 

$^-(s)        is  the  ith  normal  (or  natural)  mode  shape,  and 

T]At)  is  a  displacement  in  normal  coordinates  associated 

with  the  ith  normal  mode  shape. 
In  normal  coordinates,  the  matrix  equation  of  motion  describing  structure 
distortions  takes  the  unique  form 

[M]{n}+[C]{n}+[K]{n}=-{M^}U^g-{My}U^^g-{M^}U^g  (11) 

where  all  off-diagonal  terns  of  [M],   [C],  and  [/(]  equal  zero.     The  inertial 
terms,  integrated  over  elements  of  the  bridge,  are  again  given  by 

M^  =  /M(s)$.(s).$.(s)ds  (12) 

M^i  =  /M(s)5.(s).ids  (13) 
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My,-  =  /H{s)t^(s)-jds  (14) 

M^,-  =  /M(s)5^.(s)-kds  (15) 

It  is  particularly  advantageous  to  describe  bridge  motions  in  terms  of 
normal  coordinates  because  the  individual  equations  in  the  matrix  equation, 
equation  11,  are  uncoupled.  Specifically,  each  normal  displacement  is 
governed  by  the  single  equation  of  motion  (here  for  example,  the  ith) 

M  fi.+C-n+Kn.  =  -M  -U  -M  -U  -M  -U  (16) 

•^1  'i  1  'i  1  'i     XI  xg  yi  yg  zi  zg  ^  ' 

Typically  these  equations  are  normalized  to  have  modal  masses,  M-,  equal 

unity  (Clough  and  Penzien,  1975).  For  such  a  normalization,  equation  16 

takes  the  form 

n.+2^-w-n  +a)-n-=-Y  -*U  -y  -U  -y  -U  (17) 

'ii  '-«-i'"i'i,  lu^'i,  'xi  xg  "yi  yg  ''zi  zg  ^  '' 

in  terms  of  the  common  dynamic  response  characteristics 

(0-  =  /K-/M-  a  natural  circular  frequency  of 

vibration  (radians/sec), 
^.  =  C./2/K|M.  a  modal  viscous  damping  ratio,  and 

Yy--,  Y„i»  Y,i  renormalized  participation  factors: 

^xi^**i'  %'^^i  ^"^  '^zi^'^i  '^spectively. 
In  this  form,  a  participation  factor  is  a  measure  of  the  degree  to  which  a 
mode  of  vibration  is  excited  by  a  ground  motion  relative  to  the  response  of 
a  simple  spring-mass  oscillator  having  the  same  natural   frequency  and 
damping  ratio. 

It  is  desirable  to  describe  the  distortions  of  all  bridges  in  terms  of 
normal  modes  because  the  fewest  coordinates  are  needed  for  a  given  accuracy. 
Similarly,  for  a  given  number  of  normal  coordinates  there  are  the  fewest 
unknown  model  parameters  to  be  determined  from  the  strong-motion  records 
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because  all  off-diagonal  terms  of  [M],  [C],  and  [K]  by  definition  equal 
zero.     Unfortunately,  when  model  parameters  are  to  be  determined  from 
strong-motion  records  as  they  are  in  this  report,  bridge  distortions  cannot 
be  described  in  normal  coordinates  before  analysis  because  the  normal  mode 
shapes  ^As)  are  unknowns  that  are  to  be  determined  from  the  analysis  of  the 
strong-motion  records.     These  mode  shapes  must  be  known  to  compute  the 
normal  displacements,  Tij(t).     Again  the  normal  displacements  n,-(t)  are  not 
directly  measureable,  but  must  be  computed  by  the  inverse  of  the 
transformation 

{u(t)}=[<{>]{Ti(t)}  (18) 

where  <j)- -  is  the  displacement  of  the  jth  normal  mode  shape  in  the  ith 
physical  coordinate  direction. 

Therefore,  for  the  methods  of  analysis  described  here,  bridge 
distortions  must  be  described  in  either  physical  or  generalized  coordinates. 
The  lumped-mass  model  in  physical  coordinates  should  be  used  for  simple  span 
or  long  multi-span  bridges,  whereas  the  model  in  ^neralized  coordinates  is 
more  appropriate  for  short  continuous  bridges.     Once  model  parameters  are 
identified  from  the  strong-motion  records,  the  identified  mathematical  model 
can  be  transformed  readily  from  physical  or  generalized  coordinates 
(equations  1  or  5)  to  the  decoupled  form  in  normal  coordinates  (equation 
17).     Methods  to  make  this  transformation  from  either  physical  or 
generalized  coordinates  to  normal  coordinates  can  be  found  in  any  text  on 
structural  dynamics  (Biggs;  1964,  Clough  and  Penzien,  1975;  Hurty  and 
Rubinstein,  1964). 

If  the  reader  is  interested  in  a  more  complete  development  of  the 
equations  presented  above.  Appendix  B  should  be  consulted  or  the  referenced 
texts  on  structural  dynamics.     Gere  and  Weaver,  1965,  give  a  complete 
development  of  matrix  algebra. 
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III.  ANALYSIS  OF  ANALOG  RECORDS 

A.  Introduction 

This  section  describes  methods  of  analyzing  analog  recordings  of 
strong-motion  records  from  bridges.     For  use  in  the  procedures  described, 
the  ground  motion  and  bridge  response  values  may  be  visually  scaled  from  the 
raw  recordings.     Also,  analysis  methods  described  in  this  section  may  be 
performed  without  a  high-speed  digital   computer. 

As  mentioned  in  Section  II  there  are  two  principal  objectives  of  the 
strong-motion  record  analysis  procedures.     The  first  is  to  determine  force 
levels  throughout  the  bridge  for  the  measured  earthquake.     The  second  is  to 
deduce  from  the  ground  motion  and  bridge  response  records  the  mathematical 
model  of  the  bridge  which  best  duplicates  the  bridge  response.     Mathematical 
models  are  most  often  characterized  by  natural   frequencies,  damping  ratios, 
and  mode  shapes.     In  order  to  find  force  levels  there  is  no  need  to 
formulate  a  mathematical  model   describing  the  bridge  motion  as  complete  and 
as  complex  as  those  formulated  in  the  previous  section  (equations  1,  5,  11, 
and  17).     Instead,  force  levels  can  be  computed  from  the  recorded  bridge 
motions,  known  geometry  of  the  bridge,  and  known  mass  distribution  in  the 
bridge.     The  recorded  bridge  motions,  however,  must  be  described 
mathematically  in  terms  of  one  of  the  three  coordinate  systems  described  in 
Section  II. 

B.  Force  Levels  in  Bridges  from  Strong-Motion  Records 

It  is  assumed  that  accelerations  are  measured  as  raw  data.     This  may 
seem  to  be  overly  restrictive,  but  in  practice  nearly  all   strong-motion 
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instruments  do  record  acceleration  as   raw  data. 

A  second  more  important  assumption  is  that  the  bridge  is  properly 
instrumented  in  the  first  place.     By  proper  instrumentation  it  is  meant  that 
there  exist  a  sufficient  number  of  time-correlated  bridge  accelerations  to 
describe  all   significant  bridge  motions.     This  is  absolutely  necessary  if 
accurate  force  levels  are  to  be  measured  from  the  experimental   data  without 
the  use  of  a  mathematical   model.     Many  instrumentation  programs  have  as 
their  principal   objective  the  verification  of  a  particular  mathematical 
model.     Relatively  few  strong-motion  instruments  may  be  required  to  verify  a 
mathematical   model   adequately,  but  these  may  be  totally  inadequate  for  the 
calculation  of  force  levels  in  the  bridge  without  the  use  of  a  mathematical 
model.     Therefore,  it  is  important  to  determine  first  whether  or  not  there 
are  enough  recordings  to  perform  a  force  analysis.     If  there  are  not,  a 
force  analysis  can  still   be  made  but  only  with  the  assistance  of  a 
mathematical   model.     The  mathematical  model   is  first  verified  by  the 
experimental   data,  and  then  force  levels  are  computed  from  the  mathematical 
model . 

For  a  typical   bridge  deck  it  is     assumed  that  all   significant 
horizontal    distortions  can  be  described  in  terms  of  four  modes  of  vibration: 
a  lateral  translation  relative  to  the  ground,  a  rotation  about  a  vertical 
axis,  a  lateral   flexural   mode,  and  a  longitudinal  translation  relative  to 
the  ground.     These  four  modes  of  vibration  and  the  corresponding,  assumed, 
displacement  shapes  relative  to  the  ground  are  shown  in  figure  2,     In  order 
that  force  levels  be  found  throughout  the  bridge,  bridge  motions  in  all  of 
these  modes  of  vibration  must  be  recorded.     In  general   there  must  exist  at 
least  one  recorded  motion  for  each  mode  of  vibration  of  interest.     These 
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recorded  motions,  however,  must  be  recordings  of  independent  motions  to 
describe  the  above  modes  of  vibration  sufficiently.  For  example,  the  four 
instrument  motions  shown  in  figure  3a  may  have  been  recorded  in  hopes  of 
detecting  the  modes  of  vibration.  The  motions  U-iCt)  and  U2(t)  can  be 
combined  to  yield  lateral  translation  and  rotation  but  cannot  be  combined  to 
yield  any  information  regarding  deck  flexure.  Furthermore,  either  U^(t)  or 
U^(t)  can  be  used  to  record  the  longitudinal  deck  motion.  Since 
longitudinal  deck  extension  is  likely  to  be  minimal,  U^Ct)  and  U«(t)  will  be 
nearly  identical  and  hence  one  is  redundant.  Therefore,  there  are  only 
three  independent  motions  recorded,  and  all  four  expected  modes  of  vibration 
cannot  be  detected  by  these  three  motions.  A  more  suitable  instrumentation 
plan  is  shown  in  figure  3b.  Corresponding  to  the  assumed  deflected  shapes 
ii;.(s)shown  in  figure  2,  generalized  accelerations  can  be  expressed  in  terms 
of  the  measured  absolute  accelerations  by 

a^Ct)  =  (U^(t)+U3(t))/2-U^g(t)  (19) 

a2(t)  =  (U3(t)-U^(t))/L  (20) 

a3(t)  =  (2U2(t)-U^(t)-U3(t))/2  (21) 

a4(t)  =  U4(t)-U^g(t)  (22) 

These  equations  are  equivalent  to  the  inverse  of  the  transformation 
from  physical  to  generalized  coordinates,  equation  4.  The  generalized 
accelerations,  a-(t)(i=l ,2,3,4)  can  now  be  used  to  determine  force  levels  in 
the  bridge  associated  with  the  four  modes  of  vibration  shown  in  figure  2. 

As  a  second  example  of  the  instrumentation  required  to  compute  force 
levels  throughout  a  structure  accurately,  consider  a  bridge  deck  composed  of 
many  independent  simple  spans.  If  each  deck  is  considered  to  be  rigid  in  a 
horizontal  plane,  then  all  significant  horizontal  deck  motions  relative  to 
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the  ground  can  be  described  in  terms  of  lateral  and  longitudinal 
translation,  and  rotation  about  a  vertical  axis.  Specifically  for  each 
span,  the  instrumentation  shown  in  figure  4  may  be  used  to  describe  lateral 
deck  translation  relative  to  the  ground,  (Ui(t)  +  U2(t))/2  —  U  (t);  deck 
rotation  about  a  vertical  axis,  (U2(t)  —  Ui(t))/L;  and  longitudinal  deck 
translation  relative  to  the  ground,  Uo(t)  ~U  (t).  Force  levels  in  this 
deck  can  then  be  computed  from  these  generalized  accelerations.  Force 
levels  throughout  the  multispan  bridge  may  be  found  experimentally  only  if 
each  individual  span  is  instrumented  in  a  similar  manner. 

Force  levels  throughout  the  bridge  may  be  found  readily  from  a 
generalization  of  the  relation  F  =  ma.  At  any  instant  in  time  in  the  ith 
coordinate  direction  (whether  in  physical  coordinates,  generalized 
coordinates,  or  normal  coordinates), 

^m  '   l^cl  'hi-°  (23) 

where 

F^  is  an  inertial  force; 

F  ^-  is  a  dissipative  force;  and 

Fl,^  is  a  restoring  force. 

In  general,  the  peak  dissipative,  or  damping,  forces  are  small  with  respect 
to  the  peak  inertial  and  restoring  forces  and  occur  at  different  times 
(Clough  and  Penzien,  1975;  Newmark  and  Rosenblueth,  1971).  Therefore 

Fmil   '-^ki!  (24) 

max      max 

Maximum  restoring  forces  are  of  greatest  concern  because  they  are  related  to 
stresses,  strains,  and  possible  failure.  Since  accelerations  are  recorded  as 
raw  data,  and  since  masses  can  be  computed  quite  accurately,  maximum  inertial 
forces  are  easiest  to  compute. 
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A  general  approach  to  the  computation  of  force  levels  throughout  a 
bridge  takes  three  steps:  (1)  compute  inertlal  forces  acting  on  bridge,  (2) 
compute  reactions  to  those  Inertlal  forces,  and  (3)  compute  any  desired 
Internal  action  from  inertlal  forces  and  reactions.     Forces  will  be  computed 
for  structures  described  by  physical  as  well  as  generalized  coordinates. 
Force  levels  in  bridges  described  in  normal  coordinates  can  be  computed 
exactly  as  for  bridges  (Ascribed  by  generalized  coordinates  If  the  functions 
J-(s)  are  substituted  for  iAs)^  and  nAt)  for  a.(t)1n  the  equations  which 
follow  in  this  section. 

Inertlal  loads  acting  on  bridges  are  computed  ii»st  easily  for  lunped 
mass  models  described  in  physical  coordinates.  At  any  time,  the  inertlal 
force  acting  on  the  ith  mass  is  simp^y 

F^-(t)  =  M-U.(t)  (25) 

where  UAt)  is  the  recorded  absolute  acceleration  at  the  1th  mass.     The 
complete  inertlal  force  distribution  on  the  bridge  is  obtained  by  computing 
inertlal  forces  on  all  masses  for  the  same  Instant  in  time. 

The  computation  of  inertlal  forces  acting  on  a  bridge  whose  distortions 
are  described  in  generalized  coordinates  is  somewhat  more  difficult.    At  any 
Instant  in  time,  t,  for  any  location  along  a  bridge,  s,  the  inertlal   force 
per  unit  length  of  bridge  is  given  by 

F(s,t)  =  M(s)U(s,t)  (26) 

where  U(s,t)  is  the  absolute  acceleration  at  that  point  on  the  bridge.  If 
these  absolute  accelerations  are  expanded  in  terms  of  ground  accelerations 
and  generalized  accelerations 


U(s.t)  =1  ii;,-(s^(t)  +  Uj^g(t)i  +  Uyg(t)j  +  "zg(t)k  (27) 
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then  the  inertia!  forces  can  be  expanded  in  like  manner: 

F(s,t)  =  M(s)  ?  ii;i(s)a^.(t)+M(s)[U^g(t)i+U^g(t)j+U^g(t)k]   (28) 

The  generalized  accelerations,  a.(t),  are  readily  computed  from  the  vector 
of  bridge  distortions  in  physical  coordinates  {u(t)}  using  the  inverse  of 
equation  4,  namely 

•  {a(t)}  =  [i|;r''{u(t)}  (29) 

Inertial  forces  acting  on  the  bridge  can  then  be  computed  using  equation  28. 

The  second  step  in  computing  force  levels  throughout  a  bridge  is  to 
compute  reactions  resisting  those  inertial  forces.  These  can  be  computed 
using  any  acceptable  procedure  of  structural  engineering.  One  such 
technique  using  tlie  principle  of  virtual  work  (Clough  and  Penzien,  1975) 
will  be  discussed  in  detail.  Reactions  will  be  computed  first  for  bridges 
whose  motions  are  described  in  physical  coordinates;  and  second,  in 
generalized  coordinates. 

As  an  example,  the  horizontal  reaction  at  one  end  of  a  simply  supported 
bridge  span  will  be  computed.  This  reaction  and  the  lumped-mass  model  of 
the  span  is  shown  in  figure  5.  Let  r-(s)  be  a  vector  deformed  shape  that 
the  bridge  takes  if  the  left  support  is  deformed  a  unit  amount  corresponding 
to  the  unknown  reaction,  R^(t).  In  addition,  let  r . .  be  the  vector 
displacement  due  to  this  unit  support  displacement  at  the  jth  mass.  The 
r- (s)  and  r- •  are  also  shown  schematically  in  figure  5, 

By  the  principle  of  virtual  work  when  a  bridge,  which  is  in  static  or 
dynamic  equilibrium  under  the  action  of  a  set  of  forces,  is  subjected  to  a 
virtual  displacement,  the  total  work  (force  x  displacement)  done  by  those 
forces  equals  zero  (Clough  and  Penzien,  1975).  The  forces  acting  on  this 
example  bridge  deck  are  all  the  reactions  and  the  inertial  forces 
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Lumped  Mass  Model  And  Recorded  Accelerations 


Deformed  Shape  Due  To  Unit  Displacement 
Corresponding  To  R.j  (t; 


Figure  5.-  Definition  of  shape  functions, 
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F.(t)  =  M.O.(t)  (25) 

If  the  entire  bridge  is  subjected  to  a  single  virtual  displacement,  6r-, 

corresponding  to  the  unknown  reaction,  R-(t),  then  the  total  virtual  work  is 
given  by 

R,.(t)6r^  +jS,Hjiij(t).r.j6r.=0  (30) 

and  for  this  example  the  number  of  masses,  n,  is  five.     There  are  many  other 
non-zero  reactions  besides  the  ith  but  their  virtual   displacements  are  by 
definition  zero,  so  they  produce  no  virtual  work.     Since  6r-   is  small   but 
non-zero  it  may  be  divided  from  equation  30.     Hence 

Such  an  equation  can  be  derived  for  all  reactions,  whether  they  be  forces  or 
moments  (if  the  reaction  is  a  moment,  then  r(s)  is  the  deformed  shape  of  the 
bridge  due  to  a  unit  support  rotation  corresponding  to  the  unknown  moment). 

If  the  distortions  of  the  bridge  Are  described  in  generalized 
coordinates,  it  is  easiest  to  find  all   the  reactions  to  unit  generalized 
accelerations,  and  unit  ground  accelerations,  and  then  linearly  combine  the 
results  for  the  actual   values  of  a^-(t),  'Uw(j(t),  U     (t),  and  '^^At)  measured. 

Consider  first  the  ith  reaction  of  the  same  simply  supported  bridge 

deck  to  a  unit  generalized  acceleration,  a.(t).     Let  this  reaction  be 

R-j,-(t).     This  reaction  and  the  jth  unit  generalized  acceleration  are  shown 

schematically  in  figure  6.     For  this  unit  generalized  acceleration,  the 

inertial   force  per  unit  length  along  the  axis  of  the  bridge 

Fj(s,t)   =  M(s)4)j-(s)aj(t)  (32) 

and  because  a-(t)  =  1 ,  by  assumption, 
u 

Fj(s)  =  M(s);i;j(s)  (33) 

Again  define  r-(s)  as  the  vector  deformed  shape  that  the  bridge  takes  if  the 
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left  support  is  deformed  by  a  unit  amount  which  corresponds  to  the  unknown 

reaction,  R--.  By  the  principle  of  virtual  work,  the  total  work  produced  by 

all  inertial  forces  and  reactions,  when  subjected  to  a  virtual  displacement 

6rj:9   must  equal  zero: 

R.  .6r.  +  /M(s)ii)j(s)-r.(s)6r.ds  =  0  (34) 

Again  since  6r-   is  small  but  non-zero,  it  may  be  divided  from  equation  (34) 

and  so  after  minor  rearrangement, 

Rij  =  -/M(s)?.(s).ii;j(s)ds  (35) 

There  is  no  relation  between  the  number  of  reactions,  i,  and  the  number  of 

assumed  deformed  shapes,  j.  Equation  35  similarly  is  completely  valid  for 

statically  determinate  or  statically  indeterminate  structures.  For 

statically  indeterminate  structures,  the  only  difficulty  arises  in  the 

computation  of  the  r-(s).  Because  typical  highway  grade  separation  bridges 

rarely  have  many  degrees  of  indeterminacy,  this  does  not  seem  to  be  an 

insurmountable  problem.  In  many  cases  these  deformed  shapes  may  be 

approximated  with  suitable  accuracy. 

In  certain  cases  it  may  be  desireable  to  find  the  total  reaction  in  a 

certain  direction  due  to  a  unit  generalized  acceleration  in  the  jth 

direction.  The  principle  of  virtual  work  can  again  be  used  to  advantage. 

For  example,  the  total  reaction  in  the  x-direction  (parallel  to  the  recorded 

ground  acceleration,  0  ) will  be  found.  If  all  supports  are  displaced  by  a 

xg 

unit  amount  in  the  x-direction,  the  entire  bridge  will   undergo  a  rigid  body 
translation,  without  distortion,  in  the  x-direction  by  a  unit  amount.     This 
rigid  body  translation  can  be  described  by  the  vector  function  r(s)  =  i, 
where  i   is  a  unit  vector  in  the  x-direction. 

Again,  the  inertial   force  per  unit  length  of  bridge  due  to  a  unit 
generalized  acceleration  in  the  jth  direction  is 
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F.{s)   =  M{s)ij5j.{s)  (36) 

By  the  principle  of  virtual  work,  the  total  work  produced  by  all   inertial 
forces  and  reactions  must  equal   zero  when  subjected  to  a  virtual 
displacement.     In  calculating  the  total   reactive  force  in  the  x-di recti  on  it 
is  convenient  to  let  the  virtual  displacement  be  a  virtual   rigid  body 
translation  in  the  x-di recti  on. 

6r^(s)  =  r^(s)6r^  =  i6r^  (37) 

Therefore,  the  total  work  produced  by  all   forces  acting  on  the  bridge  is 

where  there  are  assumed  to  be  m  reactions  and  R. .  is  the  ith  reaction  in  the 
x-direction  due  to  a  unit  generalized  acceleration  in  the  jth  direction. 
The  total   reaction  in  the  x-direction  due  to  a  unit  generalized  acceleration 
in  the  jth  direction  is  given  by 


Then 


Rtjx  =  ^Ijx  "  ^2jx  ^  •   •  "  %x-  (39) 


R^j/r^  +  /M(s)Jj(s).  i6r^ds  =  0  (40) 


and  hence 

^tjx  "  -•^f^(s)ifj(s)   •  i   ds  (41) 

Similarly  for  j  and  k  unit  vectors  in  the  y-  and  z-di recti ons,  the  total 
reactive  forces  in  those  two  directions  respectively  due  to  a  unit 
generalized  acceleration  in  the  jth  direction  are 

^tjy  =  -/''Ks)Jj(s)   •  j  ds  (42) 


and 


^tjz  "  -•^^''(s)ij(s)   •   k  ds  (43) 


In  addition  to  the  reactions  produced  by  the  unit  generalized 

accelerations,  there  are  reactions  due  to  unit  rigid  body  translations  in 

the  'Uy„(t),  0  (t),  and  0  (t)  directions.  The  inertial  forces  per  unit 
xy      yg         zg 
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length  of  bridge,  produced  by  unit  ground  accelerations  U    (t)  =1, 

xg 

Uy-(t)  =1,  and  U    it]  -  1  are  respectively 

F^(s)  =  M(s)1  (44) 

ys)  =  HCs)j  (45) 

F^Cs)  =  M(s)k  (46) 

Following  the  logic  used  in  the  confutation  of  the  reaction  R- .,  the  ith 

•J 

reactions  due  to  unit  rigid  body  accelerations  in  the  U  (t),  ii  (t),  and 

U  (t)  directions  are  respectively 

^^   =  -/M(s)i-r.(s)ds  (47) 

R^y  =  -/M(s)j-r.(s)ds  (48) 

and      R^-^  =  -/M(s)k-r^(s)ds  (49) 

It  is  often  advantageous  to  know  the  total  reaction  of  a  rigid  body 

response  to  a  unit  ground  acceleration.  Following  the  logic  used  to  derive 

the  previous  reaction  equations,  and  for  T*i  =1,  the  total  reaction  in  the 

x-di recti on  of  a  rigid  body  response  to  U   =  1  is 

xg 

R^^  =  -/M(s)ds  (50) 

Similarly  the  total  rigid  body  reactions  in  the  y-  and  z-directions  due  to 

unit  accelerations  U„„  =  1  and  U,„  =  1  respectively  are 

yg       ^g 

Rty  =  -/M(s)ds  (51) 

Rtz  =  -/M(s)ds  (52) 

Total  rigid  body  reactions  normal  to  the  direction  of  the  ground 

acceleration  will  obviously  be  zero. 

For  any  instant  in  time,  the  total  ith  reaction  must  include  all  of  the 

effects  from  the  n  generalized  coordinates  and  the  three  rigid  body  ground 

translations.  Specifically, 

R^Ct)  =  R|iai(t)^  — ^^-n^(*^^^-x"xg^*>'^Vyg^*^*'^i'z"zg^*^  ^^^^ 
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Once  the  distribution  of  inertial   forces  over  the  entire  bridge  and  all 
reactions  have  been  found,  the  entire  force  field  on  the  complete  bridge, 
isolated  as  a  rigid  body,  is  known.     These  forces  are  found  using  equations 
25  and  31   if  bridge  distortions  are  described  using  physical   coordinates,  or 
equations  28  and  53  if  bridge  distortions  are  described  using  generalized 
coordinates  (or  normal  coordinates  if  $-(s)   is  substituted  for  ij^-(s),  and 
n.(t)   for  a- (t)).     Since  most  bridges  can  be  modeled  as  a  composition  of 
one-dimensional  elements,  all   internal   actions  acting  on  the  bridge  can  be 
computed  from  statics  alone,  knowing  the  reactions  and  the  distribution  of 
inertial   loads.     Specifically,  for  the  free-body  shown  in  figure  7,  the 
internal   forces  N   ,  V  ,  and  V    and  internal  moments  M  ,  M  ,  and  M    may 
readily  be  found  from  the  six  equations  of  equilibrium  on  the  free-body. 
The  reader  should  note  that  in  most  cases,  distributions  of  forces 
throughout  a  bridge  should  be  computed  using  generalized  or  normal 
coordinates,  rather  than  physical   coordinates  and  a  lumped-mass  model, 
unless  many  masses  spaced  closely  together  are  used  in  the  latter. 

Finding  a  maximum  reaction,  or  a  maximum  internal   action,  may  require  a 
trial  and  error  approach.     In  general,  however,  the  time  where  the  maximum 
occurred  should  be  easy  to  determine.     The  accelerations  can  be  measured  in 
intervals  throughout  this  region  and  tabulated  in  a  systematic  manner. 
Reactions  or  internal   actions  can  then  be  computed  in  a  standardized  manner, 
and  the  maxima  can  easily  be  identified.     These  procedures  are  demonstrated 
in  the  example  analyses. 

C.     Identification  of  the  Mathematical   Model    from  Strong-Motion  Records 

At  best,  only  a  portion  of  the  mathematical  model   can  be  expected  to  be 
identified  by  visual   analysis  of  the  strong-motion  records.     It  is  most 
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reasonable  to  obtain  model  properties  in  normal  coordinates.  In  normal 
coordinates  each  mode  of  vibration  is  characterized  by  the  equation 

n-+2£.a)-n-+a)-n.  =  -v  -U  -y  .U  -y  -U  (54) 

'i  ^1  1  h  1  h    'xi  xg  V'  yg  zi  zg  ^     ' 

where 

oj-        =  natural  circular  frequency  of  vibration 

^-        =  viscous  damping  ratio 

Yxi»Tw^»Y2-i  =  participation  factors 
The  total  relative  displacement  vector  is  a  combination  of  modal  responses 

u(s.t)  =  $^(s)n-,(t)+$2(s)n2(t)+*"+$n(s)%(t)         (55) 

where  the  ^^(s)  are  the  normal  modes  of  vibration.  For  ideal  data  we  can 

expect  to  obtain  only  a  few  natural  frequencies,  go-,  corresponding  mode 

shapes,  $^-(s),  and  corresponding  participation  factors,  y  .  y  .,   and  y  • . 
1  XI ,  y  1      z  1 

Little  information  about  the  damping  ratios  can  be  expected  to  be  found 
without  some  computational  effort  on  a  digital  computer.  Even  this  small 
amount  of  information  gained  may  be  significant.  A  brief  discussion  of  the 
behavior  of  lightly  damped  bridges  will  demonstrate  this. 

Motions  of  bridges  that  behave  in  a  linear  manner  can  be  described  as 
the  sum  of  n  independent  normal  mode  motions  (equation  55).  If  the  ground 
vibrated  continuously  with  a  frequency  equal  to  one  of  the  natural 
frequencies  of  the  bridge,  w-,  the  well-known  phenomenon  of  resonance  would 
occur;  namely,  the  bridge  would  vibrate  in  that  mode  with  an  amplitude  of 
motion  many  times  greater  than  the  ground  motion.  Naturally  occurring 
earthquakes  do  not  vibrate  at  constant  frequencies  for  any  appreciable 
length  of  time,  but  have  vibrational  energy  at  many  frequencies  with 
amplitudes  varying  with  time.  Since  the  ground  motion  is  not  steady, 
resonance  will  not  occur,  but  amplified  response  motions  at  the  natural 
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frequencies  generally  will  occur.  Not  all  modes  of  vibration  will  be 
excited  to  great  amplitudes;  usually  only  the  first  few  or  at  least  those 
with  model  frequencies  matched  by  the  frequency  content  of  the  earthquake's 
energy.  It  follows  that,  even  for  a  \/ery   large  value  of  n,  the  total  bridge 
distortion  often  may  be  well  approximated  by  a  truncated  sum,  such  as: 

u(s,t)  =  ?i(s)a^(t)+52(s)a2(t)+$3(s)a3(t)+e  (56) 

where  e  represents  insignificant  higher  modal  distortions.  If  the 
properties  of  only  a  few  significant  modes  of  vibration  are  found,  the 
significant  motion  of  the  bridge  may  be  modeled,  A  general  rule  of  thumb  is 
that  if  a  mode  responds  significantly  to  an  earthquake,  its  properties  may 
be  found.  The  danger  of  using  only  those  modes  that  are  excited 
significantly  by  one  earthquake  to  predict  the  bridge  motion  to  a  future 
earthquake,  is  that  the  first  earthquake  may  not  excite  a  particular  mode  of 
vibration  that  might  be  excited  significantly  by  a  future  earthquake.  It  is 
therefore  most  important  to  have,  in  addition  to  a  mathematical  model 
deduced  from  the  earthquake  records,  a  mathematical  model  formulated  from 
basic  principles. 

A  mode  will  vibrate  only  at  its  natural  frequency,  w-,  if  excited  at 
that  frequency,  or  if  allowed  to  vibrate  freely  without  any  input  ground 
motion.  A  steady-state  input  motion  with  frequency  oo  will  produce  a 
steady-state  modal  response  also  at  that  frequency,  co,  and  not  at  the 
natural  frequency  of  the  mode,  w- .  Fortunately,  the  modal  response  to  an 
earthquake  is  not  a  great  deal  unlike  the  response  to  a  series  of  impulses. 
Between  impulses  vibration  is  free,  and  hence  at  the  natural  frequency  w-. 
Because  of  the  erratic  arrival  time  of  impulses  (approximating  an 
earthquake)  the  modal  response  will  have  varying  amplitudes,  but  still  will 
be  predominantly  at  the  natural  frequency,  oo-.  If  the  earthquake  lasts  for 
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a  long  time  relative  to  the  natural  period  of  the  mode  (T^-=27t/co^-)  »  then  the 
ground-motion  energy  at  the  frequency  oo-  will  tend  to  amplify  that  modal 
response  on  the  order  of  two  to  four  times  the  amplitude  of  the  ground 
motion  at  a  frequency  of  co-,  depending  on  the  modal  damping.  If  the 
duration  of  the  earthquake  is  short,  that  modal  response  simply  does  not 
have  time  to  be  amplified  to  a  significant  level. 

The  accuracy  of  the  modal  properties  found  depends  upon  how  large  the 
bridge  distortions  are  relative  to  the  ground  motion.  Accel erometers 
measure  absolute  acceleration.  If  the  bridge  response  is  identical  to  the 
ground  motion,  there  are  no  bridge  distortions,  and  nothing  can  be  learned 
about  the  specifics  of  the  dynamic  mathematical  model.  However,  significant 
results  may  still  be  obtained  for  purposes  of  design.  If  there  are  no 
bridge  distortions  there  is  no  dynamic  amplification,  implying  that  an 
equivalent  static  model  may  be  used  to  model  the  earthquake  behavior  for 
this  bridge.  Reactions  may  be  computed  directly  from  the  equations, 
repeated  here, 

Ri^  =  -/M(s)r.(s)-Tds  (57) 

Riy  =  -/M(s)?i(s)-jds  (58) 

Rix  =  -/M(s)r.(s)-kds  (59) 

and  used  to  find  any  internal  actions  as  outlined  previously.  If,  on  the 
other  hand,  the  bridge  distortions  are  significant,  mode  shapes,  natural 
frequencies,  and  participation  factors  may  be  found  as  outlined  beloWo 

Often  it  can  be  assumed,  if  the  relative  response  is  large  with  respect 
to  the  ground  motion,  that  the  relative  response  is  approximately  equal  to 
the  absolute  response,  especially  at  the  dominant  natural  frequencies.  This 
assumption  will  be  made  for  this  example. 
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Dominant  frequencies  of  motion  are  first  identified  and  assumed  to 
equal  natural  frequencies  of  the  bridge.  It  is  important  to  check  the 
ground  motion  to  make  sure  that  the  dominant  frequencies  observed  in  the 
response  are  not  simply  ground  motions  transmitted  without  amplification  to 
the  bridge.  For  the  example  traces  shown  in  figure  8,  two  dominant 
frequencies  of  motion  exist  having  frequencies  o^n^Z-n/T-,   and  Wp=27T/Tp.  Note 
that  if  motions  having  a  few  dominant  frequencies  do  not  occur  consistently 
throughout  all  traces,  the  bridge  is  not  likely  to  be  responding 
significantly  in  its  natural  modes  of  vibration.  In  this  example  the  first 
two  periods  of  vibration  are  estimated  from  a  single  period  of  vibration 
each.  Errors  are  generally  minimized  if  periods  are  measured  from  as  many 
cycles  of  motion  as  possible.  In  figure  8,  the  solid  lines  represent 
recorded  motions,  including  all  modal  responses  whereas  the  dashed  lines 
have  been  drawn  as  estimates  of  the  first  mode  of  response  (lowest  frequency 
mode). 

Mode  shapes  are  obtained  by  comparing  the  relative  amplitudes  of  the 
modal  responses  for  one  instant  in  time.  Usually  that  time  is  taken  at  a 
modal  peak  to  minimize  errors.  At  time  t-.  (corresponding  to  a  peak  in  the 
lowest  frequency  mode)  the  relative  responses  for  traces  1,  2,  and  3  are 
2.0,  3.0,  and  2.0,  respectively.  These  are  shown  plotted  on  a  plan  of  the 
bridge  in  figure  9  from  which  the  continuous  normal  mode  shape  }i(s)  can  be 
estimated.  Since  mode  shapes  depict  relative  displacements,  any  convenient 
scale  may  be  used.  For  instance,  instead  of  the  values  2.0,  3.0,  2.0  the 
values  10,  15,  10;  0.667,  1.000,  0.667;  or  any  other  scale  multiple  may  be 
used.  Once  chosen  they  must  remain  unchanged  for  the  rest  of  the  analysis. 

Similarly,  at  time  to  the  second  modal  peak  values  relative  to  the 
fundamental  mode  response  may  be  measured  as  1.0,  0.0,  -1.0  for  traces  1,  2, 
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Figure  8.-  Example  accelerations. 
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Figure  9.-  Natural  modes  of  vibration. 
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and  3  respectively.  These  values  are  plotted  and  shown  in  figure  9,  from 
which  the  continuous  mode  shape  ^o^^)  ^^"  ^^   estimated. 

Assuming  the  mass  distribution  of  the  bridge  to  be  known,  the 
participation  factors  may  be  computed  from  equations  12  through  15,  expanded 
here: 

Yxi  =  /M(s)?.(s)-ids//M(s)5.(s)-$.(s)ds  (60) 

Yyi  =  /M(s)i^-(s).jds//M(s)$.(s).$.(s)ds  (61) 

Y^i  =  /M(s)$.(s)-kds//M(s);.(s).$.(s)ds  (62) 

If  in  the  above  example  i  is  the  unit  vector  in  the  direction  of  s;  j, 

perpendicular  to  s  in  the  plane  of  the  paper;  and  k,  directed  vertically  out 

of  the  paper,  then  the  estimated  mode  shapes  can  be  rewritten  as 

?l(s)  =  (|)i(s)j  (63) 

$2(s)  =  (l)2(s)j  (64) 

Because  i*j  =  i*k  =  j*k  =  0  then 

Yyi  =  /H(s)<t.,(s)ds//H(s)<t,,^(s)ds  (56) 

Yy2  =  /M(s)(f2(s)ds//M(s),f2^(s)ds  (67) 

Y^l  -\z-0  (68) 

where  the  integrals  no  longer  contain  vector  functions. 

The  damping  ratios  E,-,   and  Co  cannot  be  obtained  visually  from  the 
induced  vibrations  as  sketched.  Only  in  certain  circumstances  can  damping 
ratios  be  obtained. 

If  response  spectra  (Clough  and  Penzien,  1975;  Newmark  and  Rosenblueth, 
1971;  Wiegel ,  1970)  are  published  for  the  earthquake  and  the  site  in 
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question  (or  if  published  spectra  are  assumed  to  be  valid  for  the  site)  then 
modal  damping  ratios  may  be  found  by  comparing  measured  peak  modal 
accelerations  to  peak  values  for  various  damping  ratios.  The  damping  ratio 
at  which  the  two  peak  modal  responses  coincide,  is  assumed  to  equal  the 
actual  modal  damping  ratio.      ' 

Specifically,  a  response  spectrum  is  a  collection  of  points,  each  one 
of  which  is  the  peak  response  of  a  simple  oscillator  governed  by  the 
equation  of  motion  with  parameters  co  and  E, 

y+2E,i^y+Jy   =  -U^^  (or  U^^  or  U^^)  (69) 

subjected  to  the  single  ground  motion  input.  A  typical  set  of  response 

spectra  for  a  single  ground  motion  is  shown  in  figure  10.  If  the  modal 

response  is  excited  by  a  single  component  of  the  ground  motion,  as  in  this 

example,  then  the  damping  ratios  can  be  found  as  follows.  Let  fin  ^^  and 

I  max 

Tip  ^3  be  the  maximum  modal  accelerations  measured,  corresponding  to  the 
mode  shapes  $i(s)  and  ^o^^)*  The  maximum  corresponding,  simple 
oscillator  responses  are 

h^''z1   =  'iz  max^^y2  (^1) 

By  comparing  these  spectral  values  measured  to  the  spectral  values  computed 
at  frequencies  o),  and  cop  it  can  be  deduced  that  the  modal  damping  ratios  are 

^   =   0.05  (5%) 

^2  =  0.03  (3%) 
If  the  modal  responses  are  the  result  of  multiple  ground  motion  inputs  (more 
than  one  non-zero  participation  factor)  then  the  response  spectrum  method  of 
determining  modal  damping  ratios  cannot  be  used. 

It  should  be  noted  that  response  spectra  cannot  be  computed  easily 
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Figure  10.-  Response  spectrum  for  U 


yg' 
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without  the  aid  of  a  digital  computer.  If  a  response  spectrum  is  not 
published  by  others  from  nearby  ground  motions,  it  is  not  reasonable  to 
pursue  this  analysis  approach.  If  the  traces  are  digitized  in  order  that  a 
response  spectrum  be  computed,  general  computer  aided  methods  of  analysis  of 
the  strong-motion  records  (described  in  the  next  section)  preferably  may  be 
used  to  obtain  modal  properties. 

A  second  method  of  obtaining  damping  ratios  visually  from  analog  records 
is  possible  if  there  exists  a  distinct  free- vibration  portion  of  the  response 
at  the  end  of  the  record.  For  some  earthquakes,  particularly  near  the 
source,  the  ground  motion  ceases  shortly  after  it  begins.  If  at  that  time 
of  cessation,  the  relative  bridge  response  is  significant,  it  will  continue 
to  vibrate  freely  until  all  kinetic  energy  is  dissipated  by  the  damping 
mechanisms.  One  of  several  possible  modal  free  vibrations  is  sketched  in 
figure  11. 

The  envelope  of  a  viscously  damped  free  vibration,  shown  by  the  dashed 
line,  is  given  in  terms  of  the  frequency  and  damping  ratio  by  e  ^^  (Clough 
and  Penzien,  1975).  For  C  equal  to  some  constant,  the  ratio  of  peak  modal 


amplitudes  at  the  times  t  and  t+nT  is 

n(t)        ^__C£^ 
n(t+nT)  cg-5w(t+nT) 


n(t)        = C£^ (72) 


=  L„-_ =  e^^"^  (73) 

g-?wtg-ga3nT  ^'^^ 

If  natural   logarithms  are  taken  of  this  ratio  then 

'ln(n(t)/n(t+nT))  =  ^wnT  (74) 
or  finally 

^  =  -^^M^{t)/r){t+nl))  (75) 
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Figure  11.-  Example  free  vibration. 


43 


for  T  =  Ztt/w.  Unfortunately,  the  ground  motion  often  continues  to  be 
significant  for  a  sufficient  length  of  time  so  that  no  distinct, 
free-vibration  portion  of  the  record  exists.  It  may  be  advisable  to  test 
the  bridge  to  find  these  modal  damping  ratios.  The  bridge  may  be  perturbed 
to  produce  a  free  vibration  by  a  heavy  truck  stopping  abruptly  on  the 
bridge,  or  if  convenient  it  may  be  pulled  and  released  by  a  tractor 
(Douglas,  1976;  Newmark  and  Rosenblueth,  1971;  Wiegel ,  1970).  In  both  cases 
the  free  vibrations  recorded  may  be  analyzed  by  the  method  outlined  above. 

It  should  be  noted  that  force  levels  throughout  the  bridge  may  be 
obtained  using  normal  modes  and  normal  coordinates  exactly  as  outlined  in 
the  beginning  of  this  section.  In  place  of  the  assumed  deflected  shapes 
i|)-j(s),  the  normal  modes  $j(s)  are  used.  Similarly,  in  place  of  the 
coordinates  a^-(t),  the  normal  coordinates  rii(t)  are  used. 

D.  Example  Analyses 

In  this  section,  hypothetical  strong-motion  records  for  two  bridge 
designs  will  be  analyzed.  For  both  examples,  absolute  bridge  accelerations 
have  been  generated  theoretically  (using  the  finite-element  program  SAP  IV 
(Bathe,  Wilson,  and  Peterson,  1973))  from  actual  recorded  ground  motions  at 
the  D.  Johnson  Ranch,  Oroville,  Calif,  on  6  August  1975  (Sherburne  and  Hauge, 
1975).  These  theoretically  generated  acceleration  traces  are  then  analyzed 
as  if  they  had  been  recorded  as  raw  data. 

For  both  examples,  force  levels  are  found  as  well  as  the  dynamic 
response  characteristics.  The  two  examples  include  a  bridge  of  three  simple 
spans,  and  a  continuous  two-span  bridge.  As  will  become  evident,  the  motion 
of  a  simple  continuous  two-span  bridge  is  yery   complex.  Although  it  is 
entirely  feasible,  it  is  impractical  to  analyze  more  complex  bridges 
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visually  from  analog  records  in  as  systematic  a  manner  as  is  done  for  these 
two  examples. 

Three-Span  Simply-Supported  Bridge  -  The  subject  bridge  and  the  strong- 
motion  instrument  locations  are  shown  in  figure  12.  The  absolute 
accelerations  recorded,  corresponding  to  the  instrument  locations,  are 

U„„»U„„,U,„       ground  accelerations 
xg  yg  zg      ^ 


^T^2'^3         structure  accelerations 

Because  of  symmetry,  the  two  lateral  pier  cap  motions  are  assumed  to  be 
identical  so  only  one  need  be  recorded.  Shown  in  figure  13  are  the 
three  recorded  ground  motions,  and  in  figure  14  the  three  hypothetical 
structure  motions  Up  U29  and  Do. 

Before  the  traces  can  be  analyzed,  some  assumptions  regarding  the 
bridge  behavior  must  be  made.  It  is  assumed  that  the  three  deck  elements 
behave  as  rigid  plates,  restrained  longitudinally  at  one  end,  and  restrained 
vertically  and  laterally  at  both  ends.  It  is  expected  that  vertical 
flexural  motions  will  exist,  but  they  cannot  be  identified  with  the 
instrumentation  shown  on  figure  12.  It  is  also  assumed  that  the  vertical 
flexural  motions  are  uncoupled  from  the  lateral  and  longitudinal  motions  to 
be  analyzed.  Torsional  motions  about  the  bridge  center! ine  axis  are  also 
ignored  for  the  following  two  reasons.  First,  torsional  components  of  the 
ground  motion  are  not  likely  to  be  large.  And  second,  secondary  torsional 
motions  of  the  bridge  deck  due  to  lateral  translations  will  not  be  excited 
at  the  bents  because  the  pair  of  columns  at  each  bent  prevents  the  pier  cap 
from  rotating  about  the  centerline  axis.  Although  the  lateral  and 
longitudinal  ground  motions  are  recorded  at  a  single  point  only,  it  is 
assumed  that  both  abutments  and  all  pier  bases  experience  motions  identical 
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Figure  12.-  Three-span  simply  supported  bridge, 

geometry  and  strong-motion  instrumentation. 
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Figure  13.-  August  6,  1975  earthquake  aftershock  ground  motion  recorded  at  the 
D.  Johnson  Ranch,  Oroville,  Calif. 
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Figure  14.-  Three-span  simply  supported  bridge, structure  motions. 
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to  those  recorded. 

Because  of  the  simple  and  discontinuous  nature  of  this  bridge,  physical 
coordinates  are  used  and  masses  corresponding  to  bridge  degree  of  freedom 

are  assumed  to  be  lumped  at  the  accelerometer  (and  coordinate)  locations. 

2 
The  mass  of  a  column  cap  and  the  two  columns  is  3,863  kip-sec  /ft,  and  the 

2       2 
mass  per  foot  of  deck  is  0.303  kip-sec  /ft  .     The  masses  corresponding  to 

each  degree  of  freedom  are  shown  shaded  in  figure  15. 

(1)  Maximum  Force  Levels 

For  this  example  the  maximum  force  reactions  experienced  will  be  found. 
Since  the  degrees  of  freedom  are  both  inertially  as  well  as  elastically 
uncoupled,  maximum  force  reactions  (figure  16)  corresponding  to  each  degree 
of  freedom  are  found  simply  as  the  product  of  mass  (figure  15)  times  the 
corresponding  peak  absolute  acceleration.  The  reactions  to  be  found  are 
shown  in  figure  16.  The  mass  M  .  is  used  in  calculating  lateral  and 
vertical  reactions  at  the  ends  of  the  bridge. 

The  maximum  accelerations,  masses,  and  maximum  reactions  are  shown  on 
table  1.  Note  that  the  maximum  reactions  do  not  occur  at  the  same  instant 
but  at  many  different  times.  At  any  given  instant,  all  reactions  and 
distributed  inertial  forces  could  be  found  in  a  similar  manner.  Then  from 
equations  of  dynamic  equilibrium  for  any  free-body  all  internal  member 
forces  could  be  computed. 

(2)  Structural   Model   Parameters 

Because  of  the  simple,  uncoupled  behavior  of  this  bridge,  structural 
model  parameters  may  be  found  rather  easily  from  a  hand  analysis  of  the 
records. 

As  described  in  the  previous  section,  the  mathematical  model  is 
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Figure  15.-  Three-span  simply  supported  bridge  mass. 
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Figure  16.-  Three-span  simply  supported  bridge, 
unknown  reactions. 
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Reaction      Acce 

leration 

Max  Value 

Mass     ^ 

Force 

(fps^) 

(ill 

■p-sec  X 
ft        ^ 

(kips) 

1                        1 

^1 

3.0663 

64 

.463 

197 

661 

2 

h 

10,634 

26 

.588 

282 

.742 

3 

^3 

45.830 

34 

.163 

1565. 

59 

4 

^3 

45.830 

34 

.163 

1565 

.69 

5 

^9 

9.899 

11 

.363 

112 

.48 

6 

"^9 

9.899 

11 

.363 

112 

.48 

7 

5.439 

34 

.163 

185 

.82 

..Z9 

8 

'.^, 

5.439 

34 

.163 

185 

.81 

9 

'.^9 

5.439 

11 

.363 

61 

.81 

10 

"zq 

5.439 

n 

.363 

61 

.81 

Table  1 .-Three-Span  Simply  Supported  Bridge,  Maximum 
Reactions 
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inertially  uncoupled,  i.e.,  M- •  =  0  for  all   i  ^  j.     Due  to  the  expansion 
joints  and  structural   symmetry  of  this  bridge,  it  is  obvious  that  the  bridge 
also  is  elastically  uncoupled,  that  is,  K-.  =  0  for  all   i  ^  j.     It  is  also 
assumed  that  the  bridge  motions  are  not  coupled  by  off-diagonal  damping 
terms,  C--.     Therefore  all   bridge  distortions  may  be  modeled  by 


M^  0     0 


0    M2  0 
0     0     M' 


^^1^ 


no>+ 


C^     0     0 


0       C2  0 
0       0     C 


^1 


<n 


2 


>+ 


K^      0     0 


0       /C2  0 
0       0^- 


S' 


^M^ 


>=  -< 


Mr 


^0' 


^  \  -< 


ii  \  ^y 


[' 


U 


yg 

(76) 


where  n,-  are  relative  modal  displacements  (in  what  are  already  normal 
coordinates)  defined  by 

ni(t)  =  U^(t)-U^g(t)  (77) 

'^2(t)  =  U2(t)-Uy^(t)  (78) 

n3(t)  =  lJ3(t)-U^g(t)  (79) 


Equation  (76)  rearranged  becomes 

fi^.+2^.a).f^.^ni  =  -Txi*^xg-Yyi^g 
where  again 


03^.=  /K./M. 

and  for  this  example 

^   =   C^./2/iCjM. 

and  for  this  example 

^xl  =  "•    ^x2  = 

1 

\3-' 

^yl  =  ^   ^y2  = 

0 

'y3-' 

(80) 


(81) 


(82) 


It  is  intended  to  find  both  w-  and  5.  from  the  strong-motion  records;  and 
then  knowing  M- ,  C.  and  K-   may  be  found  from  equations  81  and  82. 
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Figure  13  indicates  that  approximately  3  seconds  into  the  earthquake, 
the  ground  motion  has  died  out,  and  the  bridge  vibrated  freely  without 
additional  significant  ground  motion  input.  Natural  frequencies  can  simply 
be  measured  by  counting  the  cycles  in  this  free  vibration  period,  and  the 
damping  can  be  estimated  from  the  logarithmic  decrement  of  the  decaying 
records.  For  t>3  seconds,  it  is  assumed  that  U   =  U   =  0  so  n-  =  IJ- 
(i=l,2,3).  Shown  in  table  2  is  a  summary  of  the  frequency  information 
obtained  from  figure  14.  Similarly  damping  ratios  are  found  by  computing 
the  logarithmic  decrements  for  t>3  seconds.  These  damping  ratios  are 
summarized  in  table  3. 

Shown  in  tables  2  and  3  are  the  actual  values  of  oj-  and  E,-   which  were 
used  to  generate  these  hypothetical  records  for  this  example  analysis.  As 
can  be  seen  the  natural  frequencies  can  be  found  with  excellent  accuracy, 
but  the  damping  ratios  with  less  accuracy.  The  free  vibration  portion  of 
these  records  are  exceptionally  clean  considering  that  they  were  created 
from  real  earthquake  motions.  Even  for  these  clean  records,  however,  the 
damping  ratio  estimates  are  not  particularly  accurate.  Estimating  dampiiiy 
from  less  clean  strong-motion  records  is  likely  to  be  highly  inaccurate. 

Since  mode  shapes  are  relations  between  structure  displacements  for  a 
given  natural  mode,  and  since  these  bridge  motions  are  naturally  decoupled, 
no  mode  shapes  exist  as  such.  More  precisely,  in  terms  of  displacements  n-i* 
rjos  Ho  the  mode  shapes  are 


1 
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0 

V  J 


i<^o}    = 


(   \ 

0 

^> 
0 


{^3}  = 


0 
0 
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Two-Span  Continuous  Bridge  -  For  this  example  force  levels  and  structural 
model   parameters  will   be  found  from  hand  analyses  from  the  strong-motion 
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Degree  of       No.   of        Time         Period      oo.  ^a^^-Moi 

Freedom  cycles  ^  ^^^"^' 


1 

10 

5.103 

.510 

12.313 

12.236 

2 

14 

4.600 

.329 

19.123 

19.052 

3 

25 

4.500 

.183 

34.282 

33.615 

Table  2.-  Three-Span  Simply  Supported  Bridge, 
Natural    Frequencies 


DOF       N0.^0f  .,„         n,„  ,,         .(,,,„,,) 


1 

9 

1.30 

.18 

0.035 

0.040 

2 

15 

6.30 

1.20 

0.018 

0.040 

3 

20 

21.00 

2.00 

0.019 

0.020 

Table  3.-  Three-Span  Simply  Supported  Bridge, 
Damping  Ratios 
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records  of  a  continuous  two-span  grade  separation  bridge.  Hypothetical 
absolute  accelerations  have  again  been  generated  at  the  proposed  instrument 
locations  and  again  will  be  analyzed  as  if  they  had  actually  been  recorded. 

The  bridge  to  be  analyzed  is  shown  schematically  in  figure  17.  Shown 
in  figure  18  are  the  instrument  locations.  The  ground-motion  time  histories 
used  are  the  same  ones  used  for  the  previous  example  and  are  shown  in  figure 
13.  The  hypothetical  bridge  motions  which  are  to  be  analyzed  are  shown  in 
figures  19,  20,  and  21. 

A  few  general  assumptions  are  again  made  with  regard  to  the  expected 
behavior  of  this  bridge.  Although  ground  motions  have  been  recorded  at  a 
single  point,  it  has  been  assumed  that  the  abutment  motions  and  the  column 
base  motions  are  identical  to  those  ground  motions  recorded.  Furthermore, 
it  is  assumed  that  the  vertical  motion  of  the  bridge  deck  directly  above  the 
central  column  equals  the  vertical  ground  motion.  The  ends  of  the  bridge 
deck  are  free  to  move  longitudinally,  are  free  to  rotate  about  a  transverse 
axis  and  are  free  to  rotate  about  a  vertical  axis,  but  are  fixed 
transversely,  vertically,  and  rotationally  about  a  longitudinal  axis.  The 
central  column  is  assumed  to  be  fixed  at  both  its  top  and  bottom. 

(1)  I^ximum  Force  Levels  -  For  this  example  the  maximum  force  reactions 
will  again  be  found  from  a  hand  analysis  of  the  hypothetical  bridge  motions. 
Those  reactions  sought  are  shown  in  figure  22.  Because  the  example  bridge 
is  continuous,  it  is  not  advisable  to  use  physical  coordinates  and  a  lumped 
mass  model;  instead  it  is  advisable  to  use  generalized  masses  and  a 
corresponding  generalized  coordinate  system  as  described  Section  II.  Such  a 
coordinate  system  can  describe  structure  distortions  accurately  with 
relatively  few  coordinates. 
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Figure  17.-  Continuous  two-span  bridge. 
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Figure  18.-  Continuous  two-span  bridge,  strong-motion 
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Figure  19.-  Continuous  two-span  bridge  structure  motions. 
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Figure  21.-  Continuous  two-span  bridge  structure  motions. 
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Figure  22.-  Continuous  two-span  bridge, desired  force  reactions, 
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It  is  assumed  that  all  structure  distortions  can  be  described 

adequately  as  a  linear  combination  of  the  assumed  distortions  shown  in 

figure  23.  (Note  that  it  is  assumed  that  all  structure  motions  can  be 

described  adequately  as  a  linear  combination  of  these  assumed  distortions 

and  the  rigid  body  translations.)  The  vectors  i,  j,  and  k  are  unit  vectors 

in  the  x-,  y-,  and  z-directions  respectively.  For  use  in  calculating 

reactions  from  the  above  distortions,  distortions  due  to  unit  displacements 

corresponding  to  the  unknown  reactions  must  be  found  or  assumed.  Those 

distortions  needed  are  shown  in  figure  24.  In  addition,  for  use  in 

computing  the  total  reactive  forces  in  the  three  mutually  perpendicular 

directions,  the  distortions  r  (x)=i,  r  (x)=j  and  r  (x)=k  are  used. 

X  y  z 

Reactions  to  unit  generalized  accelerations  are  obtained  using  equation 
35  and  a  simple  numerical   integration  procedure.     For  symmetric  iAx),  the 
central   reaction  is  found  using  equation  35.     Each  end  reaction  is  then- 
found  as  half  of  the  difference  between  the  total   reaction  and  the  central 

reaction.     For  the  mass  per  unit  length  of  deck  equal   to  0.3295 

2       2 
kips-sec  /ft  ,  and  ignoring  the  column  mass,  the  force  reactions  to  unit 

2 
values  1   ft/sec     of  the  generalized  accelerations  are  given  in  figure  25. 

Similarly  reactions  for  unit  rigid  body  motions  are  found  as  shown  in  figure 

26.     The  reactions  in  figure  26  are  found  using  equations  47,  48,  and  49. 

Now  that  all   reactions  due  to  unit  generalized  and  rigid-body 

accelerations  have  been  found,  it  remains  to  find  the  actual   generalized  and 

rigid-body  accelerations  experienced  by  the  bridge  during  the  hypothetical 

earthquake.     For  the  calculation  of  maximum  force  levels,  the  recorded  time 

histories  have  been  expanded  to  show  only  four  seconds  of  the  strongest 

motion.     These  expanded  time  histories  are  shown  in  figures  27,  28,  29,  and 

30. 
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Fiaure  27.-  Continuous  two-span  bridge,  ground  accelerations. 
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Fiaure  28.-  Continuous  two-span  bridge,  structure  accelerations. 


69 


980 


o 

+->(\J 
(O    o 

0)    CO 

o  o 

o 


■980 -J 


980 


O 
+JC\J 

n3   O 

S-    QJ 
<U    «/) 


O    O 
O^ 


-  -  ---\---v-\-- 


980 


344  0  1.934  sec 


Time  (seconds) 
Lateral  Acceleration  at  Pier  Cap,  U.(t) 


v/  V  V  \/  ^/  y  'A, 


Time  (seconds) 
Vertical  Acceleration  at  West  Edge  of  Deck  at  Pier  Cap,  u'  (t) 


Figure  29.-  Continuous  two-span  bridge,  structure  accelerations. 
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(a)  Longitudinal  Reactions 

Rigid  body  motion  in  the  x-di recti  on  and  distortions  in  the  first 
generalized  coordinate  direction  only  contribute  to  longitudinal  reactions. 
The  extreme  absolute  acceleration  in  the  U,  direction  is  given  by 

'^1 'extreme"  -232.34  cm/sec^  §  t  =  1.72  seconds 
At  t  =  1.72  seconds,  the  longitudinal  ground  acceleration  is  given  by 

'^vn  =  ^^3  cm/sec^ 
By  definition  the  corresponding  generalized  acceleration  is  obtained  from 

U-j  =  U^  +^-j(150)a^     (see  equation  27) 

where  (j;-j(150)  =  1  is  the  amplitude  of  the  first  assumed  distortion  at  the 
instrument  1  location.  Therefore 

^1  =  '^r\g  (83) 

and  at  t  =  1.72  seconds 

2 

d-i  =  -715.34  cm/sec 

From  figures  25  and  26  and  equation  83  the  total  longitudinal  reaction  is 

found  as  the  linear  combination  of  reactions  from  a-, (1.72)  and  U  (1.72). 

I  xg 

This  maximum  reaction  of  753.5  kips  is  shown  graphically  in  figure  31. 

Although,  in  this  case,  this  reaction  could  have  been  found  directly  from 

2 
U-|Ly.^^_,=  -232.34  cm/sec  ,  it  is  of  value  to  decompose  this  force  into 

components  associated  with  structural   distortions  and  rigid  body  motions. 

(b)  Lateral    (transverse)  Reactions 

Rigid  body  motion  in  the  y-direction  and  distortions  in  the  fourth  and 
fifth  generalized  coordinate  directions  contribute  to  lateral   bridge  reactions 
(Rpj  ^2'  ^^^  ^4^*     Maximum  lateral   forces  are  likely  to  occur  when  l)^  and  Dg 
are  at  or  near  their  extreme  values.     Both  have  their  extreme  values  at 
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Figure  31.-  Continuous  two-span  bridge,  maximum  longitudinal  reactions. 
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t  =  1.93  seconds  (see  figures  29  and  30) 

'J4lextreme=  "^^^-^^  ^/^ec^ 

iJ6lextre,«=  -220-^9  cm/sec^ 
At  t  =  1.93  seconds 

U   =217  cm/sec^ 
By  definition  the  absolute  accelerations  at  the  instrument  locations  4  and  6 
can  be  expanded  in  terms  of  the  generalized  coordinates  by 

U4  =  %+i|^4(150)a4+i|;3(150)a5  (84) 

Ug  =  Uyg+ip4(225)a4+i|j5(225)a5  (85) 

Using  the  measured  values  of  U.,  (jg  and  0   at  t  =  1.93  seconds,  and  the 
known  distortions,  equations  84  and  85  are  equivalent  to 

1. 0000x4(1  o93)-1.000ag(l, 93)   =  -561.33 
.707a4(lo93)+  .7070^^(1.93)   =  -437.59 
from  which 

a4(1.93)  =  -590.14  cm/sec^ 

0x5(1.93)  =  -28.81  cm/sec^ 
and  * 

U     (1.93)   =  217.00  cm/sec^ 
Since  lateral   reactions  are  presently  known  for  'dn='d^=\S^    =^    (foot/sec*^), 
lateral   reactions  to  the  above  accelerations  at  t  =  1.93  seconds  can  be 
scaled  and  combined  as  shown  in  figure  32. 
(c)     Vertical   Reactions 

Rigid  body  motion  in  the  z-di recti  on  and  distortions  in  the  second  and 
third  generalized  coordinate  directions  contribute  to  vertical   force 
reactions.     Extreme  vertical   reactions  are  likely  to  occur  when  Up  and/or  U^ 
are  extreme.     Note  that  extreme  seismic-induced  reactions  only  are  being 
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Figure  32.-  Continuous  two-span  bridge,  maximum  lateral  reactions. 
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found  and  no  attempt  is  made  to  include  dead  loads  in  this  analysis. 

The  times  and  accelerations  scaled  from  figures  27  and  28  (and  shown  in 
table  4)  are  possible  candidates  for  producing  maximum  vertical  reactions. 
By  definition,  the  structure  accelerations  can  be  expanded  in  terms  of 
ground  accelerations  and  generalized  accelerations  by 


U2  =.U^g+i|;2(75)a2+ip3(75)a3 

U3  =  ij^  +11^2(225)012+11^3(225)013 

From  figure  23 

1^2(75)  =     1.000  1^3(75)  =  1.414 

i|;2(225)   =  -1.000  11^3(225)   =  1.414 


(86) 
(87) 


so 


>    = 


1.000     1.414 
-1.000     1.414 


f        \ 

(ip 

i 

cip 
\     J 

(88) 


or  inverted  (Gere  and  Weaver,  1965) 
0.500         -0.500 
a„  0.354  -.354 


a 


>  = 


2  zg 

I  3  zgj 


(89) 


The  peak  generalized  accelerations  are  then  computed  from  equation  78,  and 
the  scaled  values  are  shown  in  table  4,  Peak  generalized  accelerations  are 
shown  in  table  5. 

The  reactions  R^,  Rg,  Ry  sought  are  shown  in  figure  33.  At  any  time, 
the  vertical  reactions  are  found  as  the  linear  combination  of  reactions  due 
to  unit  generalized  accelerations  a^  and  a-,  and  a  unit  vertical 
acceleration  U  .  Shown  in  table  6  are  the  vertical  reactions  for  the  five 
candidate  times  for  producing  extreme  vertical  reactions.  As  can  be  seen 
from  table  6,  the  extreme  seismically  induced  value  of  Rn  (-96.58  kips) 
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Accelerations   (cm/ 5 

;ec^) 

Time 

^2 

^3 

'z. 

1.72 

105.00 

-164.90 

165.79 

1.96 

-133.93 

149.03 

-85.00 

2.06 

128,33 

-168.72 

148.47 

2.76 

139.13 

-119.00 

-22.67 

2.97 

-136.73 

107.67 

17.00 

Table  4.-  Continuous  Two-Span  Bridge,  Peak 
Vertical   Accelerations 


2 
Acceleration  (fps   ) 

Time 

^2 

^3 

^zg 

1.72 

4.427 

-4.547 

5.439 

1.96 

-4.462 

2.150 

-2.789 

2.06 

4.873 

-3.918 

4.871 

2.76 

4.235 

.760 

-.744 

2.97 

-4.009 

-.732 

.558 

Table  5.-  Continuous  Two-Span  Bridge,  Peak 
Generalized  accelerations 
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Figure  33.-  Continuous  two-span  bridge,  desired  vertical  reactions, 
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Time 

Source 

Reactions  (kips) 

^7 

1.72 

^12% 

-69.64 

0 

69.64 

^•3^3 

78.39 

224.76 

78.39 

Total 

-97.68 

-342.28 

-97.68 

-88.93 

-117.52 

-50.35 

1.96 

R^2^2 

70.19 

0 

-70.19 

^ish 

-37.07 

-106.27 

-37.07 

Total 

50.09 

175.51 

50.09 

83.21 

69.24 

-57.17 

2.06 

R.2S2 

-76.65 

0 

76.65 

^•3^3 

67.55 

193.67 

67.55 

f^iz'^zg 
Total 

-87.48 

-306.53 

-87.48 

-96.58 

-112.86 

56.72 

2.76 

^•2^2 

-66.62 

0 

66.62 

^^13^3 

-13.10 

-37.57 

-13.10 

^'z^zg 
Total 

13.36 

46.82 

13.36 

-66.36 

9.25 

66.88 

2.97 

f^i2^2 

63.06 

0 

-63.06 

f^i3^3 

12.62 

36.18 

12.62 

^•z'^zg 
Total 

-10.02 

-35.11 

-10.02 

65.66 

1.07 

-60.46 

Table  6.-  Continuous  Two-Span  Bridge,  Peak 
Vertical   Reactions 
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occurs  at  t  =  2.06  sec,  of  Rg  (-117.52  kips)  occurs  at  t  =  1.72  sec,  and  of 
Ry  (66.88  kips)  occurs  at  t  =  2.76  sec.  Total  extreme  loads  must  include 
the  dead  load  of  the  structure. 

(2)  Structural  Model  Parameters 

Structural  model  parameters  can  best  be  obtained  by  hand  analysis  in 
normal  coordinates  only.  In  normal  coordinates,  the  model  is  completely 
described  by  natural  frequencies,  damping  ratios,  mode  shapes,  and 
participation  factors.  Because  of  the  many  modes  of  vibration  which 
respond,  even  for  this  yery   simple  two-span  continuous  bridge,  it  is 
difficult  to  estimate  damping  from  the  free-vibration  portions  of  the 
structure  motions.  As  a  consequence,  only  natural  frequencies  and  natural 
mode  shapes  are  expected  to  be  found  from  this  hand  analysis.  Participation 
factors  are  readily  computed  from  known  mass  distribution  and  mode  shapes 
using  equations  60,  61,  and  62. 

Natural  modes  of  vibration  are  generally  sought  in  regions  of  the 
record  where  there  is  relatively  little  input  motion.  A  natural  mode  will 
vibrate  at  the  forcing  frequency  when  forced  to  vibrate,  but  after  excited 
it  will  vibrate  at  a  unique  natural  frequency.  For  that  natural  mode  of 
vibration,  all  parts  of  the  structure  will  vibrate  in  phase  at  its  natural 
frequency.  For  a  single  mode  of  vibration  the  relative  amplitudes  of  motion 
for  any  two  points  will  remain  constant  for  all  time.  These  properties  of  a 
mode  will  be  used  to  identify  their  presence.  It  should  be  noted  that, 
precisely  speaking,  normal  modes  of  vibration  are  mathematical  entities 
only.  It  is  fortunate,  however,  that  the  distortions  of  structures  in 
general  can  be  expanded  in  terms  of  natural  modes  quite  accurately  for  some 
structures  and  reasonably  accurately  for  most  structures. 
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It  is  unrealistic  to  assume  that  natural  modes  of  vibration  will  be 
found  in  all  analyses.  It  is  particularly  unrealistic  to  assume  that  they 
will  be  found  for  large  amplitude,  highly  nonlinear,  damaging  motions. 
Unfortunately,  it  is  often  the  case  that  these  motions  are  the  most 
interesting  of  all,  and  yet  are  the  most  difficult  to  analyze  and  are  rarely 
recorded. 

The  characteristics  of  normal  modes  described  above  will  be  used  to 
identify  the  existence  of  normal  modes  in  the  analysis  to  follow.  Five 
modes  of  vibration  are  identified  from  the  hypothetical  response  time 
histories  shown  in  figures  19,  20,  and  21.  The  expanded  time  histories 
shown  in  figures  28,  29,  and  30  are  also  used  to  find  structure  model 
parameters  for  t<  4  sec.  The  identified  modes  of  vibration  are  as  follows: 
Mode  1  -  Observed  in  Up  and  U^  for  t>  8  seconds. 

It  is  an  asymmetric  vertical  flexural  mode  of  vibration. 
Mode  2  -  Observed  in  'U-, ,  Up,  and  U^  for  2  1/2  <  t  <  6  seconds.  It  is  an 
asymmetric  vertical  and  longitudinal  flexural  mode  of  vibration, 
Mode  3  -  Observed  in  U^  and  Ug  for  t  >  2  seconds.  It  is  an  assumed 
symmetric,  lateral,  flexural  mode  of  vibration  that 
probably  has  a  torsional  component  which  is  not 
observed. 
Mode  4  -  Observed  in  Up  and  U^  for  t  <  2  seconds.  It  is  a  symmetric 

vertical  mode  of  vibration. 
I^de  5  -  Observed  in  Ur  and  U^  for  t  >  2  seconds.  It  is  an  assumed 
symmetric  torsional  mode  of  vibration  that  probably 
has  a  lateral  component  which  is  not  observed, 
(a)  Mode  Number  1 

Information  needed  to  determine  properites  of  this  mode  of  vibration 
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are  obtained  from  Up  and  Uo.  Because  the  observed  motions  in  this  mode  are 

180%  out  of  phase  (asymmetric  vertical  flexure)  it  is  likely  that  there 

should  exist  some  longitudinal  component  to  the  mode,  but  none  is  observed 

in  U^. 

Natural    frequency  information  is  obtained  from  both  Up  and  U^  from  t  = 

7.3  seconds  to  IKS  seconds.     For  this  time  increment  both  U      and  U       are 

xg  zg 

assumed  to  be  negligible,  and  so  the  response  is  essentially  a  free 
vibrationo     Frequency  information  is  tabulated  in  table  7;  the  natural 
frequency  oo-,   is  computed  from  the  equation  oo-,   =  Ztt/T-,  where  T-,   =  (time  for  n 
cycles )/n. 

Mode  shape  information  is  obtained  by  comparing  the  relative  amplitudes 
of  Up  cind  U^for  free  vibrations  at  this   frequency.     From  the  average  of 
sixteen  peak  values  in  the  above  time  period,  for  U^  arbitrarily  given  the 
value  -1.000,  then  Up  =  1.025.     Mode  1   characteristics  are  summarized  in 
figure  34. 

(b)  Mode  Number  2 

Natural   frequency  and  mode  shape  information  for  this  mode  is  obtained 
from  [i-i ,  Up,  and  U^.     This  mode  is  excited  primarily  in  the  early  part  of 
the  earthquake  response,  and  so  information  is  most  accurately  obtained  from 
the  expanded  acceleration  records  shown  in  figure  28. 

Frequency  information  is  obtained  from  U-,  from  t  =  2.5  to  5  seconds. 
Seven  cycles  of  motion  are  observed  in  2.386  seconds  yielding  u^^  ^  18.434 
radians/sec. 

Mode  shape  information  is  obtained  from  the  average  of  four  peaks  from 
3  to  4  seconds  into  the  record.     It  is  assumed  that  the  motion  is 
essentially  free  vibration.     For  U-,   scaled  to  1.000,  the  relative  amplitudes 
of  Up  and  U^  in  this  mode  of  vibration  are  -2.005  and  2.130  respectively. 
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Mode  2  characteristics  are  summarized  in  figure  35. 

(c)  flode  Number  3 

Mode  3  information  is  obtained  from  the  lateral  accelerations  U^  and 
Ug.  Because  the  bridge  is  supported  by  a  single  column,  it  is  highly  likely 
that  there  should  exist  significant  torsional  motion,  but  none  is  observed 
at  this  lateral  mode  frequency.  Mode  shape  data  is  taken  from  t  =  2.5  to 
4.0  seconds  even  though  significant  lateral  ground  motion  exists  in  this 
interval.  The  ground  motion  during  this  time  has  a  predominant  frequency 
significantly  different  from  this  natural  frequency.  Therefore,  mode  shape 
information  will  be  obtained  from  averaging  many  response  peaks,  assuming 
free  vibration  behavior. 

This  mode  of  vibration  is  the  fundamental,  symmetric,  lateral  flexural 
mode  of  vibration.  This  behavior  is  anticipated,  as  there  exists  only  one 
side-span  instrumented  with  a  lateral  motion  accelerometer.  The  symmetric 
nature  of  this  mode  is  assumed. 

From  five  cycles  of  motion  from  t  =  2.5  to  4.0  seconds,  from  both  U^ 
and  U5,  the  natural  frequency  was  measured  to  be  0J3  =  23.957  rps.  At  this 
natural  frequency  both  U4  and  Ug  were  in  phase,  and  an  average  of  eleven 
peaks  of  motion  yielded  the  relative  amplitudes  (U4  arbitrarily  scaled  to 
unity)  of  U4  =  1.000  and  Ug  =  0.862.  Mode  3  characteristics  are  summarized 
in  figure  36. 

(d)  Mode  Number  4 

flode  4  was  anticipated  and  is  a  symmetric,  vertical  flexural  mode  of 
the  deck.  Motion  in  the  mode  is  observed  in  U2  and  U3  only  for  t<  2 
seconds.  Evidently,  a  small  vertical  earthquake  motion  preceded  the 
significant  ground  motion  and  excited,  to  a  small  extent,  this  symmetric 
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vertical  mode  of  vibration.  Undoubtedly  motion  in  this  mode  exists 
throughout  U2  and  U^  but  is  overshadov/ed  by  the  large  vertical  components 
of  Modes  1  and  2. 

Frequency  and  mode  shape  information  is  obtained  from  three  cycles  of 
motion  before  t  =  2  sec.  The  natural  frequency  was  found  to  be  W4  =  17.527 
rps.  Both  U2  and  U3  are  in  phase  in  this  mode  and  are  approximately  equal 
in  amplitude.  No  more  accurate  quantitative  measurements  are  justified  from 
the  data.  Mode  4  characteristics  are  summarized  in  figure  37. 

(e)  Mode  Number  5 

Mode  5  is  a  torsional  mode  of  vibration  observed  in  the  eccentrically 
located  vertical  accelerations  U5  and  U7.  Such  motions  have  \/ery   high 
accelerations,  but  because  of  the  high  natural  frequency  of  this  mode, 
torsional  displacements  may  not  be  large.  Frequency  and  mode  shape 
information  is  obtained  for  t  =  2.7  to  4.0  seconds.  The  torsional  component 
is  obtained  as  {U^-l}     )/c   and  (U7-U3)/c  where  c  is  half  of  the  deck  width. 
For  t  >  2.5  seconds  the  motion  is  essentially  free,  so  U^q  can  be  ignored 
relative  to  U5,  but  U3  must  still  be  removed  from  U7. 

Frequency  data  are  obtained  from  nine  cycles  of  motion  between  t  =  2.7 
and  4.0  seconds  from  [J5  and  U7,  and  are  summarized  in  table  8.  Mode 
shape  data  were  taken  from  the  average  of  eighteen  peak  values  after  the 
vertical  motion  U3  had  been  subtracted  from  U7  (see  dashed  line  on 
U7  trace  on  figure  30).   Torsional  motion  at  mid-side-span  is  in  phase  with 
torsional  motion  at  the  center  of  the  bridge.  Although  motion  on  one  side 
only  was  recorded,  it  is  assumed  that  this  mode  is  symmetric.  Furthermore, 
although  it  is  not  measurable,  it  is  assumed  that  there  must  exist  some 
lateral  component  to  this  mode  (it  is  also  assumed  that  there  must  exist 
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Instru. 

No.  of  Cycles 

Time  (sec) 

w] (rad/sec) 

2 
3 

8 
8 

4.381 

4.427 

Average 

11.474 
11.354 

=  11.414 

Table  7.-  Continuous  Two-Span  Bridge,  Mode  1 
Frequency  Surraiiary 


Instru. 

No.  of  Cycles 

Time 

W5(rad/sec) 

5 

7 

9 
9 

1.283 
1.292 
Average 

44.075 
43.768 

=  43.922 

Table  8.-  Continuous  Two-Span  Bridge,  Mode  5 
Frequency  Summary 
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wj  =  11.414  rps 


Figure  34.-  Continuous  two-span  bridge,  mode  1  characteristics. 
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tOg  =  18.434  rps 


Figure  35.-  Continuous  two-span  bridge,  mode  2  characteristics. 
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CO  =  23.957  rps 


Figure  36.-  Continuous  two-span  bridge,  mode  3  characteristics. 
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some  torsional  component  to  the  lateral  mode).  For  a  rotation  at  the  center 
of  the  bridge  of  1,000,  the  rotation  at  mid-side-span  is  measured  to  be 
0.794.  Mode  5  characteristics  are  summarized  in  figure  38. 

(f)  Discussion 

It  is  of  value  to  compare  the  structure  model  parameters  which  have 
been  found  from  this  hand  analysis  to  those  known  parameters  which  were  used 
to  generate  these  hypothetical  records.  A  summary  of  the  actual  structure 
model  parameters  are  shown  in  figure  39.  As  can  be  seen  the  natural 
frequencies  were  found  quite  accurately.  On  the  other  hand,  mode  shapes 
were  found  only  in  a  qualitative  sense,  with  numerical  accuracy  for  a  few 
components.  In  order  that  mode  shapes  be  found  with  any  quantitative 
accuracy,  even  for  this  simple  example,  the  records  must  be  analyzed  on  a 
digital  computer. 

Note  that  damping  estimates  were  not  even  attempted  for  this  example. 
In  certain  cases  numerical  accuracy  may  be  achieved.  In  other  cases, 
damping  estimates  may  be  highly  erroneous.  Because  of  the  great 
inconsistency  which  may  exist  in  damping  estimates  obtained  from  hand 
analyses  (except  for  the  simplest  of  examples  like  the  triple  simple  span 
bridge)  it  is  not  recommended  to  find  damping  ratios  from  hand  analyses. 

These  two  examples,  as  well  as  one  additional  example,  will  be 
reanalyzed  with  the  aid  of  a  high-speed  digital  computer  in  the  next 
section. 

One  last  comment  is  in  order.  Note  that  at  most,  mode  shapes  may  be 
found  with  numerical  accuracy  at  a  few  points  where  instruments  exist. 
Deformations  between  instruments  can  only  be  assumed.  Unless  bridges  are 
yery   heavily  instrumented,  mode  shapes  are  not  likely  to  be  found  with  great 
accuracy.  The  results  from  these  analyses  are  best  used  to  verify 
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(0^  =  17.527  rps 
Figure  37.-  Continuous  two-span  bridge,  mode  4  characteristics. 


oog  =  43.914  rps 


Figure  38.-  Continuous  two-span  bridge,  mode  5  characteristics, 
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Figure  39.-  Continuous  two-span  bridge,  actual  modal  characteristics, 
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theoretical  behavior  of  a  few  points.  Experimentally  obtained  structure 
models  will  rarely  have  sufficient  detail  to  be  used  by  themselves  for 
complete  earthquake  analyses. 
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IV.     ANALYSIS  OF  DIGITIZED  STRONG-MOTION  RECORDS  WITH  THE 
ASSISTANCE  OF  A  HIGH-SPEED  COMPUTER 

A.     Introduction 

The  objectives  of  this  section  are  the  same  as  those  for  the  previous 
section:     to  compute  force  levels  throughout  the  bridge  for  a  given 
earthquake,  and  to  deduce  what  mathematical  model   would  have  best 
duplicated  the  measured  response  to  the  measured  ground  motion.     Again, 
this  second  objective  is  the  more  important. 

Once  deflected  shapes  and  corresponding  coordinate  amplitudes  have 
been  obtained  by  any  means,  force  levels  throughout  the  bridge  may  be 
computed  as  outlined  in  the  previous  section.     Again  those  force  levels  may 
be  computed  from  assumed  shapes  ip^-(s)  and  corresponding  coordinate 
accelerations  a-(t),  or  computed  normal   mode  shapes  i-(s)  and  corresponding 
modal   accelerations  ri^(t).     The  remainder  of  this  section  will   be  devoted 
to  the  second  objective;  to  identify  mathematical   model   properties  of  the 
bridge  from  the  measured  strong-motion  records. 

Identification  of  a  best-fit  mathematical  model   from  strong-motion 
records  (called  system  or  parameter  identification)   is  still   yery  much  in 
the  developmental   stage,  particularly  for  bridges  because  their  natural 
frequencies  of  vibration  are  so  often  closely  spaced.     Traditionally, 
mathematical   models  of  the  dynamic  bridge  behavior  are  formulated, 
theoretical   bridge  responses  are  compared  to  the  measured  responses,  and 
finally  the  theoretical  model   is  evaluated  by  comparing  theoretical   and 
measured  responses.     If  the  fit  is  good,  nothing  more  need  be  done,  but  if 
the  fit  is  poor,  many  or  all   of  the  model   parameters  must  be  changed  in 
hopes  of  improving  the  fit.     Unfortunately,  even  for  the  simplest  bridge,  a 
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mathematical  model  may  include  hundreds  if  not  thousands  of  parameters  that 
may  be  varied.  Since  the  possible  combinations  of  changed  parameters  is 
enormous,  system  identification  by  this  unsystematic  trial  and  error 
procedure  is  extremely  tedious  if  not  entirely  impractical.  Since  it  is 
not  within  the  scope  of  this  report  to  discuss  mathematical  modeling 
(Biggs,  1964;  Clough  and  Penzien,  1975,  Hurty  and  Rubinstein,  1964),  system 
identification  by  modifying  an  a  priori  model  will  not  be  considered.  In 
any  event,  this  approach  has  not  yet  been  demonstrated  to  be  the  best; 
therefore  in  this  report  no  a  priori  model  will  be  formulated. 

There  are  two  fundamental  methods  of  identifying  the  mathematic  model 
from  strong-motion  records  when  no  a  priori  model  is  formulated:  frequency 
domain  methods  and  time  domain  methods.  Traditionally  in  the  fields  of 
electrical  and  mechanical  engineering,  the  frequency  domain  method  has  been 
used.  Specifically,  Fourier  transforms  (Clough  and  Penzien,  1975;  Hurty 
and  Rubinstein,  1964)  of  the  ground  motion  and  the  bridge  motion  are 
computed.  The  bridge  motion  transform  is  then  divided  by  the  ground-motion 
transform  to  yield  a  complex  transfer  function,  H(oj),  or  frequency  response 
function.  The  amplitude  of  a  typical  transfer  function  is  shown  in  figure 
40.  The  transfer  function  amplitude  between  a  point  on  the  bridge  and  the 
ground  is  the  ratio  of  peak  bridge  response  to  peak  ground  motion  for  a 
steady  harmonic  ground  motion  with  constant  frequency,  co.  For  input  motion 
frequencies  near  a  natural  frequency  of  the  bridge,  resonance  will  occur 
causing  large  modal  motions.  For  input  motion  frequencies  away  from 
natural  frequencies  the  bridge  may  act  as  a  vibration  isolator  and  filter 
out  those  input  motions  entirely. 

From  a  transfer  function,  natural  frequencies  can  be  found  as  the 
locations  of  the  dominant  peaks;  damping  ratios  can  be  found. from  the  peak 
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Figure  40.-  Transfer  function  amplitude. 
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widths;  and  mode  shapes  can  be  found  from  the  relative  record  amplitudes  at 
a  specific  natural  frequency  (Clough  and  Penzien,  1975;  Hurty  and 
Rubinstein,  1964).  Once  mode  shapes  have  been  found,  participation  factors 
can  be  computed  from  equations  60,  61,  and  62  from  the  bridge  geometry  and 
mass  distribution,  M(s). 

Although  this  procedure  sounds  wery   simple  once  computer  subroutines 
are  available  to  compute  Fourier  transforms,  it  is  not  particularly 
accurate  when  used  to  analyze  earthquake  induced  bridge  motions.  Inherent 
in  any  frequency  domain  method  using  Fourier  transforms  is  the  assumption 
that  the  input  and  output  motions  are  periodic  indefinitely;  namely  motion 
begins  at  t=-°°  and  continues  to  t=°°.  Obviously  earthquakes  are  not 
periodic,  and  most  cannot  even  be  approximated  well  as  being  periodic.  The 
assumption  of  periodicity  in  the  response  implies  that  some  sort  of 
steady-state  response  exists.  In  actuality,  by  the  time  that  significant 
ground  motions  have  ceased,  bridge  motions  may  have  just  begun  to  grow 
towards  a  steady-state  response.  Consequently,  when  a  transfer  function 
between  input  and  output  motions  is  found,  peaks  at  natural  frequencies  are 
deceptively  small.  As  a  result  modal  peaks  may  not  even  be  identifiable 
against  background  "noise",  and  when  identifiable  may  yield  erroneous 
values  of  damping. 

Because  of  the  many  uncertainties  in  the  results  obtained  from 
frequency  domain  procedures,  they  are  not  recommended  for  general  use  in 
analyzing  strong-motion  records  for  highway  bridges.  In  certain 
circumstances,  natural  frequencies  may  be  identified  from  Fourier 
transforms  of  absolute  bridge  motions,  but  not  reliably.  Other 
publications  (Biggs,  1964;  Clough  and  Penzien,  1975;  Hurty  and  Rubinstein, 
1964;  Newmark  and  Rosenblueth,  1971;  Wiegel ,  1970)  have  further  information 
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regarding  frequency  domain  analysis  procedures. 

Most  reliable  strong -mot ion  record  analysis  procedures  seem  to  be  time 
domain  procedures.  Recent  summaries  of  system  identification  procedures 
are  given  in  Hart  and  Yao  (1976)  and  Ibanez  (1977).  The  methods  described 
were  developed  to  analyze  building  structures,  and  may  be  used  with  only 
some  success  for  highway  bridges.  Consequently,  a  time  domain  analysis 
procedure  has  been  developed  to  analyze  highway  bridge  motions  from  a 
general  least-squares  fitting  procedure  described  by  Kozin  and  Kozin 
(1968.) 


b.  Time  Domain  Analysis  Procedure 

A  time  domain  analysis  procedure  which  has  been  developed  primarily 
for  the  analysis  of  highway  bridges  will  be  described  here.  Although 
highway  bridges  are  simple  in  outward  appearance,  their  dynamic  behavior  is 
in  general  very   complex  because  natural  frequencies  of  vibration  may  be 
closely  spaced,  and  modes  of  vibration  may  be  excited  by  all  components  of 
the  ground  motion.  The  procedure  described  here  is  unaffected  by  these 
complexities. 

The  objective  of  this  system  identification  analysis  is  to  describe 
completely,  from  the  recorded  motions,  the  mathematical  model  that 
duplicates  the  actual  motion  best.  The  mathematical  model  to  be  identified 
can  be  either  in  physical  coordinates  associated  with  lumped  masses 

[M]{U}+[C]{U}+  [K]{U}  =  -{f\>%-{My>Uyg-^:Mz^U^g       (90) 

or  in  generalized  coordinates  associated  with  a  set  of  assumed  distorted 
shapes 
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[m]{a}+[c]{a}+[l<]{a}  =  -^n^>Uxg"("iy>UygOn^}U^g         (91) 

In  either  case,  the  ground  accelerations  U  ,  U  .  and  U   are  recorded  as 

'^y     jy  ^y 

such  and  need  no  further  processing  before  analysis. 

The  following  paragraphs  describe  the  preprocessing  required  to 
identify  the  mathematical  model  in  physical  coordinates  described  by 
equation  90. 

Relative  motions,  or  distortions,  {u},   {u},   {u}  must  first  be  computed 
from  the  set  of  recorded  absolute  bridge  accelerations  {U}.     A  bridge 
acceleration  relative  to  the  ground,   u-(t),  is  found  simply  as  the 
recorded  bridge  acceleration,  less  the  component  of  the  recorded  ground 
acceleration  in  the  direction  of  that  transducer.     Relative  bridge 
velocities,  uAt) ,  and  displacements,  u^(t),  must  then  be  computed  from  the 
relative  acceleration.     These  time  histories  may  be  computed  by  any 
technique,  but  typically  the  accelerations  are  first  digitized  and  assumed 
to  be  linear  between  points.     These  piecewise  linear  functions  are  then 
integrated  numerically  to  yield  the  corresponding  velocities  and 
displacements  (Conte,  1965).     Often,  after  one  and  particularly  after  two 
integrations,  low-frequency  noise  becomes  obvious.     For  ground  motions  with 
frequencies  that  are  low  with  respect  to  the  lowest  natural   frequency  of 
the  bridge,  the  bridge  will   respond  essentially  as  a  rigid  body  with 
negligible  distortions.     Consequently,  large-amplitude  and  low-frequency 
components  in  the  relative  velocities  and  displacements  must  be  errors  and 
can  therefore  justifiably  be  filtered  from  the  time  histories,     A 
reasonable  lower  limit  to  the  frequencies  passed  by  the  filter  is  wi/4, 
where  w-,   is  the  estimated  lowest  natural   frequency  for  the  bridge.     The 
time  histories  u.(t),  u-(t),  and  ii. (t)  are  now  in  suitable  form  for 
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analysis. 

It  is  assumed  that  the  bridge  geometry  and  mass  distribution  can  be 

obtained  accurately  from  material  weights  and  the  bridge  drawings.  The 

rule  for  lumping  masses  differs  for  each  bridge,  and  is  demonstrated 

graphically  in  the  example  analyses  which  follow.  The  participation  factor 

M  ■  simply  equals  the  corresponding  mass  M.  if  the  ith  transducer  is 

directed  parallel  to  U   or  equals  zero  if  it  is  perpendicular  to  U  .  The 

xg  xg 

participation  factors  M  .  and  M  .  are  similarly  defined.  If  the  ith  mass 

y  I       z  I 

is  a  rotational  inertia,  then  all  corresponding  participation  factors  M  . , 
M  •,  and  M  •  equal  zero. 

Equation  90  can  now  be  written  in  a  more  convenient  form  for  analysis: 

[C]{u}+[K]{u}  =  {B}  (92) 

where  the  vector  {B}  is  a  vector  of  known  time  histories  defined  by 

{B}  =  -[M]{U}-{M/u^g-{M^}U^^-{M^}U^g  (93) 

and  the  damping  coefficients,  C--,  and  the  stiffness  coefficients,  K--,  are 
the  unknowns  to  be  determined  from  this  identification  procedure. 

This  completes  the  preprocessing  required  to  identify  the  mathematical 
model  in  physical  coordinates  (equation  90).  The  following  paragraphs 
describe  the  required  preprocessing  to  identify  the  mathematical  model  in 
generalized  coordinates  (equation  91), 

In  order  to  compute  the  appropriate  masses  and  time  histories  in 
generalized  coordinates,  assumed  distorted  shapes,  ii;-(s),  must  be  selected. 
One  distorted  shape  is  allowed  for  each  independent  bridge  motion  recorded. 
Although  the  success  of  this  analysis  procedure  does  not  depend  upon  it, 
the  accuracy  of  the  analysis  is  improved  if  the  first  n  assumed  shapes 
if,-(s),  i=l,  ...»  n.  most  closely  approximate  the  first  n  most  significant 
normal  mode  shapes  $.(s)  i=l,  ...,  n.  Time  histories  associated  with  the 
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assumed  distorted  shapes  are  computed  as  follows. 

Let  ^p^-   be  the  amplitude  of  the  jth  assumed  shape  at  the  ith 
transducer  location  in  the  ith  transducer  direction.  There  will  exist  a 
coordinate  transformation  between  the  vector  of  relative  bridge 
accelerations,  {u(t)},  defined  earlier,  and  the  vector  of  generalized 
accelerations,  {a(t)},  given  by 

{u(t)}  =  [ti;]{a(t)}  (94) 

As  the  assumed  distorted  shapes  must  be  linearly  independent,  the  inverse 
transformation  exists: 

{a(t)}  =  [i!;r''{u(t)}  (95)  . 

yielding  the  generalized  accelerations  immediately.  The  generalized 
velocities,  {a(t)},  and  generalized  displacements,  {a(t)},  may  be  computed 
from  the  generalized  accelerations,  {a(t)},  by  the  numerical  integration 
procedure  outlined  earlier. 

Given  the  mass  distribution,  M(s),  and  the  assumed  shapes,  $.(s),  the 
masses  and  participation  factors  needed  for  analysis  are  computed  again  by 

niij  =  /M(s);i;.(s).;i;j.(s)ds  (96) 


m 


XI 


,-  =  /M(s);i;.(s).ids  (97) 


myi  =  /M(s);i;.(s).jds  (98) 

m^i  =  /M(s);i;.(s).kds  (99) 

Equation  91  can  now  be  written  in  the  more  convenient  form  for 
analysis  (similar  to  equation  92): 

[c]{a}+[k]{a}  =  {b}  (100) 
where  {b}  is  another  vector  of  known  time  histories 

{b}  =  -[m]{a}-{m/u^g-{my}Uyg-{m^}U^g  (101) 

and  the  damping  coefficients,  c..,  and  the  stiffness  coefficients,  k-., 

'J  J 
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are  the  unknowns  in  generalized  coordinates  to  be  identified  from  this 
identification  procedure. 

This  completes  the  preprocessing  required  to  identify  the  mathematical 
model   in  generalized  coordinates   (equation  91).     The  procedure  to  identify 
the  coefficients  in  equation  92  is  obviously  identical   to  the  procedure 
required  to  identify  the  unknown  coefficients  in  equation  lOO.     The  latter 
equation,  in  generalized  coordinates,  will  be  used  in  the  description  of 
this  procedure. 

The  unknown  coefficients  are  found  by  a  linear  least-squares  fitting 
procedure  in  the  time  domain.     Specifically,  the  left-hand  side  of  equation 
(100)   is  fitted  in  a  least-squares  sense  to  the  right-hand  side  as  a 
function  of  the  unknown  coefficients  c • .  and  k-.. 

Equation  (100)   is  a  system  of  n  equations  linear  with  respect  to  the 
unknown  coefficients.     The  total   squared  error,  E,  is  defined  to  be  the  sum 
of  the  squared  errors  between  left-hand  and  right-hand  sides  of  all 
equations: 

1=1      t-j 

where  the  integral  is  over  a  time  interval  of  any  size,  located  anywhere  in 

the  total  record,  i.e.,  the  initial  conditions  for  the  time  interval 

analyzed  are  not  restricted  to  be  zero. 

The  squared  error  is  minimum  or  maximum  for  c^-  •  and  k^  •  that  satisfy 

8E/3C..  =  0  (103) 

u 

9E/8k^.j  =  0  (104) 

for  all  independent  c^  •  and  k... 
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Equations  (103)  and  (104)  when  expanded  yield  a  system  of  linear 
equations  with  respect  to  the  unknowns  c-  •  and  k...  The  exact  form  of  the 
equations  depends  upon  the  choice  of  independent  unknowns.  In  general,  it 
is  assumed  that  c-.  =  c-  and  k.  .  =  k...  These  off-diagonal  terms  may  be 
evaluated  independently,  however,  or  other  dependence  amongst  the  unknowns 
may  be  assumed,  if  justified.  The  resulting  system  of  linear  equations  may 
then  be  readily  solved  for  the  unknown  coefficients  c,--  and  k-.  using  any 
accepted  procedure  (Conte,  1965;  Gere  and  Weaver,  1965). 

The  actual  fitting  procedure  is  performed  on  a  high-speed  digital 
computer.  Therefore,  recorded  motions  are  analyzed  as  time  series. 
Integrals,  which  result  when  equations  103  and  104  are  expanded,  are 
replaced  by  summations.  Listings  and  descriptions  of  all  computer  programs 
necessary  to  perform  the  various  phases  of  this  analysis  are  included  in 
Appendix  C.  Use  of  these  programs  is  demonstrated  in  the  example  analyses. 

C.  Transformation  to  Normal  Coordinates 

Once  the  records  are  analyzed,  the  best-fit  linear  mathematical  model 
of  the  dynamic  bridge  behavior,  either  equation  90  in  physical  coordinates 
or  equation  91  in  generalized  coordinates,  is  completely  defined.  However, 
as  discussed  in  Section  II,  this  mathematical  model  is  not  in  its  most 
useful  form.  The  model  is  often  most  informative  when  expressed  in  normal 
coordinates.  By  making  this  transformation  no  new  information  is  obtained, 
but  what  has  been  obtained  is  simply  changed  into  a  more  useful  form.  The 
transformation  from  generalized  to  normal  coordinates  will  be  used  to 
demonstrate  the  procedures;  the  transformation  from  physical  to  normal 
coordinates  is  identical. 

Related  to  equation  91  is  the  eigenvalue  problem  (Biggs,  1964;  Clough 
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and  Penzien  1975;  Hurty  and  Rubinstein,  1964) 

[[k]-a)^[m]]{e}  ={0}  (105) 

where  w     (the  natural   frequencies  squared)  are  the  eigenvalues  and  {9}  are 
the  corresponding  eigenvectors  (natural  mode  shapes  of  the  bridge  in  terms 
of  generalized  coordinates).     Equation  105  will   have  a  solution  (non-unique) 
if  and  only  if  the  determinant 

|[k]-a)^[m]|=  0  (106) 

Equation  106  when  expanded  is  an  nth  order  polynomial  which  can  be  solved 
for  the  unknown  natural   frequencies,  oo.     For  each  natural   frequency  w- 
found,  equation  105  can  then  be  solved  for  a  corresponding  mode  shape  {9^-}. 
Since  the  solution  to  equation  105  is  non-unique  on  satisfying  equation 
106,  the  mode  shapes  can  be  found  only  in  terms  of  some  arbitrary  constant. 
Often  the  maximum  term  in  {9-}  is  arbitrarily  set  to  unity,  but  any 
normalization  scheme  is  acceptable. 

The  mode  shapes  found  in  this  manner  have  the  desirable  property  of 
di agonal i zing  both  the  mass  and  stiffness  matrices.     Specifically,  for  the 
mode  shapes  arranged  as  columns  in  the  transformation  matrix,   [9], 

[ef  [m][9]  =  [M]  (107) 

[ef  [k][9]  =  [K]  (108) 

where  both  [M]  and  [K]  are  diagonal   matrices. 

The  matrix  [9]  is  in  fact  the  desired  transformation  matrix  defining 
the  transformation  from  generalized  to  normal  coordinates.     For  the 
transformation  defined  by 

{a(t)}  =   [9]{n(t)}  (109) 

equation  91  is  transformed  to 

[M]{n}+[C]{n}+[/<]{n}  =  -^VV^^^V^^^z^'^zg  ^™^ 

where  [M]  and  [K]  are  defined  by  equations  107  and  108  and 
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[c]  =  [ef[c][e]  (111) 

{M^}  =  [ef{m^}  (112) 

{My}  =  [efim^}  (113) 

{M^}  =  [ef  {mp  (114) 

The  mass  and  stiffnesses  are  diagonal   by  definition  of  this 
trans fomiati on.     The  damping  matrix  [C]  may  or  may  not  be  diagonal i zed. 
Often  this  determination  of  the  character  of  [C]  is  in  itself  an  objective 
of  the  analysis,  whether  or  not  a  specific  modal   damping  ratio  does  exist. 
Again  if  it  is  assumed  that  [C]  is  diagonalized,  or  if  the  off-diagonal 
terms  are  small   and  can  be  ignored,  then  equation  110  premulti plied  by 
[M]"    yields  n  uncoupled  modal  equations  of  motion 

where  ? • ,  w • ,  y   ■,  y   • ,  and  y   .   are  defined  from  equation  17. 

It  is  to  be  noted  that  the  mode  shapes,  {6}-,  are  found  in  terms  of 
generalized  coordinates.     It  is  more  useful   to  know  the  mode  shapes,  $^-(s), 
in  physical   coordinates.     Specifically  these  are  given  by 

5^.(s)  =  e^.5^(s)+e2^.$2(^^'^'"'^^ni^n^^^  '         ^^^^^ 

Had  the  transformation  been  from  physical  to  normal  coordinates  directly, 
the  mode  shapes  {6-}  would  have  been  in  physical  coordinates  already,  and 
no  further  expansion  with  equation  116  would  have  been  necessary. 

This  concludes  the  description  of  how  strong-motion  records  may  be 
analyzed  with  the  aid  of  a  high-speed  digital  computer.     Many  extensions  to 
this  analysis  method  are  possible,  including  the  assumption  of  nonlinear 
behavior,  but  a  description  of  these  extensions  is  beyond  the  scope  of  this 
manual . 
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D.     Example  Analyses 

In  all   of  the  example  analyses  which  follow,  the  motions  of  the  highway 
bridges  to  an  actual  earthquake  are  simulated  mathematically  on  a 
high-speed,  digital   computer.     For  all  examples,  the  earthquake  ground 
motions   used  were  the  three  ground  accelerations  recorded  at  D,  Johnson 
Ranch,  Oroville,  Calif.,  for  the  aftershock  recorded  6  August  1975 
(Sherburne  and  Hauge,  1975).     Plots  of  the  three  recorded  ground 
accelerations  are  shown  in  figure  13. 

As  a  first  step  in  the  development  of  the  simulated  responses, 
lumped-mass  models  of  all  examples  were  made  and  their  natural   frequencies, 
natural  mode  shapes,  and  participation  factors  were  obtained  using  the 
structural  analysis  program  SAP  IV  (Bathe,  Wilson,  and  Peterson,  1973). 
Modal   damping  ratios  were  assumed,  and  simulated  bridge  accelerations  were 
computed  at  hypothetical   instrument  locations  using  the  modal   superposition 
technique  (Biggs,  1964;  Clough  and  Penzien,  1975;  Hurty  and  Rubinstein, 
1964).     These  absolute  bridge  accelerations  were  then  used  in  the  analyses 
which  follow  as  though  they  had  been  recorded  by  strong-motion  instruments o 

Of  the  three  examples  which  follow,  the  continuous  two-span  bridge  and 
the  triple  simple  span  bridge  are  the  same  bridges  subjected  to  the  same 
ground  motions  as  the  two  examples  in  Section  III.     The  three-span 
continuous  bridge  example  is  introduced  in  this  section. 

Two-Span  Continuous  Bridge  -  The  subject  bridge  is  a  simple  continuous 
two-span  highway  grade  separation  bridge  with  geometry  and  section 
properties  as  shown  in  figure  17.  It  is  a  typical   reinforced  concrete  box 
girder  bridge  supported  by  a  single  central   column.     Both  ends  of  the  bridge 
deck  are  restrained  laterally  and  free  to  move  longitudinally.     The  column 
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is  fixed  at  its  base  and  rigidly  cast  into  a  stiff  transverse  beam  in  the 
bridge  deck. 

The  instrumentation  required  to  identify  the  significant  distortions  of 
the  bridge  is  shown  in  figure  41.  This  instrumentation  is  minimal  and  is 
expected  to  be  in  keeping  with  what  has  been  used  or  what  should  at  least  be 
used  on  such  a  simple  bridge.   Again  it  should  be  noted  that  the  abutment 
and  column  base  motions  have  been  assumed  to  be  identical.  Furthermore,  no 
instrumentation  is  provided  to  measure  relative  bridge-abutment  motions  at 
either  of  the  two  ends.  The  instrumentation  shown  is  designed  principally 
for  the  identification  of  the  best-fit  mathematical  model.  The  mathematical 
model  of  the  bridge  motion  to  be  identified  cannot  be  expected  to  be  more 
than  approximate  without  additional  instrumentation.  The  expected 
distortions  recorded  by  this  instrumentation  are  the  following:  (1) 
longitudinal  deck  motionj  (2)  antisymmetric  vertical  flexure;  (3)  symmetric 
vertical  flexure;  (4)  lateral  flexural  motion;  and  (5)  longitudinal 
torsional  motion.  There  is  one,  and  only  one,  instrument  to  record  each  of 
these  distortions.  Since  the  spanwise  distributions  of  these  five 
distortions  are  not  measured,  they  must  be  assumed.  For  this  reason  the 
results  expected  from  this  analysis  will  be  approximate,  unless  of  course 
the  assumed  spanwise  distributions  of  the  distortions  happen  to  be  exact. 

Before  analyzing  these  bridge  motions,  the  data  must  be  processed  as 
described  in  the  previous  section.  Let  Ui(i=l ,2,3,4,5)  be  the  five  absolute 
accelerations  recorded  on  the  bridge. 

The  absolute  accelerations  must  first  be  transformed  into  accelerations 
relative  to  the  ground.  If  all  bridge  instruments  are  aligned  with  at  least 
one  ground  instrument,  this  transformation  is  achieved  simply  through 
subtraction.  The  relative  bridge  accelerations,  Uj  (i=l ,2,3,4,5),  are  then 
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Because  masses  are  difficult  to  estimate  for  coordinates  corresponding 
to  these  instrument  motions,  it  is  advantageous,  as  described  in  earlier 
sections,  to  transform  these  relative  accelerations  to  generalized 
accelerations  in  coordinates  corresponding  to  the  five  principal 
distortions.     If  T,  j,  k,  are  unit  vectors  corresponding  to  the  IJ     ,  l)     , 
and  U      accelerations  respectively,  the  five  deck  distortions  are  assumed  to 
be 

longitudinal  motion  'i^l(x)  =  i 

Antisymmetric  Vertical   Flexure      H>oi>^)  =  [sin(7rx/150)]k 
Symmetric  Vertical   Flexure  ii^3(x)  =  [sin(7Tx/(300)+sin(7Tx/100)]k 

^4(x).  =  [sin(7Tx/300)]j  > 

i^{x)  =  [Sin(7Tx/300)]i 
The  first  four  of  these  distorted  shapes  are  shown  in  figure  23  and  the 
fifth  has  the  identical  distribution  with  the  fourth.     For  a-(t)  equal   to 
the  generalized  accelerations  corresponding  to  these  assumed  distortions, 
the  transformation  from  relative  accelerations  to  generalized  accelerations 
is  given  by  equation  122.     The  term  Tprr  equals  c*  sin(1507T/300) ,  where 
c  equals  18  feet,  half  the  deck  width. 


Lateral   Flexure 
Longitudinal  Torsion 
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Equation  122  is  equation  94  expanded  for  this  example.     The  generalized 
accelerations  are  readily  obtained  from  the  relative  accelerations  by 
inverting  the  transformation  above.     Therefore, 

"1.000        0  0  0  0 

0  .500     -.500         0  0 

0  .354       .354         0  0         { 

0  0  0         1.000         0 

0  0  0  0 


< 


a2(t) 
a3(t)  \   = 
'd^it) 
\d^{t)J 


U2(t) 
U3(t)  ) 


(123) 


056  J  \u^{t)J 


The  generalized  accelerations  must  be  integrated  once  to  yield  generalized 
velocities  a-(t)  and  again  to  yield  generalized  displacements  a^-(t).  As 
shown  in  figure  42,  linear  acceleration  approximation  is  made  (acceleration 
is  assumed  to  be  linear  between  points)  in  the  integration  computations. 
Integrating  the  expression  in  figure  42  two  times  yields 

a(t)  =  %-^%^H%+'[-%)^^/^^^  (124) 

a(t)  =  a^+a^t+a^t^/2+(a^^^-a^)t^/6At  (125) 

Evaluating  these  expressions  at  t=At  yields  the  next  increments  in  velocity 
and  displacement 

(126) 


%+i  =  %^(%-nS^^^/2 


%+l  =  %-^V*+(%/3+%-,i/6)At 
At  the  start  of  each  record  ai=0  and  ai=0,  by  assumption. 


(127) 
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Figure  41.-  Continuous  two-span  bridge,  instrumentation  plan. 


a(t)  -  an+ Jt — ^ 


Figure  42.-  Linear  acceleration  approximation. 
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Numerically  integrating  real    records  two  times  often  amplifies 
tremendously  the  low-frequency  errors  in  the  recorded  accelerations.     Since 
all   signals  in  the  generalized  velocities  and  displacements  having 
frequencies  low  with  respect  to  the  lowest  natural   frequency  (<— !   or  so)  must  be 
errors,  they  may  justifiably  be  filtered  from  the  records.     From  the 
relative  accelerations,  the  lowest  natural   frequency  of  this  bridge  is 
estimated  to  be  2,0  cps  (12.566  rps).     Therefore  all   signals  in  the 
velocities  and  displacements  having  frequencies  less  than  0.5  cps 
(3.142  rps)  are  removed.     A  phaseless,  Fast  Fourier  Transform  filter  could 
be  used  to  remove  these  errors  (Taoka,  1976).     Since  the  time  histories 
analyzed  here  are  artificially  generated,  there  are  no  low  frequency  errors 
after  the  numerical   integrations,  and  therefore  the  velocities  and 
displacements  have  not  been  filtered. 

Masses  for  this  bridge  are  calculated  from  the  drawings  and  the  unit 
weight  of  concrete,  taken  as  150  pcf.     Therefore,  they  will   not  be 
considered  as  unknowns,  but  as  knowns.     As  outlined  earlier,  the 
generalized  masses  corresponding  to  the  generalized  coordinate  system  are 
given  by 

m^j  =  /o°%x);i;.(x)-ii;j(x)dx  (128) 

and  the  participation  factors  by 

'"xi   ^  /o^°M(x);f.(x)-idx  (129) 

%i  =  -/■o^°M(x)ii(x)-jdx  (130) 

m^.   =  /^'^%x);i;.(x)-kdx  (131) 

for  assumed  deck  distortions  1  through  4.  The  mass  of  the  column  has  been 
ignored  in  the  calculations  relative  to  the  mass  of  the  deck.  Finally, 
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"^55  ""  -^0     ^(x)%(x)*%(x)dx 


(132) 


Since  ii)r(x)   is  a  rotation,  the  appropriate  inertia  to  use  in  the  above 
integral   is  a  rotational   inertia  per  unit  length  of  bridge,  I(x),     Since 
the  ground  rotations  are  assumed  to  equal   zero,  all   participation  factors 
for  this  mode,  m  n  m  g*  m  5,  will   be  zero.     All   off-diagonal  mass  terms  with 
respect  to  the  torsional  mode,  and  all   three  participation  factors  for  this 
mode  are  zero.     These  generalized  masses  and  participation  factors  expanded 
for  this  example  bridge  are,  for  M(x)=0.3295  kip-sec^/ft^  and  I{x)  =  50.939 


kip-sec  , 
[m] 
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0         49.42         0 
0  0         98.85 
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49.42  0 
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V  0   ; 


i^y}    = 


(    0     ^    {m^} 
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V    0    ) 

The  final   preprocessing  required  before  analysis  is  to  combine  all 
inertial  time  histories  into  a  single  vector 


f  0  ^ 
0 

83. 90) 

0 
V    0    ) 


{b(t)}  =  -[m]{a{t)}-{m^}U^g(t)-{my}Uyg(t)-{m^}U^g(t) 

Specifically, 

b^(t)  =  -98.85a^(t)-98.85U^g(t) 


(133) 
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b2(t)  =  -49.42a2(t) 
b3(t)  =  -98.85a3(t)-83.90U^  (t) 
b^Ct)  =  -49.42a4(t)-62.93Uyg(t) 
bgCt)  =  -7641. ag(t) 
All  input  quantities  are  in  suitable  form  for  analysis,  namely 

^ii^i""  ••*  '"^isV^ii^'i"*"  ••*  ""hs^s  =  h 


^is^'i'"  •••  ^^ssV^is^i""  •••  ^h5%  =  ^5 


(134) 


where  the  c--  and  k^--  (both  assumed  symmetric)  are  the  unknown  damping 
and  stiffness  coefficients  to  be  found. 

Although  the  analysis  may  be  performed  at  this  stage,  numerical  errors 
may  be  reduced  if  motions  known  to  be  decoupled  are  in  fact  decoupled 
before  analysis.  This  is  true  because  even  though  the  motions  are 
decoupled  (motions  in  one  direction  entirely  independent  from  motions  in 
another  direction)  the  actual  time  histories  representing  those  motions 
may  be  correlated.  Any  correlations  existing  between  uncoupled  motions 
will  lead  to  numerical  errors  in  the  analysis. 

This  example  problem  may  justifiably  be  separated  into  three  uncoupled 
problems,  due  to  the  symmetry  in  the  structure.  The  three  uncoupled 
problems  become  the  following:  (1)  longitudinal  motion  and  antisymmetric 
vertical  flexure;  (2)  symmetric  vertical  flexure;  and  (3)  lateral  motion 
and  longitudinal  torsion.  Mathematically  this  division  corresponds  to 
dividing  the  above  system  of  equations,  equation  134,  into  the  following 


110 


sets: 


c^„ai+C220'2^'^l?'^l  ^'^22*^2  ~  '^2 


(135) 


for  longitudinal  motion  and  antisymmetric  vertical   flexure; 

c^oOto+k^^a^  =  ^o  (136) 

for  symmetric  vertical   flexure;  and 


I 

^45V^55^5''''45°^4""'55°^5  =  ^5 


(137) 


for  lateral  motion  and  longitudinal   torsion.     All   other  coefficients  c-. 

and  k--  must  equal   zero  due  to  symmetry  in  the  bridge. 

Subroutine  FIT  was  used  to  identify  the  best-fit  values  of  the  unknown 

damping  coefficients,  c--  and  stiffness  coefficients,  k...     Results  of 

the  analyses,  corresponding  to  equations  135,  136,  and  137  are  as  follows: 

69.6ai-5.48ao+27,900aT-5530ao  =  -98.85aT-98.85U„ 

12  12  I  xg  (^3gj 

-5.48a^+23.4a2-5530a-|+8400a2  =  -49.42a2 
and 

69.8a3+31,300a3  =  -98.85613-83. 90U^  (139) 

and 

49.4a.+338(^(-+39,400a.+369,000at:  =  -49.42a/, -62. 93*d^ 

338a4+l  6, 8000^5+369 ,000a4+l  6, 600  jOOOag  =  -7641  ctg 

Once  again  these  equations  of  motion  may  be  recombined  in  the  more  usual 
form 


[m]{a}+[c]{a}+[k]{a}  =  -{m^}U^^-{m^}[i^^-{m^}[i^^  (141) 


where 
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[m] 


98.85  0  0 
0  49.42  0 
0     0    98.85 


0 

0~ 

0 

0 

0 

0 

49.42 

0 

0 

7641 

[c] 


69.6  -5.48  0     0     0 

-5.48  23.4  0     0     0 

0     0  69.8    0     0 

0     0  0  49.4    338. 

0     0  0  338.  16,800 


and 


[k] 


27,900  -5530    0 
-5530   8400    0 


0       31 ,300 


0 
0 
0 


0       39 ,400     369 ,000 
0     369,000  16,600,000 


W       -\ 


/'98.85^ 
0 

0  ){m  >  =  <; 


f    0    ^ 
0 


f  0  ^ 

0 


0 

V   0   ; 


0     >     {m  }     =  <83.90> 


62.93 

V  0  ; 


0 

V  0  ; 


At  this  point  the  mathematical  model   of  the  bridge  has  been  entirely 
defined  from  the  analysis  of  the  strong-motion  records.     The  model,  however, 
is  not  in  its  most  useful   form.     It  is  much  more  useful   for  comparative 
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purposes  to  transform  these  equations  of  motion  into  uncoupled  equations  of 
motion  in  terms  of  normal  coordinates.  At  this  time  we  do  not  know  if 
such  a  transformation  will  exist,  but  we  can  find  a  coordinate  transfor- 
mation which  will  diagonalize  the  mass  and  stiffness  matrices  (Biggs,  1964; 
Clough  and  Penzien,  1975;  Hurty  and  Rubinstein,  1964);  there  is  no 
guarantee,  however,  that  the  transformation  will  also  diagonalize  the 
damping  matrix. 

A  transformation  matrix  is  found  which  transforms  the  above  equations  of 
motion  from  generalized  coordinates  {a}  to  normal  coordinates  (n) 

{a}  =  [e]{n}  (142) 

by  solving  the  related  eigenvalue  problem  (Biggs,  1964;  Clough  and  Penzien, 
1975;  Hurty  and  Rubinstein,  1964): 

[[k]-a)^[m]]{e}  =  {0}  (143) 

The  transformation  matrix  (matrix  of  eigenvectors  arranged  as  columns) 
found  is 

[6] 


1.00 

1.00 

0 

0 

0 

2.74 

-.730 

0 

0 

0 

0 

0 

.707 

0 

0 

0     0     0    1.00   1.00 

0     0     0  -.0302   .214 
The  equations  of  motion,  equation  141,  may  then  be  transformed  into  normal 
equations  of  motion  in  the  following  manner: 


[m]{a}+[c]{a}+[k]{a}  =  -{mx>%-^  VU^g-{m^}U^g 

[m][e]{fi}+[c][e]{n}+[k][e]{n}  =  -K>  V^^>V^'^z^'jzg 

[ef[m][e]{n}+[ef[c][e]{n}+[ef[k][e]{n}  = 

-[ef  {mxliJxg-Ee]"''  m^  Uyg-CejTlm^lU^g 


(144) 
(145) 

(146) 
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For  the  obvious  definitions  this  last  equation  can  be  rewritten  as 


[M]{n}+[c]{n}+[K]{n}  =  -^M^>%-%>Uyg-^Mz}U-,g 


(147) 


where 


[M]      = 

469.9 

0 

0 

0 

0 

0 

125.2 

0 

0 

0 

0 

0 

49.42 

0 

0 

0 

0 

0 

56.39 

0 

0 

0 

0 

0 

399.3 

[C]       = 

215.3 

11.78 

0 

0 

0 

11.78 

90.07 

0 

0 

0 

0 

0 

34o89 

0 

0 

0 

0 

0 

44.31 

2.930 

0 

0 

0 

2.930 

963.4 

m    = 

60,660 

0 

0 

0 

0 

0 

40,450 

0 

0 

0 

0 

0 

15,650 

0 

0 

0 

0 

0 

32,250 

0 

0 

0 

0 

0 

957,500 

{M,}   = 

^98.85^ 

(M^}  = 

^0     ^ 

{M^}  = 

'    °     ] 

98.85 

0 

0 

< 

0 

>             I 

0 

> 

< 

59.32> 

0 

62.93 

0 

I     0     , 

^2.93, 

<    0    , 

As  can  be  seen,  the  off-diagonal  damping  terms  are  small  but  are  not 
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negligible.     The  off-diagonal   damping  forces,  when  compared  to  inertial   and 
restoring  forces,  are  negligible,  and  so  the  off-diagonal   damping  terms  in 
the  end  will  be  ignored. 

The  system  of  equations  of  motion  in  normal   coordinates  can  be  rewritten 
as  five  uncoupled  single-degree-of-freedom  equations  of  motion  (each 
equation  has  been  divided  by  its  mass): 

ny0.458ri^+129.1n^   =  -0.210U  (148) 

fi2+0.719n2'*"223.1n2  =  -0.790U  (149) 

n3+0.706n3+316.7n3  =  -1.200U^  (150) 

.      n'^+0. 786^4+571.9114  =  -l-llSU  (151) 

fi5+2.413n5+2398.n5  =  -0.158U  (152) 

Each  equation  is  compared  to  the  standard  form  (equation  17) 

mm^^^      '-      -T^D^g-Tyiiyg-Y^iJ^g  (l53) 

where  again 

0)   =   natural  circular  frequency  of  vibration  (radians/second), 

^   =   damping  ratio,  and 

Y   =   participation  factor 

For  the  natural  period,  T=27T/a),  the  final  results  are  found  and  are 
given  in  table  9.  The  natural  modes  of  vibration  are  obtained  from  the 
eigenvectors  (columns  in  the  transformation  matrix  [6]),  and  the  assumed 
shapes  of  distortion,  i-(x).  Specifically,  the  ith  mode  shape  in  physical 
coordinates  is  given  by 

$^(x)  =  e-,.$^(x)+e2^.$2^^^^  •"•  ■^^5i'i^5^^)  (^^^^ 

The  natural  modes  of  vibration  found  from  this  identification  procedure  are 
shown  in  figure  43. 
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Figure  43.-  Continuous  two-span  bridge,  mode  shapes  found  from  analysis, 
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Figure  43  (continued).-  Continuous  two-span  bridge,  mode  shapes  found  from 

analysis. 
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This  completes  the  analysis  for  the  best-fit  mathematical  model  for  this 
bridge  from  the  strong-motion  records  alone.  Now  that  the  modal  properties 
have  been  identified,  forces  in  various  elements  of  the  bridge  may  be  found 
in  exactly  the  same  manner  as  was  described  in  Section  III.  For  comparison, 
the  exact  modal  properties  are  given  in  table  10,  and  the  exact  mode  shapes 
are  shown  In  figure  44,  These  properties  were  used  to  generate  the 
hypothetical  responses  used  to  demonstrate  this  analysis. 

In  comparing  the  tabulated  values  identified  from  the  time  histories  to 
the  actual  values  used  to  create  the  theoretical  time  histories,  one  can 
see  that  the  values  identified  are  quite  acceptable  for  the  first  three 
modes.  On  the  other  hand,  the  values  identified  for  modes  4  and  5  are 
somewhat  in  error.  These  discrepancies  are  due  to  the  large  errors  in  the 
shapes  of  the  assumed  lateral  and  torsional  distortions.  Better  results 
may  have  resulted  had  the  assumed  distortions  more  closely  resembled  the 
actual  modes  of  vibration.  These  differences,  however,  would  never  have 
been  detected  without  further  instrumentation.  An  additional  transverse 
and  vertical  accelerometer  at  midspan  would  most  likely  improve  the 
accuracy  of  the  fourth  and  fifth  modal  properties  found. 

Although  the  fifth  mode  shape  identified  looks  particularly  erroneous,  in 
fact  it  is  quite  reasonable.  Mode  5  is  primarily  a  torsional  mode. 
Although  the  numerical  amplitudes  of  rotation  are  comparable  (actually 
slightly  smaller)  than  the  numerical  values  of  lateral  translation  for  the 
actual  mode  shape,  the  vertical  translations  at  the  bridge  deck  edges  due 
to  the  torsional  displacements  are  in  fact  \/ery   large.  Specifically  for  a 
rotation  of  unity,  the  vertical  displacement  at  the  deck  edges  is  ±18.000. 
A  comparison  of  the  identified  and  actual  mode  5  deck  rotations  is  fairly 
reasonable.  It  can  therefore  be  concluded  that  the  significant  deck 
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Mode 

co(rps) 

T(sec) 

?(%) 

Yx 

Yv 

^z 

1 

11.362 

.553 

2.0 

.210 

0 

0 

2 

17.975 

.350 

2.0 

.790 

0 

0 

3 

17.796 

.353 

2.0 

0 

0 

1.200 

4 

23.914 

.263 

1.6 

0 

1.116 

0 

5 

48.969 

.128 

2.5 

0 

.034 

0 

Table  9.-  Continuous  Two-Span  Bridge  Model, 
Parameters  Identified 


Mode 

a3(rps) 

T(sec) 

m) 

^x 

^y 

^z 

1 

11.410 

.551 

2.0 

.189 

0 

0 

2 

18.340 

.343 

2.0 

.780 

0 

0 

3 

17.770 

.354 

2.0 

0 

0 

1.208 

4 

23.880 

.263 

2.0 

0 

1.062 

0 

5 

46.840 

.134 

2.0 

0 

.027 

0 

Table  10.-  Continuous  Two-Span  Bridge  ,  Exact 
Model  Parameters 
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Figure  44.-  Continuous  two-span  bridge,  actual  mode  shapes. 
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Figure  44  (continued).-  Continuous  two-span  bridge,  actual  mode  shanes 
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displacements,  even  for  mode  5,  have  been  found  reasonably  accurately. 
Three-Span  Simply  Supported  Bridge 

The  subject  bridge  for  this  example  is  a  straight  highway  grade 
separation  bridge  with  three  simple  spans.  The  geometry  and  properties  of 
the  bridge  as  well  as  the  accelerometer  locations  are  shown  in  figure  12. 
The  instruments  are  on  the  pier  caps  and  are  assumed  to  sense  the  motion 
of  the  pinned  ends  of  the  simple  spans.  All  spans  are  assumed  to  be 
restrained  laterally. 

Bridge  accelerations  u-(t)  relative  to  the  ground  are  first  computed 
from  the  recorded  bridge  accelerations  U-(t)  and  the  ground  accelerations. 
Specifically, 

u^(t)  =  U^(t)-U^g(t)  (155) 

U2(t)  =  U2(t)-Uxg(t)  (156) 

U3(t)  =  (i3(t)-Uyg(t)  (157) 

U4(t)  =  U4(t)-U^g(t)  (158) 

Since  masses  corresponding  to  each  recorded  motion  may  easily  be  estimated 
from  the  bridge  plans,  there  is  no  need  to  transform  the  relative 
accelerations  to  a  generalized  coordinate  system.  Therefore  parameters  of 
a  lumped  mass  model  in  physical  coordinates  are  to  be  identified. 

As  described  before,  the  relative  accelerations  u^.(t)  are  now 
numerically  integrated  to  yield  the  relative  bridge  velocities  u^-(t)  and 
relative  displacements  u^-(t).  Again,  if  these  had  been  computed  from 
actually  recorded  bridge  accelerations,  low-frequency  errors  would 
undoubtably  exist  in  the  velocities  and  displacements.  These  errors  must 
be  filtered  from  the  bridge  velocities  and  displacements  as  suggested  in 
the  first  example o  Since  no  low- frequency  errors  exist  in  these 
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artificially  generated  time  histories,  they  have  not  been  filtered;  the 
time  histories  are  already  in  suitable  form  for  analysis. 

In  addition  to  these  time  histories,  bridge  masses  and  participation 
factors  are  needed.     It  is  assumed  that  it  is  sufficiently  accurate  to  lump 
bridge  deck  masses  at  the  pier  caps  corresponding  to  the  accelerometer 
coordinates.     Note  that  lateral   flexural  motion  of  the  individual   spans 
have  been  ignored  relative  to  the  rigid  body  translations.   In  any  event, 
such  flexural   distortions  cannot  be  detected  with  this  instrumentation. 
The  masses  corresponding  to  coordinates  at  the  instrument  locations  have 
been  found  for  this  example  and  are  shown  graphically  in  figure  15  (M^ 
is  equal  to  M^)* 

[M]  =     64.463       0  0  0 

0         26.588       0  0 

0  0         34.163       0 

0  0  0         34.163 

Since  the  masses  are  assumed  to  be  lumped  at  the  degrees  of  freedom,  the 
participation  factors  equal   the  corresponding  masses  and  are  given  by 


{M,}     = 


64.463 
26.588 

0 

0 


i^y}       - 


\ 


0 

0 
34.163 
34.163 


/ 


\ 


J 


There  are  no  vertical   degrees  of  freedom,  and  the  horizontal  degrees  of 
freedom  corresponding  to  each  accelerometer  direction  are  not  excited  by 
vertical   ground  motion;  therefore  all   M  •  equal  zero. 

For  the  time  histories  defined  by 

B^Ct)   =  -64.463[uyt)+U^g(t)] 
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B2(t)  =  -26.588[Li2(t)+U^g(t)] 

B3(t)  =  -34.163[U3(t)+Uyg(t)] 
B^Ct)  =  -34.163[u4(t)+Uyg(t)] 
the  system  of  equations  to  be  analyzed  by  this  method  is 
Cll^l^  •••  +C^4U4+K^^u^+  ...  +K^4U4  =  B^ 


(159) 


Ci4Ui+---     +C^4V'^14^1"'  ••*  ■"^44"4  ^  ^4 
Again,  this  system  of  equations  can  be  simplified  by  anticipating  structural 
behavior.     By  assumption  there  is  no  inertial   coupling  (no  off-diagonal 
mass  terms),  and  by  inspection,  there  will   be  no  stiffness  coupling  (a 
force  applied  in  any  individual  accelerometer  direction  will   produce  a 
distortion  only  in  that  individual  direction).     Furthermore,  because  the 
example  bridge  really  consists  of  three  independent  simple  spans  there  is 
not  likely  to  be  any  significant  damping  coupling.     Therefore,  all 
off-diagonal   damping  and  stiffness  coefficients  can  reasonably  be  assumed 
to  equal  zero,  and  the  system  of  equations  to  be  identified  reduces  to  the 
following  four  uncoupled  si ngle-degree-of- freedom  equations  of  motion: 

C-j^u^+K^^u^   =  B^  (160) 

CpoUp+KppUp  =  Bp  (161) 

C33^3^S3"3  =  ^3  ^^^2) 

C44U4+K44U4  =84  (163) 

Subroutine  FIT  was  used  to  identify  the  unknown  coefficients,  and  yielded, 
after  a  slight  rearrangement  of  terms  in  the  B^-(t), 

64.463uy63.1u^+9650u^   =  -64.463U  (164) 

26.588u2+40.5u2+9650u2  =  -26.588U  (165) 
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34.163U3+45.9U3+38,600U3  =  -34.163U  (166) 

34.163u4+45.9u4+38,600u4  =  -34.163U  (167) 

By  comparing  these  single-degree-of-freedom  equations  of  motion,  equations 
164  to  167,  to  the  standard  form,  equation  153,  the  model   parameters  are 
identified  and  shown  in  table  11.     For  comparison,  the  exact  values  which 
were  used  to  generate  the  artificial   bridge  motions  are  given  in  table  12, 
In  essence,  the  results  are  identical. 


Three-Span  Continuous  Bridge 

The  subject  bridge  for  this  example  is  a  continuous  box  girder  bridge  of 
three  spans.     The  geometry  of  this  reinforced  concrete  bridge  is  shown  in 
figure  45. 

Because  the  column  heights  vary,  the  response  of  this  bridge  is  not 
expected  to  be  symmetrical.     Also,  because  the  ends  of  the  bridge  are 
torsi onally  restrained,  as  well   as  at  the  column  pairs,  no  significant 
longitudinal   torsional  deck  motion  is  expected.     Therefore,  adequate 
instrumentation  to  describe  all  significant  structure  distortions  is  given 
in  figure  46.     All   structure  instruments  are  located  at  mid-span  positions. 

If  U-(t)(i=l  ,*",7)  are  the  recorded  absolute  accelerations,  relative 
structure  accelerations  are  computed  as 

ui(t)  =  0^(t)-Uyg(t)  (168) 

U2(t)-Uyg(t)  (169) 

U3(t)-Uyg(t)  (170) 

U4(t)-U     (t)  (171) 


U2(t)  = 

U3(t)  = 

M^)  =  -4^-'    -xg' 

U^t)  -  U5(t)-U^  ft) 


(172) 
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Mode 

co(rps) 

T(sec) 

?(%) 

^x 

^y 

1 
2 
3 

4 

12.235 
19.051 
33.614 
33.614 

.514 
.330 
.187 
.187 

4.0 
4.0 
2.0 
2.0 

1.000 
1.000 

0 

0 

0 

0 
1.000 
1.000 

Table  11.-  Three-Span  Simply  Supported  Bridge, 
Model  Parameters  Identified 


Mode   w(rps)  T(sec)  E,{%)  y  y„ 

X       J' 


1 

12.236 

.513 

4.0 

1.000 

0 

2 

19.052 

.330 

4.0 

1.000 

0 

3 

33.615 

.187 

2.0 

0 

1.000 

4 

33.615 

.187 

2.0 

0 

1.000 

Table  12.-  Three-Span  Simply  Supported  Bridge, 
Exact  Model  Parameters 
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Figure  45.-  Continuous  three-span  bridge. 
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Figure  46.-  Continuous  three-span  bridge,  instrumentation  plan. 
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Ug(t)  =  Ug(t)-U^g(t)  (173) 

U7(t)  =  \}y{t)-\j^^{t)  (174) 

As  described  previously,  relative  velocities  u-  and  relative  displacements 
u-  are  numerically  calculated  from  the  relative  accelerations  computed  above, 
Again,  if  the  data  were  real,  some  filtering  would  undoubtedly  be  required 
to  remove  amplified  low-frequency  errors  in  the  computed  velocities  and 
displacements. 

For  this  example,  a  lumped  mass  approach  is  taken  and  coordinates 
describing  structure  distortions  (relative  motions)  corresponding  to 
accelerometer  directions  and  locations  are  assumed.  All  axial  distortions 
of  the  columns  and  bridge  deck  are  ignored. 

Again  the  mass  matrix  is  computed  from  bridge  geometry  and  material 
properties,  Masses  associated  with  the  three  lateral  degrees  of  freedom 
are  assumed  to  equal  one-third  of  the  total  bridge  mass;  the  mass 
associated  with  the  longitudinal  degree  of  freedom  is  assumed  to  equal  the 
total  bridge  mass;  and  masses  associated  with  the  three  vertical  degrees 
of  freedom  are  assumed  to  equal  half  of  each  span  mass,  Masses  are  chosen 
so  that  the  kinetic  energy  in  the  lumped-mass  system  approximately  equals 
the  kinetic  energy  in  the  actual  bridge,  when  vibrating.  The  mass  matrix 
and  vectors  of  participation  factors  become 
[M]  = 


39,54 

0 

0 

0 

0 

0 

0 

0 

39,54 

0 

0 

0 

0 

0 

0 

0 

39,54 

0 

0 

0 

0 

0 

0 

0 

118.62 

0 

0 

0 

0 

0 

0 

0 

19,77 

0 

0 

0 

0 

0 

0 

0 

19.77 

0 

0 

0 

0 

0 

0 

0 

19.77 
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«x>  = 


/ 


0 
0 
0 
118.62 
0 
0 
0 


«y>= 


/ 


39.54 
39.54 
39.54 

0 
0 


{M^}  = 


0 
0 
0 

0      ) 
19.77 

19.77 

19.77 


7 


Lumped-mass  models  are  always  inertially  uncoupled,  and  so  all  off-diagonal 
mass  terms  by  definition  equal  zero. 

For  the  time  histories  defined  by 

B^Ct)  =  -39.54  [ui(t)+U^g(t)] 

B2(t)  =  -39.54   [U2(t)+Uyg(t)] 

B3(t)  =  -39.54  [U3(t)+Uyg(t)] 

B^Ct)  =  -118.62[u4(t)+U^g(t)] 

B5(t)  =  -19.77[u5(t)+U^.(t)] 

Bg(t)  =  -19.77[ug(t)+U^g(t)] 

B/t)   =  -19.77[u7(t)+U     (t)] 
the  system  of  equations  describing  the  bridge  motions  becomes 


C„u,f  - 


S7"l^ 


+C,,U7+K,,u,+  •••  +K,,U,  =  B, 


+C„U7«^^u,+  -  +K„U7  =  B; 


(175) 


Numerical  errors  in  the  solution  for  the  unknown  parameters  C^-  and  ^,. 
again  are  greatly  reduced  if  the  motions  known  to  be  uncoupled  are  analyzed 
as  such.  It  is  assumed,  from  the  geometry  of  the  bridge  and  anticipated 
structural  behavior,  that  the  lateral  or  transverse  response  will  be 
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independent  of  the  longitudinal  and  vertical  motion.  Therefore,  the 
equations  of  motion  describing  the  bridge  motions  can  be  separated  into  the 
following  two  systems: 


Cll  V^12V^13^3^^1l"l^^l2V^13'^3  =  ^1 
^12^1^42'^2"'^23^3^^12"l"'^22V^23"3  =  ^2 
Cl3"l^^23^2-'S3^3^^13"l^^23V^33"3  =  ^3 


(176) 


and 


(177) 


C44V^45^5^^46^*6'-^47^7-"^44V^45"5^^46V^47"7  =  ^4 

C45V^55^^^56%-'S7^-*'S5V^55"5'-^56"6^'^57"7  =  ^5 

S6V^56^5''^66V^67^7"'^46VS6"5'"^66"6'"^67"7  =  ^6 
C47V^57VS7V^77"7"'S7V^57S^^67"6'"'^77"7  =  ^7 
Notice  that  the  damping  matrix  and  stiffness  matrix  are  assumed  to  be 
symmetric;  and,  by  the  assumption  that  these  two  systems  are  decoupled,  all 
other  C..  and  K- .  not  shown  are  equal  to  zero. 

•J  *  J 

Subroutine  FIT  was  used  to  find  the  best-fit  values  of  C-.  and  K- .  for 
the  two  systems  above.     Regrouped  back  into  a  single  system  of  equations 
[M]{u}+[C]{u}+[K]{u}  =  -(Mx^%"^^^^g"^^z^^'zg  ^^^^^ 

where  [M],  {H\,   {K,}>  and  {M  }  are  given  previously,  the  best  fit  damping 
X  y  u 

and  stiffness  matrices  found  are 
[C]  = 


64.1 

-24.5 

6.82 

0 

0 

0 

0 

-24.5 

61.4 

-4.74 

0 

0 

0 

0 

6.82 

-4.74 

53.7 

0 

0 

0 

0 

0 

0 

0 

54.1 

.288 

6.16 

-7.35 

0 

0 

0 

.288 

14.3 

-6.18 

-5.41 

0 

0 

0 

6.16 

-6.18 

.193 

-.800 

0 

0 

0 

-7.35 

-5.41 

-.800 

1.82 
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[K]  = 


0 
0 
0 
0 


,200 

0 

0 

0 

0 

,000 

0 

0 

0 

0 

,900 

0 

0 

0 

0 

0 

24,800 

-1330 

-2770 

2500 

0 

-1330 

10,800 

4510 

291 

0 

-2770 

4510 

17,300 

4140 

0 

2500 

291 

4140 

11,100 

At  this  point  the  mathematical  model   describing  the  significant  bridge 
distortions  is  completely  defined.     Again  it  is  not  in  its  most  useful 
form.     Mathematical  models  are  most  usefully  described  in  terms  of  normal 
coordinates  characterized  by  natural   frequencies,  natural  mode  shapes, 
damping  ratios,  and  modal  participation  factors. 


By  solving  the  related  eigenvalue  problem 
[[K]-a)^[H]]{({)}  =  {0} 
the  following  natural   frequencies  are  found 
oj^       =  17.7144     rps 

CO. 


(179) 


=  36.0280     rps 

=  64.1405     rps 

=  13.7408     rps 

=  19.5835     rps 

=  23.4614     rps 

=  32.6969     rps 

The  first  three  are  lateral  modes  of  vibration  and  the  following  four  are 
coupled  longitudinal   and  vertical  modes  of  vibration. 

Corresponding  to  each  natural   frequency  is  a  natural  mode  of  vibration. 
These  natural  mode  shapes,  in  physical   coordinates,  are  arranged  as  columns 


'2 

^5 

CO -7 
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in  the  transformation  matrix  M  as  follows: 


M   = 


.7280  -.7429  -.9934 

1.000  .2648  1.000 

.2760  1.0000  -1.0029 

0  0  0 

0  0  0 

0  0  0 

0  0  0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1.0000       1.0000         .0912       -.0227 

-.0344       8.6506       -.7880         .4511 

.3845     -6.2609       -.0520       1.0000 


-.5535       5.9418       1.0000 


.4200 


The  mass  matrix,  damping  matrix,  and  stiffness  matrix  are  transformed 
into  corresponding  mass,  damping  and  stiffness  matrices  in  normal 
coordinates  by 

[M]  =  [(Df  MU];   [C]  =  [<Df  [C][(j)];   [K]  =   [c^f  [K][(|,] 
By  definition,  all  off -diagonal   terms  of  the  mass  and  stiffness  matrices 
are  zero,  but  they  are  not  necessarily  zero  for  the  damping  matrix.     The 
transformed  matrices  become 
[M]  = 


63.508 

0 

0 

0 

0 

0 

0 

0 

64.135 

0 

0 

0 

0 

0 

0 

0 

118.329 

0 

0 

0 

0 

0 

0 

0 

127.623 

0 

0 

0 

0 

0 

0 

0 

3071.007 

0 

0 

0 

0 

0 

0 

0 

33.086 

0 

0 

0 

0 

0 

0 

0 

27.342 
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[C]  = 


63.915 

8.371 

3.278 

0 

0 

0 

0 

8.371 

90.378 

29.180 

0 

0 

0 

0 

3.278 

29.180 

250.442 

0 

0 

0 

0 

0 

0 

0 

67.855 

-29.404 

-3.644 

2.141 

0 

0 

0 

-29.404 

1209.331 

140.761 

-9.892 

0 

0 

0 

-3.644 

-140.761 

17.812 

-.387 

0 

0 

0 

2.141 

-9.892 

-.387 

-4.991 

m  = 


19,928. 

0 

0 

0 

0 

0 

0 

0 

83,248. 

0 

0 

0 

0 

0 

0 

0 

486,810. 

0 

0 

0 

0 

0 

0 

0 

24 

,094. 

0 

0 

0 

0 

0 

0 

0 

1,177, 

767. 

0 

0 

0 

0 

0 

0 

0 

18,213. 

0 

0 

0 

0 

0 

0 

0 

29,231 

The  vectors  of  participation  factors  become 


/ 


\ 


^v  -{ 


0 

0 

0 
118.62     ) 
118.62 

10.818 
,  -2.693 


^  79. 238 "^ 


{M,}  =  < 


20.636 
-39.394 
0 
0 
0 
0 


{M^} 


0 

0 

0 

=    {    -4.021 

164.714 

3.163 

^  36.992 


> 


Off-diagonal  damping  terms  are  in  general   not  negligible  but  are  small   in 
most  cases.     However,  those  off-diagonal  damping  terms  are  ignored,  and  it 
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is  assumed  that  the  equations  of  motion  in  normal  coordinates  can  be 
separated  without  significant  error.  Errors  will  be  evaluated  after  the 
analysis. 

Ignoring  the  off-diagonal  damping  terms  the  modal  equations  of 
motion  become 

nyi.006n^+313.8n^  =  -1.248U^  (180) 

fi2+1.409n2+1298.Ti2  =  --322*0^3  (181) 

n3+2.117ri3+4114.n3  =   -3330^^  (182) 

fi4+.532i^4+188.8n4  =  -.929U^g+.032U2g  (183) 

fi5+.394n5+383.5n5  =  -'0390     -.0540^  (184) 

fi5+.538fig+550.5n5  =  -.327U^  -.096U^  (185) 

n7-.183n7+1069.n7  =   .099U     -1.3530^  (186) 

Comparing  these  equations  to 

fi+2?wf|+a)  n  =  -Y  U     -Y  U     -Y  U  (187) 

i^x  xg  ^y  yg   'z  zg 

yields  the  modal   parameters  found  by  this  analysis  shown  in  table  13. 

The  mode  shapes  identified  are  shown  in  figure  47.     For  comparison,  the 

actual  modal   parameters  used  to  generate  the  artificial   time  histories  are 

shown  in  table  14,  and  the  exact  mode  shapes  are  shown  in  figure  48. 

From  this  example  analysis  it  is  apparent  that  frequencies, 

participation  factors,  and  mode  shapes  can  be  found  with  acceptable 

accuracy  using  this  technique.     Errors  in  the  damping  ratios  are 

significant  for  those  modes  of  vibration  not  significantly  excited  (note 

in  particular  the  negative  damping  ratio  found  for  the  seventh  mode  of 

vibration).     These  errors  in  damping  ratios  are  not  likely,  however,  to 
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Mode 

w(rps) 

e(%) 

Yx 

'y 

^z 

1 

17.714 

2.84 

0 

1.248 

0 

2 

36.028 

1.96 

0 

.322 

0 

3 

64.141 

1.65 

0 

-.333 

0 

4 

13.741 

1.94 

.929 

0 

-.032 

5 

19.584 

1.01 

.039 

0 

.054 

6 

23.461 

1.15 

.327 

0 

.096 

7 

32.697 

-.28 

-.099 

0 

1.353 

Table  13  -  Three-Span  Continuous  Bridge  Model 
Parameters  Identified 


Mode 

aj(rps) 

a%) 

^x 

^y 

^z 

1 

17.780 

2.00 

0 

1.1888 

0 

2 

36.040 

2.00 

0 

.3291 

0 

3 

64.310 

2.00 

0 

-.4002 

0 

4 

13.730 

2.00 

.9380 

0 

-.0446 

5 

19.660 

2.00 

.0413 

0 

.0639 

6 

23.360 

2.00 

.3240 

0 

.0546 

7 

31.740 

2.00 

-.6340 

0 

1.6910 

Table  14.-  Three-Span  Continuous  Bridge  Exact 
Model   Parameters 
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Figure  47.-  Continuous  three-span  bridge,  mode  shapes  found  from  analysis. 
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Figure  47  (continued).-  Continuous  three-span  bridge,  mode  shapes  found  from 

analysis. 
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Figure  48.-  Continuous  three-span  bridge,  actual  mode  shapes. 


139 


z 

A 


i>x   — 


Mode  4 


Z 
A 


^>X 


Mode  5 


Z 
A 


-[>x 


z 

A 


H>x 


Mode  6 


Mode  7 


Figure  48  (continued).-  Continuous  three-span  bridge,  actual  mode  shapes, 
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produce  significant  errors  in  response.  Although  response  amplitudes  of 
resonating  systems  are  highly  dependent  upon  damping  ratios,  response 
amplitudes  of  seismically  induced  motions  generally  are  not  highly 
sensitive  to  variations  in  damping  within  the  range  of  values  to  be  found 
in  practice.  One  general  rule  of  thumb  in  time  history  analyses  is  that 
if  the  response  is  particularly  sensitive  to  a  parameter,  that  parameter 
will  be  found  with  accuracy. 
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V.   CONCLUDING  REMARKS 

This  report  has  discussed  methods  of  analyzing  stronq-motion  records 
obtained  from  highway  bridges.  There  have  been  two  objectives  of  the 
analyses:  (1)  obtain  maximum  force  levels  throughout  the  bridge,  and  (2) 
identify  a  mathematical  model  which  best  describes  the  observed  dynamic 
behavior.  There  are  few  uncertainties  regarding  the  calculation  of  maximum 
force  levels,  and  few  assumptions  needed  to  be  made.  On  the  other  hand, 
many  uncertainties  are  associated  with  the  identification  of  the  best-fit 
mathematical  model.  Some  of  those  uncertainties  and  additional  limitations 
are  discussed  here. 

This  report  has  one  procedure  to  identify  one  mathematical  model  of 
observed  dynamic  behavior.  Two  fundamental  assumptions  have  been  made  about 
that  model:   (1)  the  model  is  linear,  and  (2)  all  support  motions  are 
identical  and  in  phase.  These  assumptions  may  seem  overly  restrictive  since 
it  is  well  known  that  bridge  response  can  be  highly  non-linear,  and  support 
motions  can  vary  significantly,  particularly  if  the  foundation  material 
differs  from  one  support  to  another.  Unfortunately,  the  ability  to  identify 
non-linear  models  with  variable  support  motions  lags  far  behind  the  ability 
to  formulate  such  models  and  generate  acceptable  response.  In  fact,  it  is 
almost  beyond  the  state  of  the  art  to  identify  linear  mathematical  models 
with  identical  support  motions  when  the  natural  frequencies  of  the  bridge  are 
closely  spaced.  Unfortunately,  this  is  often  the  case  for  highway  bridges. 
The  procedure  described  in  this  report  will  yield  adequate  results  for 
structures  that  do  or  do  not  have  closely  spaced  frequencies.   In  fact,  it 
was  developed  with  that  purpose  in  mind. 

Any  non-linear  model  can  be  interpreted  as  a  linear  model  with  time- 
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dependent  parameters.  Since  the  identification  procedure  described  in  this 
report  can  identify  a  best-fit  linear  model  over  a  small  time  slice, 
variations  in  the  best-fit  linear  models  identified  can  be  analyzed  as 
non-linearities  in  the  bridge  behavior.  Bridge  motions  with  variable 
support  motions  can  also  be  analyzed  with  this  identification  procedure. 
In  many  cases,  embankment  motions  could  be  interpreted  not  as  input  motions 
to  the  bridge  but  as  response  motions  to  bedrock  motion.  If  it  is 
appropriate  to  assume  that  column  base  motions  are  identical  to  bedrock 
motions  under  an  embankment,  then  embankment  motions  can  be  included 
with  other  bridge  motions  as  response  motions,  and  can  be  analyzed  using 
the  procedure  described  in  this  manual.  Bridge  motions  produced  from 
ground  motions  propogating  from  one  support  to  another,  however,  cannot 
be  analyzed  using  this  procedure. 
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APPENDIX  A 
GLOSSARY  OF  TERMS 
{B(t)}  vector  of  known  inertial  time  histories  in  physical 

coordinates 
{b(t)}  vector  of  known  inertial  time  histories  in  generalized 

coordinates 
C  constant 

[C]  damping  matrix  in  physical  coordinates 

[C]  damping  matrix  in  normal  coordinates 

[c]  damping  matrix  in  generalized  coordinates 

c  half  of  deck  width 

E  squared  error 

~  ■f'h 

F^(t)  force  on  i   lumped  mass 

~  th 

F-(s,t)  Force  per  unit  length  for  the  j       generalized  (normal) 

J 

shape 
F  •  damping  force  in  physical  coordinates 

F|^-  elastic  force  in  physical  coordinates 

F^  inertial  force  in  physical  coordinates 

f  •  damping  force  in  generalized  coordinates 

CI 

f.  .  elastic  force  in  generalized  coordinates 

f^.  inertial  force  in  generalized  coordinates 

mi  -^ 

I(x)  rotational  inertial  per  unit  length 

i ,  j  indeces 

T,j,k  unit  vectors  in  the  Cartesian  coordinate  directions 

x,y,z  respectively 
[K]  stiffness  matrix  in  physical  coordinates 
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[K]  stiffness  matrix  in  normal  coordinates 

[k]  stiffness  matrix  in  generalized  coordinates 

L  length  of  bridge 

M(x)  mass  per  unit  length 

M  ,M  ,M  internal  bending  moments  in  bridge  deck 

A        J         Z. 

[M]  mass  matrix  in  physical  coordinates 

{My},{M  },{M  }        vectors  of  participation  factors  in  physical 

coordinates 

[M]  mass  matrix  in  normal   coordinates 

{M  },{M^}j{M  }  *       vectors  of  participation  factors  in  normal 

A         y         ^ 

coordinates 
[m]  mass  matrix  in  generalized  coordinates 

vectors  of  participation  factors  in  generalized 

coordinates 
internal  axial   force  in  bridge  deck 
number  of  lumped  masses 
generalized  force 
i       reaction 
i       reaction  to  a  unit  j       generalized  (normal) 

acceleration 
i       reaction  to  unit  rigid-body  accelerations  in 

the  X-,  y-,  and  z-  directions  respectively 
total   reaction  in  the  x-,  y-,  and  z-  directions  due 

to  unit  rigid-body  accelerations  in  the  x-,  y-, 

and  z-  directions  respectively 


{m^  >{m}  ,{m^} 

\ 

n 

Qi 

Ri(t) 

^ij 

^ix'^iy'^-z 

^tx'V'^^tz 
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'^tix»'^tiv''^tjz        total  reaction  in  the  x-,  y-,  and  z-  directions 

respectively  due  to  a  unit  jth  generalized 

(normal)  acceleration 
rj(s)  vector  of  distortions  due  to  a  unit  displacement 

corresponding  to  the  ith  reaction 
Sg((jo)  spectral  acceleration 

s  coordinate  along  the  principal  axis  of  the  bridge 

T  natural  period  (=27t/w) 

t  temporal  coordinate 

{Ll(t)}  vector  of  absolute  bridge  displacements 

Uj^q(t),Uyq(t),U  (t)    ground  accelerations  in  x-,  y-,  and  z-  directions 

respectively 
{u(t)}  vector  of  bridge  distortions  relative  to  the  ground 

u(s,t)  vector  distortion  relative  to  the  ground  of  point  on 

bridge  axis 

Vv,V„  internal  shears  on  bridge  deck 

^    y 

x,y,z  Cartesian  coordinates 

{a(t)}  vector  of  bridge  distortions  in  generalized 

coordinates 
{ Y^},{  Yy},{  y^}        vectors  of  participation  factors  in  normal  coordinates 

when  masses  are  normalized  to  unity 
6  virtual  displacement 

{n(t)}  vector  of  bridge  distortions  in  normal  coordinates 

[e]  transformation  matrix  from  normal  to  generalized 

coordinates 
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iO  vector  of  modal  viscous  damping  ratios 

[(J)]  trans fonnati on  matrix  from  normal  to  physical 

coordinates 

$^(s)  ith  normal  mode  shape 

M  transformation  matrix  from  generalized  to  physical 

coordinates 

J-Cs)  ith  generalized  shape 

{o)}  vector  of  natural  circular  frequencies  of  vibration 

[  ]  brackets  denoting  a  matrix  quantity 

[  ]  transpose  of  a  matrix 

{  }  braces  denoting  a  vector  quantity 

Z  summation 

(•)=d/dt  differentiation  with  respect  to  time 

(~)  symbol  of  a  vector  quantity 
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APPFNDIX  B 

iNTRoniicTinN  in  stroctiiral  oynamics  of  bridges 

Before  any  mathematical  model  of  dynamic  bridge  behavior  can  be  made,  a 
coordinate  system  must  be  established  which  can  describe  the  bridge  motions 
accurately.  The  number  and  location  of  the  coordinates  to  describe  the 
bridge  motion  is  purely  arbitrary.  To  describe  the  bridge  motion  exactly 
requires  an  infinite  number  of  coordinates.   If  the  coordinates  are  chosen 
judiciously,  however,  only  a  few  coordinates  are  necessary  to  describe  the 
complete  bridge  motion  with  acceptable  accuracy.  In  the  case  of  a  curved 
highway  grade  separation  bridge,  for  example,  there  are  many  possible 
coordinate  systems  that  may  be  chosen  to  describe  the  horizontal  motion  of 
the  bridge.  If  vertical  and  twisting  motion  of  the  bridge  deck  are  ignored 
(such  motions  may  be  wery   important  but  are  not  considered  here  in  order  to 
simplify  the  discussion),  then  the  bridge  deck  motions  and  ground  motions 
can  be  described  by  the  coordinates  shown  in  figure  49,  for  a  total  of 
twenty  coordinates.  The  number  of  required  coordinates  can  further  be 
reduced  if  it  is  possible  to  make  some  additional  assumptions.  For  most 
bridges,  or  isolated  segments  of  bridges,  it  is  reasonable  to  assume  that 
bridges  are  inextensible  in  their  longitudinal  directions  (the  s  direction 
for  the  bridge  shown). 

In  many  cases,  it  is  also  reasonable  to  assume  that  the  ground  motion 
exciting  the  bridge  is  identical  under  each  support.  This  assumption  is 
reasonable  for  short  span  bridges  on  firm  foundation  material.  For  long 
span  bridges  (such  as  the  San  Francisco  Golden  Gate  Bridge),  it  is  of  course 
unreasonable  to  assume  that  the  earthquake  motion  at  one  end  is  identical  to 
the  earthquake  motion  at  the  other  end.   It  may  also  be  erroneous  to  make 
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Plan  View 


Deck  Instrumentation 


Ground  Instrumentation 


Figure  49.-  Example  curved  bridge  with  complete  instrumentation, 
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this  assumption  for  a  simple  highway  bridge  if  a  central  pier  rests  on  firm 
rock,  while  the  end  abutments  rest  on  fill.  In  any  case  the  mathematical 
modeling  assumptions  are  made  primarily  for  mathematical  convenience,  and 
each  should  be  evaluated  from  the  analysis  of  the  strong-motion  records. 

When  the  two  assumptions  described  above  are  made,  the  bridge  motion 
may  be  described  accurately  by  the  reduced  coordinates  shown  in  figure  50. 
In  this  case,  the  bridge  and  ground  motions  are  described  by  a  total  of  only 
eight  coordinates.  The  five  lateral  motion  coordinates  at  the  bridge  deck 
describe  accurately  the  motions  of  those  five  points  only.  The  lateral 
motions  of  the  deck  must  be  assumed  for  all  locations  between  coordinate 
points.  For  typical  highway  bridges,  significant  motions  involving  more 
than  a  half  wave  length  per  span  are  not  expected. 

Instead  of  approximating  the  bridge  displacements  with  the 
displacements  of  a  few  selected  points,  the  bridge  displacements  can  be 
approximated  by  a  linear  combination  of  assumed  independent  shape  functions. 
For  the  assumed  shapes  shown  in  figure  51  and  corresponding  coordinates 
a-iCt),  apCt),  and  a^lt),  the  total  lateral  bridge  deck  displacement  relative 
to  the  ground  can  be  given  by 

u(s,t)  =  $^(s)a^(t)-h52(s)a2(t)-hf3(s)a3(t)  (B-1) 

The  displacements  in  these  coordinates  a-(t)  are  not  in  themselves 
measurable  quantities  but  are  related  to  measurable  quantities  u^-(t)  by  the 
assumed  deflected  shapes  ijj-(s).  The  bridge  is  not  expected  to  vibrate 
freely,  to  any  significant  extent,  in  a  mode  like  JoCs),  but  this  assumed 
shape  should  be  included  nonetheless  to  improve  the  description  of  the 
actual  lateral  bridge  deck  motion,  otherwise  described  entirely  by  the 
assumed  shapes  ^-,{s)   and  ifoCs). 
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Plan  View 


Deck  Instrumentation 


Ground  Instrumentation 


Figure  50.-  Example  curved  bridge  with  reduced  instrumentation. 
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ii(s) 


h(^) 


^lis) 


Figure  51.-  Example  curved  bridge,  assumed  deflected  shapes. 
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In  addition  to  the  three  coordinates  that  describe  the  lateral  bridge 
motion,  a  single  longitudinal  coordinate  and  two  ground  motion  coordinates 
are  also  required.  It  is  unlikely  that  fewer  coordinates  could  describe  the 
horizontal  bridge  motion  with  equal  accuracy.  It  is  also  very   important 
that  significant  bridge  motions,  yielding  significant  force  levels  are  not 
lost  when  coordinates  are  hastily  removed  for  the  sake  of  simplicity  and 
reduced  cost.  Once  the  information  is  lost,  it  can  never  be  regained  from 
the  remaining  coordinates. 

The  mathematical  model  used  to  describe  the  dynamic  behavior  of  a 
bridge  consists  of  a  set  of  equations  of  motion.  There  will  exist  one 
equation  for  each  unknown  bridge  displacement.  For  the  coordinate  system 
described  above  there  will  exist  four  equations  of  motion.  Each  equation  of 
motion  is  nothing  more  than  an  equation  of  dynamic  equilibrium  including  all 
forces  acting  on  the  bridge  in  that  particular  coordinate  "direction". 
Forces  acting  in  a  coordinate  direction  that  is  not  directly  measurable  will 
be  defined  later. 

Specifically,  there  will  exist  inertial  forces,  damping  forces  (forces 
that  damp  the  motion  by  dissipating  energy),  elastic  forces,  and  externally 
applied  forces.  For  the  ith  coordinate  direction,  the  equation  of  dynamic 
equilibrium  can  be  written 


f  .  +  f  .  +  f ,  .  =  f  .  (B-2) 

mi    ci    ki    pi  ^   ' 


where 


fm 

= 

inertial   force. 

fc1 

= 

damping  force. 

fki 

= 

restoring  force,  and 

fpi 

= 

applied  force. 
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A  generally  accurate  assumption  for  non-damaginn  motion  is  that  the 
above  inertial,  damping,  and  elastic  forces  are  proportional  to  the 
coordinate  accelerations,  velocities,  and  displacements  respectively.  This 
assumption  of  linear  behavior  is  one  of  the  fundamental  assumptions  to  be 
checked  by  the  experimental  data.  The  assumption  of  linear  behavior  is  not 
made  primarily  because  linear  models  are  best  but  because  linear  models  can 
be  solved  in  closed  form  using  general  solution  methods.  The  same  cannot  be 
said  for  non-linear  equations.  Whether  or  not  a  linear  equation  is  the  best 
to  simulate  bridge  motion  is  to  be  evaluated  by  the  experimental  data. 

Like  u(s,t)  the  displacements  ai(t)  are  defined  as  displacements 
relative  to  the  ground.  Relative  coordinates  are  important  because  force 
levels  in  bridges  are  proportional  to  relative  displacements;  damping  forces 
are  proportional  to  relative  velocities.  Inertial  forces,  on  the  other 
hand,  are  proportional  to  absolute  accelerations. 

In  terms  of  constants  of  proportionality  the  inertial  force  in  the  ith 
coordinate  direction  can  be  expanded  in  the  most  general,  linear  manner  to 

^mi  =  mu'^]+   .-.  "iinC^n+%iUxg+'^yiUyg+f^zi'^zg  (B-3) 

where  U^q,  Uyq,  and  U^q  are  the  three  mutually  perpendicular  ground  accel- 
erations (the  definitions  and  significance  of  m-|-|,  ...,  m^p,  mxi»  myi,  and 
m2i  are  discussed  below).  Similarly,  by  definition,  the  damping  force  in 
the  ith  direction  is  given  by 

fdi  =  c-j]a-i+c-j2a2  •••  "^c-jnap  (B-4) 

and  the  elastic  force 

^ki  =  kTia-|+ki2a2+  •••  ^^in^n  ^^"^^ 

where 

aj  =  dttj/dt  (B-6) 
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a-  =  d  a./dt 

J  J 


(B-7) 


Earthquakes  do  not  apply  forces  directly  to  bridges,  so  for  all  i, 
f  .=0,  All  of  the  forces  induced  by  earthquakes  initially  are  inertial. 
Then,  only  after  relative  motion  between  the  bridge  and  ground  occurs,  are 
damping  and  elastic  forces  developed. 

The  ith  of  n  equations  of  motion  can  then  be  written  as 
m^.^a^+  ...  +m.^a^+c.^a^+  ...  +c.^a^+k.^a^+  ...  +k.^a^  ^ 


-"^xi^xg-^i^g-"^ziUzg 


(B-8) 


For  bookkeeping  reasons,  these  n  equations  of  motion  can  be  regrouped  in 
matrix  form  (Gere  &  Weaver,  1965). 
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In  simpler  form  this  can  be  written  as 

[m]{a}+[c]{a}+[k]{a}  =  -Cnix>Uxg^%>Uyg<"iz^°^zg 


(B-9) 


(B-10) 
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The  physical  significance  of  the  terms  m-.,  c-  .,  and  k-.  can  best  be 
obtained  from  equations  B-3,  B-4,  and  B-5.   Specifically  the  term  m- •  is 
the  force  in  the  i  coordinate  direction  due  to  a  unit  acceleration  in  the  j 
direction  only.  The  term  c • ■  is  the  force  in  the  i  coordinate  direction  due 
to  a  unit  velocity  in  the  j  direction  only.  The  term  k- .  is  the  force  in 
the  i  coordinate  direction  due  to  a  unit  displacement  in  the  j  direction 
only.  Similarly  the  terms  m^ • ,  m  • ,  and  m_.  are  forces  in  the  i  coordinate 
direction  due  to  unit  ground  accelerations  U  ,  U  ,  and  U   respectively, 
with  all  relative  bridge  motions  fixed. 

The  numerical  evaluation  of  these  terms  is  possible  from  the 
definitions  above  and  with  use  of  the  theorem  of  virtual  work  (Clough  and 
Penzien,  1975).  Only  the  inertial  terms  will  be  evaluated  using  this 
theorem  as  they  are  the  only  terms  that  need  to  be  evaluated  from  the 
weights  and  dimensions  of  the  bridge  for  the  analysis  procedures  that  have 
been  described  in  this  report.  The  damping  and  elastic  terms,  in  one  form 
or  another,  are  those  modeling  parameters  that  are  to  be  identified  from  the 
bridge  records.  Formal  procedures  for  their  theoretical  formulation  are 
available  in  any  of  many  texts  on  structural  dynamics  (Biggs,  1964,  Clough 
and  Penzien,  1975,  Hurty  and  Rubinstein,  1964). 

The  theorem  of  virtual  work  is  stated  as  follows  (Clough  and  Penzien, 
1975):  "If  a  structure  which  is  in  static  or  dynamic  equilibrium  under  the 
action  of  a  set  of  forces  is  subjected  to  a  virtual  displacement,  the  total 
work  done  by  those  forces  equals  zero." 

For  the  determination  of  m- .  consider  the  bridge  subjected  to  a  unit 
acceleration  in  the  j  coordinate  direction  only.  All  velocities  and 
displacements  are  assumed  to  equal  zero.  Considering  now  the  bridge  to  be 
composed  of  many  line  elements  (including  the  deck  as  well  as  piers),  the 
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acceleration  at  any  point  on  the  bridge  is  given  by 

d^u(s,t)/dt^  =  ii)j(s)aj(t)  =  ^.{s)  (B-11) 

where  the  last  equality  is  due  to  the  assumption  that  a.(t)  =  1.  The 
•'acceleration  u  and  assumed  deflected  shape  ^jj   are  in  general  vector 
quantities  having  components  in  directions  other  than  normal  to  the 
longitudinal  axis  s.  For  any  element  of  the  bridge.  As,  the  inertial  force 
due  to  that  unit  acceleration  is  given  by 

AF(s)  =  M(s);i)j(s)As  (B-12) 

where  M(s)  is  the  mass  per  unit  length  of  bridge  element. 

Let  (k  be  a  force  corresponding  to  the  k  coordinate.  For  a  set  of 

virtual  displacements  6uk  =  ij;.6a,(k=l,  2 n)  the  total  work  produced  by 

the  \^\   plus  the  total  work  produced  by  the  inertial  forces  undergoing  the 
virtual  displacements  must  add  to  zero  by  virtue  of  the  theorem  of  virtual 
work: 

Q^6a^+  ...  +Q.6a.+  ...  +Q^6a^+/M(s)ii;j(s) '60^  (s,t)  ds 

+...+/M(s)ii;j(s)'6Uj(s,t)ds+...+/M(s)ii;j(s)-6u^(s,t)ds  =  0    (B-13) 
where  the  integrations  are  over  the  entire  bridge.  There  are  likely  to 
exist  many  other  forces,  such  as  reactions,  that  result  from  the  unit 
acceleration.  However,  since  motions  of  the  bridge  are  described  completely 
by  the  coordinates,  a.  ,  at  least  to  within  an  acceptable  error,  none  of 
those  reactions  or  other  forces  will  produce  work. 

If  all  but  the  ith  virtual  displacements  equal  zero,  then 

-Q^-6a.  =  /M(s)6u.(s,t)>j(s)ds  (B-14) 

Since  6u.(s,t)  =  ^,-(s)6a.  and  since  6a-  is  non-zero,  6a.  can  be  divided  from 
both  sides  of  equation  B-14  to  yield 

-Q^.  =  /M(s)$.(s)-{j;.(s)ds  (B-15) 
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Instead  of  considering  the  force  in  the  i  coordinate  direction  resisting  the 
acceleration,  consider  the  equivalent  force  in  the  i  coordinate  direction 
produced  by  the  unit  acceleration.  By  definition  this  force  is  m^.,  and 
therefore 

m^j  =  /M(s)5;.(s)-;i;j(s)ds  (B-16) 

In  a  similar  manner  the  inertial  terms  m  •,  m  •,  and  m  •  can  be  found. 

XI   y 1       zi 

Assume  that  U  =1 ,  while  holding  U   and  U   and  all  relative  bridge 

accelerations  fixed.  For  T  equal  to  the  unit  vector  in  the  U   direction, 

xg 

an  elemental  inertial  force  due  to  this  unit  ground  acceleration  is 

AF(s)  =  M(s)Tas  (B-17) 

If  the  bridge  is  subjected  to  the  single  virtual  displacement  6u.(s,t),  then 

from  the  theorem  of  virtual  work, 

Q^6a^.+/M(s)6u^.(s,t)*Tds  =  0  (B-18) 

or  again, 

q.+/M(s)j;.(s)-ids  =  0  (B-19) 

By  definition,  m  .  is  the  force  in  the  i  coordinate  direction  produced  by  a 

A  I 

unit  acceleration  U     ;  therefore 

xg 

\^  =  /M(s);i;.(s)Mds  (B-20) 

In  an  analogous  manner, 

Tiyi  =  /M(s)ij;.(s)-jds  (B-21) 

and 

m^i  =  /M(s)f.(s)-kds  (B-22) 

where  J  and  k  are  unit  vectors  in  the  directions  of  0   and  U  , 
respectively. 

Equation  B-9  rewritten  is  again 
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[m]{a}+[c]{a}+[k]{a}  =  -i\}\-i^y}\^-i^/%  (B-23) 

where  the  inertial  terms  are  defined  by  equations  B-16,  B-20,  B-21 ,  and  B-22 
in  terms  of  the  assumed  deflected  shapes  liJ-(s),  In  practice  the  terms  k-. 
can  also  be  found  in  a  similar  manner.  The  terms  c-.  are  rarely  found  in 
this  form. 

In  general,  all  elements  of  the  mass  matrix  [m],  the  damping  matrix 
[c],  and  the  stiffness  matrix  [k]  are  non-zero,  making  the  solution 
difficult.  There  does  exist  a  special  set  of  deflected  shapes,  $^(s), 
called  normal  mode  shapes,  and  corresponding  coordinates  n^Ct)  that  greatly 
simplify  the  solution  of  equation  B-10.  Instead  of  expanding  the 
displacement 

u(s,t)  =  $^(s)a^(t)+  ...  +;i;^(s)a^(t)  (B-24) 

expand  it  in  terms  of  normal  modes: 

u(s,t)  =  $^(s)n^(t)+  ...  +$^(s)n^(t)  (B-25) 

Then  all  off  diagonal  masses 

frJ^j  =  /M(s)J.(s).J(s)ds  (B-26) 

equal  zero  for  ij^j,  as  well  as  all  off-diagonal  terms  of  [C]  and  [K],  It 
can  be  shown  that  for  all  real,  symmetric,  positive  definite  mass  and 
stiffness  matrices  for  stable  structural  systems  a  set  of  normal  modes  can 
in  fact  be  found  such  that  all  M-.  =  0  and  K.-.  =   0  for  i  ^  j  (Biggs,  1964; 
Clough  and  Penzien,  1975;  Hurty  and  Rubinstein,  1964),  In  general,  the 
off-diagonal  terms  C-.  may  not  be  zero  but  quite  often  are  very   small.  In 
such  cases  it  is  usual  to  assume  that  they  are  zero  so  that  all  matrices 
[M]»  [C],  and  [K]  are  di agonal i zed  by  the  coordinate  selection. 

Standard  methods  to  find  these  normal  coordinates  are  outlined  in  texts 
on  structural  dynamics  (Biggs,  1964;  Clough  and  Penzien,  1975;  Hurty  and 
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Rubinstein,  1964)  and  will  not  be  considered  further  here.     To  analyze  the 
vibration  records  it  is  enough  to  know  that  normal  modes  exist,  and  there  is 
no  special  need  to  know  how  to  find  them  analytically. 

For  the  specific  set  of  normal  coordinates,  equation  B-10  is  written 

[M]{n}+[C]{n}+[/(]{n}  =  -m/u^^-U^/Uy^-m^}^^^^  (B-27) 

where  [M],  [C],  and  [K]  are  diagonal  matrices  (having  all  off-diagonal  terms 
equal   to  zero).     Because  of  the  special   diagonal   form  of  all  matrices,  the  n 
equations  of  motion  can  be  written  as  n  uncoupled  si ngle-degree-of- freedom 
simple-oscillator  equations  of  motion,  the  ith  of  which  is 

M.n.+C.n.+K.Ti-   =  -M   .0     -M   .0     -M   -U  (B-28) 

I'l     1  'i     1  'i  XI  xg    yi  yg    zi  zg  \"  ^  j 

where  the  inertial  terms  are  defined  by 

M^  =  /M(s)J.(s)-J.(s)ds  (B-29) 

Mj^^  =  /M(s)J.(s)-Tds  (B-30) 

Myi  =  /M(s)$.(s)-jds  (B-31) 

and 

M^i  =  /M(s)$.(s)-kds  (B-32) 

Again,  all  integrals  are  over  all  elements  of  the  bridge. 

If  equation  B-28  is  divided  by  M-  then  this  modal  equation  of  motion 

can  be  written  in  the  alternate  form 

n.+2r.(o.r|.+(o?n.  =  -y  -U  "Y  -U  -y  U  (B-33) 

h  ^in  h  ^v'l    'xi  xg  'yi  yg  ''z  zg 

where  the  obvious  definitions  are  made 
^^.  =  C^/2/KTg7  (damping  ratio) 
w-  =  /fC . /M .   (natural   undamped  circular  frequency  of  vibration) 

'xi         xr    1 

Y  .  =  ^y^/W-j  (participation  factors) 
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Yz1  =  «zi/"i 

The  0)-  have  the  physical  significance  that  a)-=2Trf.  where  f-  is  the 

approximate  natural  frequency  of  vibration  (in  Hz)  of  the  ith  mode.  It  is  a 

property  of  the  mode  that  if  the  mode  could  be  perturbed  alone,  it  would 

vibrate  solely  at  that  frequency.  The  damping  ratio  C,-  ^Iso  has  the 

physical  significance  of  equaling  the  ratio  of  damping  in  the  bridge 

(capability  of  dissipating  energy)  to  the  value  of  critical  damping. 

Critical  damping  is  defined  as  the  minimum  damping  for  which  a  perturbed 

shape  will  return  to  its  initial  shape  exponentially  rather  than  in  an 

oscillatory  manner.  In  practice,  most  bridges  are  very  lightly  damped 

having  E,<0,^ . 

Finally  it  can  be  said  that  a  mathematical  model  of  a  bridge,  if 

formulated  under  the  assumptions  presented  here,  can  be  characterized  by  the 

natural  frequencies  w- ,  the  modal  damping  ratios  ^.,  and  the  participation 

factors  Yv-i»  Yui »  ^"^  "Y^T*  The  participation  factors  are  solely  functions 

of  the  mass  distribution  M(s)  and  mode  shapes  $-(s),  the  usual  third  modal 

parameter.  The  normal  mode  formulation  of  the  problem  is  powerful  because 

typical  bridge  responses  may  be  yery   accurately  described  in  terms  of  a  very 

few  normal  coordinates. 
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APPENDIX  C 
COMPUTER  PROGRAM  LISTINGS 
Subroutine  FIT 

Subroutine  FIT  solves  for  the  best-fit  values  of  A(I,J)  such  that 
the  left  side  of 

[A]{F(t)}  =  {B(t)} 
fits  the  right  side  in  a  least-squares  sense.     {F(t)}  and  {B(t)} 
are  known  time  histories. 
OUTPUT 

A(NN,NM)  Matrix  of  unknown  coefficients  to  be  found 

[A]  =  [[C][K]] 


INPUT 


B(NN»NP) 
F(NM,NP) 
IJ(NN»NM) 


NE 

NM 

NN 

NP 

NS 

NF 

P(NE,N£) 

Q(NE) 

R(NE) 


Vector  of  known  time  histories 

Vector  of  known  time  histories 

Array  of  unknown  coefficient  indices: 

IJ(K,L)  is  the  unknown  coefficient  index  for 

the  term  A(K,L) 

Number  of  independent  unknown  coefficients 

in  [A] 

Twice  NN 

Number  of  system  degrees  of  freedom 

Number  of  points  in  known  time  histories 

Starting  index  of  known  time  histories 

Finishing  index  of  known  time  histories 

Scratch  array 

Uo. 

Uo. 
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subroutine  fit(a,b,f,p,q,r,ij,ne,nm,nn,np,ns,nf,err) 

dimension  a(nn,nm)  ,b(nn,np) ,f(nm,np) ,ij(nn,nm) 

dimension  p(ne ,ne) ,q(ne) ,r ( ne) 

do  10  i  =  1  ,ne 

q(i)=0. 

do  10  j=1  ,ne 
10  p(i,j)=0. 

do  50  1  =  1  ,ne 

do  50  i  =  1  ,nn 

do  40  j=1 ,nm 

if  (ij(i,j)  .ne.l)  go  to  40 

do  20  k=1 ,nm 

m=ij( i  ,k) 

do  20  n=ns,nf 
20  p(l,m)=p(l,m)+f(k,n)«f( j,n) 

do  30  n=ns,nf 
30  q(l)=q(l)+b(i,n)»f(j,n) 
40  continue 
50  continue 

call  solve(p,q ,r ,ne) 

do  60  i  =  1  ,nn 

do  60  j=1  ,nm 

k=ij(i,j) 
60  a(i,j)=r(k) 

err=0 . 

cbb=0. 

do  80  i=1  ,nn 

do  80  n=ns,nf 

dum=b(i,n) 

cbb=cbb+dum*dum 

do  70  j=1  ,nm 
70  dum=dum-a(i , j)*f( j ,n) 

err=err+dum*dum 
80  continue  ' 

err=sqrt ( err/ebb) 

return 

end 


165 


Subroutine  SOLVE 

Subroutine  SOLVE  solves  the  system  of  n  linear  equations 
[A]{X}  =  {B}  for  the  vector  of  unknowns,  {X},  using  the  Gaussian 
elimination  procedure.  The  matrix  [A]  and  vector  {B}  are  destroyed 
in  the  process. 
OUTPUT 

Vector  of  unknowns 


X(N) 
INPUT 

A(N,N) 

B(N) 

N 


Matrix  of  coefficients 

Vector  of  knowns 

Number  of  equations  to  be  solved 
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subroutine  solve( a ,b ,x ,n) 

dimension  a(n  ,n)  ,b(n) ,x (n) 
c  subroutine  solve  solves  the  set  of  linear  equations  of  order  n 
c  [a][x]=[b] 

do  50  k=2,n 

m=k-1 

nn  =  0 
5  if(a(m,m))  40,10,40 
10  nn=nn+1 

if(nn.gt.n)  go  to  200 

do  20  j  =  1  ,n 

ad=a(m,j) 

a(m, j )  =  a(m+nn  ,  j  ) 
20  a(m+nn , j )=ad 

bd=b(m) 

b(m)=b(m+nn) 

b(m+nn)=bd 

go  to  5 
40  do  50  i  =  l<,n 

c=a( i  ,m) 

b(i)=b(i)-b(m)*c/a(m,m) 

do  50  j=m,n 

a(i,j)=a(i,j)-a(m,j)*c/a(m,m) 
50  continue 

x(n)=b(n)/a(n ,n) 

l  =  n-1 

do  70  k=1 ,1 

m=n-k 

dum=0. 

do  60  j=m,l 
60  dum=dum+a(m, j+1 )*x( j+1 ) 
70  x(m)= (b(m)-dum)/a(m,m) 

go  to  30C- 
200  write(6,210) 

210  format(//,"  singular  matrix") 
300  continue 

return 

end 
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APPENDIX  D 
CONVERSION  FACTORS  TO  SI  METRIC  UNITS 


inches  (in) 
inches  (in) 
inches  (in) 

feet  (ft) 
yards  (yd) 
miles  (mi) 

degrees  (°) 


meters  (m) 
centimeters  (cm) 
millimeters  (mm) 

meters  (m) 
meters  (m) 
kilometers  (km) 

radians  (rad) 

hectares  (ha)  ^ 
cubic  meters  (m^) 
cubic  meters  (m  ) 
liters  (1) 

kilograms  (kg) 
kilograms  (kg) 

newtons  (N) 
newtons  per  sq  m 

(N/m2) 
newtons  per  sq  m 

(N/m2) 

joules  (J) 
watts  (W) 
joules  (J) 


Conversion  Factor 

0.0254 
2.54 
25.4 

0.305 
0.914 
1.609 

0.0174 


0.405 
1233 
3.79  X  10 
3.79 


-3 


acres  (acre) 
acre-feet  (acre- ft) 
gallons  (gal) 
gallons  (gal) 

pounds  (lb) 

tons  (ton,  2000  lb) 

pound  force  (Ibf) 
pounds  per  sq  in 

(psi) 
pounds  per  sq  ft 
(psf) 

foot-pounds  (ft-lb) 
horsepowers  (hp) 
British  thermal   units 

(Btu) 
British  thermal  units 

(Btu) 

Definitions 

newton  -  force  that  will  give  a  1-kg  mass  an  acceleration  of  1  m/s' 

joule  -  work  done  by  a2force  of  1  N  over  a  displacement  of  1  m 

1  newton  per  sq  m  (N/m  )  =  1  pascal 

1   kilogram  force  (kgf)  =  9.807  N  ^ 

1   gravity  acceleration  (g)  =  9.807  m/s 

1  are  (a)  =  100  m^      p 

1  hectare  (ha)  =  10,000  m 

1  kip  (kip)  =  1000  lb 


0.4536 
907.2 

4.448 
6895 

47.88 


1.356 
746 
1055 


kilowatt-hours  (kWh)    2.93  x  10 


-4 
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FEDERALLY  COORDINATED   PROGRAM  OF  HIGHWAY 
RESEARCH  AND  DEVELOPMENT    (FCP) 


The  Offices  of  Research  and  Development  of  the 
Federal  Highway  Administration  are  responsible 
for  a  broad  program  of  research  with  resources 
including  its  own  staff,  contract  programs,  and  a 
Federal-Aid  program  which  is  conducted  by  or 
through  the  State  highway  departments  and  which 
also  finances  the  National  Cooperative  Highway 
Research  Program  managed  by  the  Transportation 
Research  Board.  The  Federally  Coordinated  Pro- 
gram of  Highway  Research  and  Development 
(FCP)  is  a  carefully  selected  group  of  projects 
aimed  at  urgent,  national  problems,  which  concen- 
trates these  resources  on  these  problems  to  obtain 
timely  solutions.  Virtually  all  of  the  available 
funds  and  staff  resources  are  a  part  of  the  FCP, 
together  with  as  much  of  the  Federal-aid  research 
funds  of  the  States  and  the  NCHRP  resources  as 
the  States  agree  to  devote  to  these  projects.* 


FCP   Category   Descriptions 

1.  Improved    Highway   Design   and   Opera- 
tion for  Safety 

Safety  R&D  addresses  problems  connected  with 
the  responsibilities  of  the  Federal  Highway 
Administration  under  the  Highway  Safety  Act 
and  includes  investigation  of  appropriate  design 
standards,  roadside  hardware,  signing,  and 
physical  and  scientific  data  for  the  formulation 
of  improved  safety  regulations. 

2.  Reduction    of    Traffic    Congestion    and 
Improved  Operational  Efficiency 

Traffic  R&D  is  concerned  with  increasing  the 
operational  efficiency  of  existing  highways  by 
advancing  technology,  by  improving  designs  for 
existing  as  well  as  new  facilities,  and  by  keep- 
ing the  demand-capacity  relationship  in  better 
balance  through  traffic  management  techniques 
such  as  bus  and  carpool  preferential  treatment, 
motorist  information,  and  rerouting  of  traffic. 


*  The  complete  7- volume  official  statement  of  the  FCP  is 
.ivailnble  from  the  National  Technical  Information  Service 
(NTIS),  Springfield,  Virginia  22161  (Order  No.  PB  242057. 
price  .f45  postpaid).  Sinfrle  copies  of  the  introductory 
volume  are  obtainable  without  charge  from  Program 
.\nalysi.s  (HRD-2),  Offices  of  Research  and  Development, 
Federal    Highway    Administr.Ttion,    Washington,    D.C.    20590. 


3.  Environmental  Considerations  in  High- 
way Design,  Location,  Construction,  and 
Operation 

Environmental  R&D  is  directed  toward  identify- 
ing and  evaluating  highway  elements  which 
affect  the  quality  of  the  human  environment. 
The  ultimate  goals  are  reduction  of  adverse  high- 
way and  traffic  impacts,  and  protection  and 
enhancement  of  the  environment. 

4.  Improved  Materials  Utilization  and  Dura- 
bility 

Materials  R&D  is  concerned  with  expanding  the 
knowledge  of  materials  properties  and  technology 
to  fully  utilize  available  naturally  occurring 
materials,  to  develop  extender  or  substitute  ma- 
terials for  materials  in  short  supply,  and  to 
devise  procedures  for  converting  industrial  and 
other  wastes  into  useful  highway  products. 
These  activities  are  all  directed  toward. the  com- 
mon goals  of  lowering  the  cost  of  highway 
construction  and  extending  the  period  of  main- 
tenance-free operation. 

5.  Improved  Design  to  Reduce  Costs,  Extend 
Life  Expectancy,  and  Insure  Structural 
Safety 

Structural  R&D  is  concerned  with  furthering  the 
latest  technological  advances  in  structural  de- 
signs, fabrication  processes,  and  construction 
techniques,  to  provide  safe,  efficient  highways 
at  reasonable  cost. 

6.  Prototype  Development  and  Implementa- 
tion of  Research 

This  category  is  concerned  with  developing  and 
transferring  "research  and  technology  into  prac- 
tice, or,  as.  it  has  been  commonly  identified, 
"technology  transfer." 

7.  Improved  Technology  for  Highway  Main- 
tenance 

Maintenance  R&D  objectives  include  the  develop- 
ment and  application  of  new  technology  to  im- 
prove management,  to  augment  the  utilization 
of  resources,  and  to  increase  operational  efficiency 
and  safety  in  the  maintenance  of  highway 
facilities. 
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